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PREFACE 


This  technical  report  consists  of  publications,  manuscripts,  and 

supplemental  information  related  to  the  scattering  of  light  by  bubbles 

in  liquids  and  in  glass.  The  emphasis  of  the  report  is  on  the  scattering 

by  air  bubbles  in  water  into  the  angular  regions  where  diffraction  makes 

essential  corrections  to  elementary  ray  optics.  For  bubbles  these  regions 

are  backward  scattering,  forward  scattering,  and  critical  angle  scattering. 

Most  of  the  previous  literature  on  light  scattering  is  concerned  with 

drop-like  objects  where  the  refractive  index  of  the  ^catterer  exceeds 

that  of  the  surroundings.  The  present  results  should  be  useful  for  the 

1  2 

optical  characterization  of  microbubbles  at  sea  and  in  water  tunnels  , 
or  for  the  prediction  of  the  optical  properties  of  bubbly  media.  They 
should  also  be  useful  for  1 aser-Doppler  anemometry  with  microbubbles  and 
the  detection  of  microbubbles  in  glass. 

Certain  phenomena  described  here  should  be  present  for  other  spherical 
or  nearly  spherical  scatterers  where  the  refractive  index  is  less  than  that 
of  the  surroundings,  for  example:  air  bubbles  in  ice,  water  inclusions  in 
fused  quartz  (such  as  those  present  in  varieties  of  opal  glass),  ice  spheres 
in  water,  and  superheated  drops  in  a  host  liquid.  The  coarse  structure  in 
the  critical  scattering  region  can  be  present  even  if  the  scatterer  is  not 
spherical .  v 

The  first  paper  summarizes  experiments,  models,  and  results  of  Mie 
computations  for  both  backward  and  critical -angle  scattering.  Theoretical 
results  presented  indicate  that  axial  focussing  enhances  the  scattering  by 
an  air  bubble  in  water  into  the  near  forward  and  near  backward  regions.  This 
enhancement  for  the  forward  region,  (forward  optical  glory),  is  more  signifi- 
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cant  for  depolarized  than  for  polarized  scattering  since  the  latter  is 
dominated  by  ordinary  forward  diffraction.  At  present  the  experimental 
observations  of  axial  focussing  from  single  bubbles  are  limited  to  the 
backward  glory  of  air  bubbles  1'n  a  viscous  polymer  liquid.  (These  are 
described  in  detail  in  the  seventh  paper).  The  first  paper  also  discusses 
the  scattering  efficiency  of  bubbles  and  the  relevance  of  far-field  compu¬ 
tations  to  near-field  observations  of  bubbles. 

The  second  paper  oescribes  the  first  observations  of,  and  model  for, 
the  scattering  in  the  critical  region  which  is  near  83  degrees  to.’'  an  air 
bubble  in  water.  This  paper  is  included  here  for  completeness  though  it 
was  not  supported  by  the  present  ONR  contract.  It  was  supported  in  part 
by  ONR  Contract  N00014-76-C-0527  (R.  E.  Apfel  principal  investigator). 

The  third  paper  describes  an  improved  model  (which  includes  both 
diffraction  and  interference)  for  the  coarse  structure  near  the  critical 
scattering  region. 

The  fourth  paper  is  the  first  to  be  published  which  correctly  gives 
the  results  of  Mie  theory  for  angular  structure  in  the  scattering  by  bubbles. 

This  paper  substantiates  the  model  of  the  coarse  structure  developed  in  the 
third  paper.  The  paper  was  excerpted  from  the  M;  ster's  Degree  Thesis  of 
D.  L.  Kingsbury.  A  supplement  attached  here  gives  examples  which  were  not 
published  for  reasons  of  brevity. 

The  fifth  paper  gives  Mie  theory  and  model  results  for  air  bubbles  in 
fused  silica  glass.  The  paper  is  the  first  to  discuss  the  Brewster  scattering 
angle  for  bubbles  and  the  first  to  give  valid  Mie  results  for  bubbles  over 
the  entire  (0  to  180  degrees)  range  of  scattering  angles. 

The  sixth  paper  describes  computer  codes  used  for  some  of  the  Mie  and 
model  computations.  The  codes  shown  here  have  been  adapted  fov’  use  on  an  HP  1000 
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minicomputer  system.  This  paper  is  excerpted  from  the  Master's  Degree 

Thesis  of  D.  L.  Kingsbury.  For  a  more  complete  description  of  the  Mie 

3 

scattering  algorithm,  consult  the  paper  by  Wiscombe  which  is  the  basis 
of  these  codes. 

The  seventh  paper  describes  the  first  observations  of  backscatteri ng 
from  air  bubbles  in  liquids.  It  also  describes  a  physical -optics  approx¬ 
imation  for  the  backward  axial  focussing  of  scattered  light.  This  model 

shows  that  the  intensity  of  individual  axially  focussed  rays  are  propor- 
3 

tional  to  a  ,  while  the  simply  reflected  scattering  is  proportional  to 
2 

a  where  a  is  the  bubble  radius.  (The  total  backscattering  is  not 

3 

simply  proportional  to  a  due  to  the  interference  of  various  paths; 

see  paper  1).  Polarization  and  quasi -periodic  properties  of  the  glory  are 

discussed. 

The  approximations  and  physical  models  described  here  for  bubbles 
should  also  be  useful  for  certain  cases  in  the  scattering  of  sound  from 

4 

fluid  or  elastic  spheres.  This  has  been  verified  by  Marston  and  Kingsbury 

for  scattering  in  the  critical  region  from  fluid  spheres.  Backward  axial 

focussing  or  "acoustic  glory"  has  also  been  proposed  for  elastic  spheres 
5 

in  water. 

Care  should  be  taken  when  applying  the  results  of  this  report  to 
light  scattering  by  mircobubbles  in  seawater.  Our  computations  assume 
that  the  bubbles  are  spherical  and  that  the  refractive  index  of  the 
scatterer  is  homogeneous.  There  is  a  natural  tendancy  for  the  gas  in 
microbubbles  to  disolve  into  the  surrounding  liquid. ^  There  is  some 
recent  evidence  that  microbubbles  in  seawater  can  be  stabilized  by 
substances  sorbed  onto  their  surface.^  The  presences  of  a  sorbed  film 
or  small  deviations  from  sphericity  may  alter  the  details  of  the 


6 


I\J 


scattering  patterns  described  here,  but  it  should  not  obliterate  the 
existence  of  axial  focusing  and  the  coarse  and  fine  structures. 

Philip  L.  Marston 
Principal  Investigator 
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Paper  No.  1 

Light  scattering  by  bubbles  in  liquids:  Mie  theory,  physical -optics 
approximations,  and  experiments  (P.  L.  Marston,  D.  S.  Langley, 
and  D.  L.  Kingsbury)  to  be  published  in  the  journal  Applied 
Scientific  Research  in  the  proceedinos  of  the  IUTAM  Sumposium 
on  the  Mechanics  and  Physics  of  Bubbles  in  Liquids.  This 
paper  was  presented  at  the  Symposium  (June,  1981). 


Abstract 

Angular  structures  in  the  far-field  scattering  from  bubbles  are 
observed  and  modeled.  Mie  theory  supports  a  model  of  diffraction  and 
interference  near  the  critical  scattering  angle.  A  new  expression  for 
the  angular  spacing  of  fine  structure  is  derived.  Photographs  of  scat¬ 
tering  show  some  of  the  predicted  features.  Application  of  these  struc¬ 
tures  to  bubble  sizing  and  detection  are  summarized  and  the  theoretical 
extinction  coefficient  in  water  is  plotted. 

Mie  computations  for  bubbles  in  water  also  reveal  backward  and  for¬ 
ward  glory  effects.  These  are  partially  manifested  as  cross-polarized 
scattering.  Observed  scattering  from  bubbles  in  ^he  near  bacKward  di¬ 
rection  is  found  to  have  a  strong  cross-polarized  component. 


8 


1 


1.  Introduction 

Equations  for  the  scattering  of  plane  electromagnetic  waves  by  a  di¬ 
electric  sphere  were  given  by  Mie  [1]  in  1908  and  the  resulting  features 
of  the  angular  scattering  pattern  of  drops  are  well  known  [2-4].  The  Mie 
solution,  though  exact,  does  not  give  insight  into  cither  the  scattering 
process,  or  changes  in  the  pattern  resulting  from  changes  in  shape,  re¬ 
fractive  index,  c.  profile  of  the  incident  wave  fronts.  Reviews  of  the 
literature  on  light  scattering  [2-4]  reveal  a  paucity  of  information  about 
the  scattering  pattern  of  bubbles  in  liquids  where  the  refractive  index 

of  the  scatterer  n.  is  less  than  that  of  the  surroundings  n  .  Conse- 

l  o 

quently  we  have  begun  [5-10]  a  systematic  study  of  the  scattering  of  light 
by  gas  bubbles  in  liquids.  Aspects  of  the  study  are:  (1)  the  computation 
of  Mie  scattering;  (2)  the  development  of  simple  physical  models  which  give 
insight;  and  (3)  observations  of  features  in  the  scattering  which  differ 
significantly  from  both  the  scattering  by  drops  and  the  scattering  pre¬ 
dicted  by  geometrical  optics  [11].  In  this  paper  we  summarize  the  main 
features  of  the  scattering  with  an  emphasis  on  the  critical  [5-7]  and  back- 
scatter  [10]  regions.  New  experiments  and  applications  will  be  described. 

These  are  the  first  detailed  observations  of  scattering  by  bubbles. 

Mie's  solution  [1-3]  to  the  problem  of  the  scattering  efficiency  and 
pattern  of  a  dielectric  sphere  is  a  function  of  the  ratio  m  -  n^/nQ  .  It 
is  usually  expressed  as  a  function  of  x  ::  ka  -  2 n.-i/A  where  a  is  the  sphere 
radius,  k  and  A^  are  the  wavenumber  and  wavelength  in  the  outer  dielectric. 
Optical  sources  are  typically  characterized  by  the  wavelength  in  a  vacuum 
we  note  that  k  =  2nno/Av.  For  an  air  filled  bubble,  n.  =  1.00029 
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and  it  is  usually  a  good  approximation  to  take  in  -  1/n  .  Host  of  this 
paper  will  deal  with  scattering  by  spherical  bubbles  with  plane  incident 
waves.  We  will  gain  some  insight  however,  into  which  features  of  the 
scattering  pattern  should  be  sensitive  to  deviations  from  sphericity. 

2.  Scattering  Efficiency 

The  efficiency  factor  Q  is  the  ratio  of  the  total  scattering 

SCd 

2 

cross  section  to  the  geometric  cross  section  ira  .  When  n.  is  real 

valued,  Q  also  gives  tfie  extinction  efficient.'/  [2-4,12].  Figure  1 
sea 

compares  Qrca  for  drops  of  water  in  air  with  m  -  4/3  with  that  for  a 

bubble  in  water  with  rn  -  3/4.  The  computations  were  performed  using  a 

slightly  modified  version  [9]  of  Wiscombe's  MIEVO  Mie  scattering  algorithm 

[13].  A  table  of  Q  for  bubbles  exists  [14]  which  is  consistent  with 

Fig.  1.  For  a  fixed  Ay,  and  a  given  value  of  x  ,  the  drop's  radius  is 

larger  than  the  bubble's  radius  by  a  factor  of  4/3.  For  the  case  of 

light  from  a  He-Ne  laser,  \  =  0.6328  pm;  for  bubbles  in  water,  x  of  10, 

100,  1000,  and  10  000  give  a  of  0.75,  7.55,  75.5,  and  755  pm,  respectively. 

The  salient  feature  of  Fig.  1  is  that  for  drops  Q  exhibits  a 

sea 

fine  "ripple"  structure  [3,12]  but  that  our  calculations  of  Q  for 

sea 

bubbles  do  not  reveal  a  ripple  structure.  For  both  drops  and  bubbles,  Q 

sea 

has  a  broad  undulation  with  a  quasi-period  Ax  ~  ir/|m  -  lj.  This  approxi¬ 
mation,  which  has  been  derived  from  the  theory  of  "anomalous  diffraction" 
[2,3],  appears  to  be  useful  for  both  bubbles  and  drops.  The  ripple  struc¬ 
ture  present  for  drops  is  due  to  optical  resonances  [12]  which  are  attributed, 
ir,  part,  to  internal  surface  waves.  The  absence  of  such  structure  for 
bubbles  is  probably  because  m<l  does  not  favor  the  entrapment  of  internal 
surface  waves. 


10 


3 


3.  Critical  and  Brewster  Angle  Scattering 


For  scattering  by  drops,  diffraction  is  important  for  the  description 
of  the  forward,  backward ,  and  rainbow  regions  [2,3].  For  bubbles,  there  is 
no  longer  a  rainbcw;  however,  a  new  region  appears  known  as  the  critical 
scattering  region  [5],  Diffraction  is  important  in  this  region  because  of 
an  abrupt  change  in  the  amplitude  of  the  reflected  wave  as  the  local  angle 
of  incidence  0  changes  from  6  <  0c  for  small  impact  parameters  to  0  >  0c 
for  large  ones.  Here  0c  =  arcsin(in'^)  which  is  the  critical  angle  for  a 
plane  surface.  Figure  2  illustrates  several  ray  paths  which  lead  to  a 
scattering  angle  <J>  (the  deviation  from  the  forward  direction)  of  50°. 


The  number  on  the  left  specifies  the  number  of  internal  chords 
denote  the  angle  of  incidence  and  refraction  of  the  pth  ray. 


Pi  0  and  Pp 
The  reflected 


ray  (which  has  p  -  0  )  has  a  scattering  angle  <f>  =  n  -  26_.  For  an  air 

u 

bubble  in  water  with  m  =  3/4,  the  critical  scattering  angle  is  <t  =  it  -  20  = 

c  c 

82.82°.  Geometric  optics  [5,  6,  11]  predicts  that  | d I . /dd> |  ■*  «=  as  ch  ap- 

J 

proachcs  $  from  above  4  .  Here  I.  is  the  normalized  scattered  intensity 
c  c  j 

defined  as  follows:  the  actual  j-polarized  intensity  at  a  distance  R.  >>  a 


from  the  bubble's  center  is  the  incident  j-polarized  intensity  multiplied  by 
2 

I.  (a/R)  / 4.  For  the  electric  vector  perpendicular  to  the  scattering  plane, 

J 

j  =  1,  for  the  parallel  case,  j  =  2.  This  normalization  is  appropriate  for 


bubbles  since  geometric  optics  predicts  tnat  if  the  intensity  of  the  p  =  0 

reflection  could  be  considered  by  itself  I .  (<j>  <  <j>  )  =  1  due  to  total  reflection. 

J  c 

Our  observations  [5],  model  [5,  6],  and  Mie  computations  [7-9]  demonstrate 

that  instead  of  a  divergence  of  |dl./d<j>|  ,  the  rise  in  I.  is  spread  out  over 

J  J 

the  region  | -  4>c |  -  SI  where  Cl,  =  arcsin  [0.8  (l-m2)_Js(X  /a)!*]  *  144x_iideg 
for  water.  The  model  makes  use  of  a  physical  optics  approximation  which 
is  to  (a)  use  ray  optics 
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along  with  the  reflection  coefficients  of  a  plane  surface  to  compute  the 
complex  amplitudes  of  virtual  waves  which  simulate  the  scattering,  and  (b) 
use  Fraunhofer's  approximation  to  compute  the  diffraction  of  the  virtual 
wave  to  the  far  field  where  R  >>  xa.  Step  (b)  yields  the  improvements 
over  the  geometrical  predictions  [11].  This  procedure  is  analogous  to 
Airy's  model  of  diffraction  neat'  the  rainbow  [2]  .  Additional  approximations 
used  in  the  model  are:  only  the  two  most  intense  virtual  waves  for  <t  -  <p 
are  included  (p  =  0  and  1);  and  the  divergence  of  the  derivative  of  the 
reflectivity  is  simulated  by  truncating  the  p  =  0  reflection  when  Oq  <  0C  • 

To  compare  the  model  results  with  Mie  theory,*  we  consider  examples  of 
bubbles  in  water  which  complement  those  previously  published  [7].  Figure  3 
shows  both  theories  for  x  =  100  and  j  =  1.  It  shows  that  there  is  a  coarse 
structure  in  the  Mie  scattering  for  $  <  which  has  a  quasi -period  <.  SI  ~  14.4". 
This  structure  is  described  by  the  model  exceot  in  the  near  forward  direction 
(0  £  20°). 

The  Mie  result  also  has  a  superposed  fine  structure  where  the  magnitude 

of  the  quasi-period  <  \Ja  rad.  =  360°x~^  .  Near  6  ,  this  fine  structure 

0  c 

arises  primarily  from  the  interference  of  the  p  =  0  wave  with  the  wave  due 
to  the  p  -  21  ray  in  Fig.  2.  Its  quasi-period  .nay  be  estimated  from  the 
lateral  separation  bQ  +  b?,  of  virtual  point  sources  which  would  simulate  the 
scattering  at  a  given  .  As  shown  in  f  ig.  2,  this  separation  is  the  sum  of 
the  impact  parameters  for  these  rays.  Standard  relations  for  the  far- 
f*  i  el  cl  interference  applied  to  these  sources  gives: 

nf(<|>)  =  arcsinUo/(bQ  +  , ) ]  (1) 


*The  conversion  from  the  Mie  amplitudes  3-  [2,:i  13]  to  the 

UlSjl/x)2. 


I  .  is 
J 
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i  mu  i  i  i  i  ii  ngm 


where  b^/a  =  sinO^  »  °q  =  U  ~  <}>  =  Tr  +  2  ( 0  ^ ,  -  2p^ , )  an1!  from  Snell's 

law,  msinpp  =  sinO^  .  Lq.  (1)  is  only  an  approximation  because  it  fails 
to  include  other  (such  as  p  =  1)  virtual  waves  and  it  does  not  completely 
account  for  the  longitudinal  spacing  of  virtual  sources  and  for  how  those 
sources  vary  with  <|>.  With  <{.■  =  <J>  and  ni  =  3/4,  b^  +  b^,  ~  a/0.825  and 
fir  -  arcsin  (5.18  x  ' ) •  This  gives  2.97°  at  x  •-  100  which  is  in  reasonable 
agreement  wi th  fig.  3 . 

Figure  4  shows  the  Mie  and  phys ica I -optics  results  for  x  -  10  000  and 
j  =  1.  The  fine  structure  quasi-period  is  greatly  reduced;  fiq.  (1)  gives 
8fUc)  ~  0.0297°  which  agrees  with  the  Mie  result  of  0.03’.  The  amplitudes 
of  the  coarse  undulations  (n  ~  1.4°)  decrease  slightly  with  decreasing  d  until 
4)  ~  75°.  A  graph  of  1^  for  |>  •:  75°  shows  coarse  undulations  increasing 
in  amplitude  with  decreasing  <;■.  The  explanation  is  that  for  x  =  10  000, 
with  <f>  ^  75°  the  coarse  structure  is  primarily  due  to  diffraction  of  the 
p  =  0  wave  [5]  but  for  <}>  ^  75°  it  is  primarily  from  the  interference  of  the 
p  =  0  wave  with  p  >  0  waves  [6].  Plots  [7]  of  for  this  x  also  show 
this,  transition  and  that  fine  structure  is  significantly  weaker  when  j  =  2. 

Figure  3  does  not  show  the  transition  because  the  diffraction  region  (•}>  -  <|>  ,<  Cl) 

overlaps  the  region  where  interference  with  the  p  =  1  wave  is  significant. 

As  Qq  approaches  the  Brewster  angle,  0^  -  arctan(m),  the  reflectivity 
of  the  p  =  0,  j  =  2  polarized  ray  vanishes  [2].  Consequently,  according  to 
geometric  optics,  there  is  no  contribution  to  from  this  ray  at  the  Brewster 

scattering  angle  f ..  =  it  -  2up  -  2  arctan(m  ^ ) .  Because  the  reflectivity  varies 

t)  b 

slowly  near  <pR,  diffraction  is  less  important  than  near  >}>  .  A  coarse  minimum 
is  evident  in  the  Mie  near  .p  for  bubbles  in  glass  [8]  and  for  bubbles 
in  water  where  <t>g  ~  105.3°.  When  x  drops  below  5,  its  location  shifts  toward 
90°;  the  scattering  pattern  approaches  of  o  dipole  radiator  predicted  by 

Rayleigh  scattering  theory  [2,  3,  8]. 
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4.  Observations  of  Critical  Angle  Scattering 

Previous  observations  [5]  were  limited  to  the  4>  region  where  the 
coarse  undulations  decreased  with  decreasing  -ji  .  We  have  a  new  apparatus 
similar  to  that  in  [5]  except  that  it  permits  observations  with  4)  down  to 
70°.  As  in  [5],  bubbles  were  attached  to  a  vertical  needle  in  distilled 
water.  They  were  illuminated  with  a  plane  wave  from  a  He  -  Tie  laser  with 
-  632.8  mil.  They  were  photographed  via  a  window  with  a  camera  focused 
on  infinity  to  place  the  film  in  the  far  field.  Because  of  a  technical  problem 
we  have  not  obtained  precise  direct  measurements  of  ■;  in  the  (horizontal) 
scattering  plane  along  with  the  photographs,  however,  we  have  demonstrated 
that  )  from  Eq.  (1)  gives  an  a  which  is  consistent  with  the  observed 

and  modeled  coarse  structure.  Thus  the  values  quoted  in  F'g.  5  were 
determined  from  the  coarse  and  fine  structures  and  not.  from  direct  observations 
Figure  5  serves  more  to  illustrate  phenomena  described  in  Sec.  3  than  to  rig¬ 
orously  test  the  models  with  real  bubbles.  The  angle  scale  is  for  ?.  t}> 

change  in  water.  The  horizontal  displacement  is  not  exactly  linear  in  if>  due 
to  refractive  corrections  at  the  water-window-aii  interfaces.  In  fig.  5(a), 
the  broad  bright  region  on  the  left  is  the  first  coarse  maximum  with  <},  <  . 

On  the  right,  the  coarse  undulations  increase  in  amplitude  as  <J>  decreases, 
which  is  the  modeled  behavior.  Figure  5(b)  is  for  a  smaller  bubble  (verified 
with  direct  observations).  The  coarse  undulations  are  spread  out  and  the  fine 
structure  is  clearly  visible  with  -■  0.07°. 

5.  Backward  and  Forward  Glory 

Backscattering  from  drops  is  known  to  be  enhanced  by  the  axial  focusing 
of  certain  rays  [2,  5].  A  complete  description  of  this  "glory"  for  water  drops 
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is  complicated  due  to  the  ner-ssity  of  including  surface  waves  [15].  We 

have  observed  and  modeled  bur.kscatlerincj  from  air  bubbles  [10]  and  find  that 

some  aspects  can  be  explained  with  a  physical -optics  approximation  that  does 

not  include  surface  waves,  for  the  following  discussion,  it  is  convenient 

to  define  y  =  it  -  <J>;  y  -  0  corresponds  to  exact  backscattering . 

Computed  Mie  intensities  show  that  backscattering  from  bubbles  may  be 

significantly  larger  than  that  expected  from  reflection  of  the  p  =  0  ray. 

For  y  =  0,  the  plane  of  scattering  is  no  longer  defined  and  =  I. 

! igu'c  6  shows  several  broad  peaks  in  I  which  arj  significantly  larger  than 

the  normalized  intensity  of  the  p  =  0  reflection  which,  by  itself,  is  (m  -l)'/ 

2 

(m  +  1 )  =  0.02.  Other  axial  (e.g.  p  2)  rays  are  too  weak  to  explain  the 

magni tude  of  I . 

For  polarized  incident  light,  scattering  near  y  -  0  may  be  described  in 

part  with  (non-cross}  polarized  1^  and  cross-polarized  1  ^  normalized  in- 
(2) 

tensities.  For  r  ,  the  plane  of  polarization  is  rotated  by  90°  from  incident 

polarization;  symmetry  gives  =  0)  -  0  wiiile  l^(y  -  0)  =  I.(y  =  0). 

in 

For  y  ^  0,  the  I  depend  an  both  y  and  the  angle  C  which  the  scattering 
plane  make,  with  the  incident  electric  vector.*  Figure  7(a)  illustrates  the 
computed  scattering.  Geometric  models  of  I^(y  ~  0)  predict  that  the  con¬ 
tributions  of  off-axis  (p  -  3  and  4)  rays  diverge  as  y  -►  0  because  of  a  factor 
which  accounts  for  focusing.  This  divergence  is  present  in  previous  ray  models 
[11],  but  was  not  discussed;  it  is  also  present  for  drops  [2]  for  certain  m  . 
Because  of  this  divergence,  diffraction  provides  an  essential  correction  to 
ray  optics. 


*The  conversions  from  the  complex  Mie  amplitudes  S.  to  the  I  ^  (r'<  1 0° , 

.2,2  ■  .(2)  .  .2,  2 


are  I 


(1 


5  -  S  00^2',  i  / x  and  I 


S  sin2f, !  /x  where  S' 


S-j  1  5^-  Evidently  the  dependence  of  Eg.  (2)  is  of  the  correct  type. 
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The  scattering  due  to  the  p  -  3  ray  in  Fig.  /(a)  (insert)  is  a  back- 
facing  toroidal  wave  front  which  appears  to  originate  at  a  virtual  ring¬ 
like  source  known  as  the  focal  circle  in  the  analogous  [2]  p  ?.  scattering 
from  drops.  This  source  at  point  F  is  ringlike  because  the  figure  may  be 
rotated  around  the  optic  axis  through  the  center  C.  For  largo  x,  the 
stationary-phase  approximation  of  far-field  diffraction  integrals  gives  the 
following  proportionalities  for  otf-axis  (p  >2)  contributions  to  the  1^ 
taken  separately: 

\  9 

I(p)  cx  xt(c1  +  c2)  J0(u)  4  (c]  -  c2)  d,,(u)  cos2;JL  I 


l{pj  xUc1  -  c2)  J2(j)  sinZ;.]? 


where  the  c.  are  combined  Fresnel  reflectivities  and  transmissivities  of 


the  j-polarized  fields,  u  =  kb^siny  ,  and  bp 


is  the  impact  parameter 


of  the  exactly  backscettered  pth  ray.  Because  l|  j  <•>  >.  ,  the  unnnrmal  i zed 

intensity  is  proportional  to  k4  while  that  for  the  p  =  0  ray  is  proportional 

2 

to  ,/  .  Consequently,  off-oxis  waves  are  dominant  when  <i  is  sufficiently 

large.  It  is  evident  from  Fig.  b  and  7(a)  that  the  modeling  of  I'1  ^  may  re¬ 
quire  summation  of  electric  fields  from  several  virtual  sources.  (In  physical 
optics,  intensities  do  not  add  but  fields  do.)  Interference  of  various  ring¬ 
like  and  axial  virtual  sources  depends  on  x  and  could  lead  to  the  broad  struc¬ 
tures  in  Fig.  6.  Detailed  computation  of  ij^j  (y  -■  0)  give  magnitudes  suf¬ 
ficient  to  explain  the  Hie  scattering.  Due  to  the  Bessel  functions  in  Fq.  (2), 
the  actual  scattering  should  be  peaked  at  or  near  y  0. 

We  have  observed  backscattering  from  single  air  bubbles  in  a  dimethyl- 
si  loxane-polymcr  liquid  which  were  nearly  immobilized  by  viscosity.  In  these 
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experiments ,  in'1  =  1.403,  =  632.8  nm,  a  -  0. 3  -  0.8  mm  and-  x  -  4000  - 

n  000.  The  far-field  cross-polarised  intensity  had  a  dependence  on  y  and 
similar  to  that  predicted  by  Eq.  (2b)  with  p  -  3.  Equation  (2b)  predicts 
that  when  both  siny  ^  y  and  u  :•>  1,  the  minima  in  1^  should  be  spaced 
by  Ay  ~  it  (kb^)'1  radians  where  b'  -  0.44/a.  The  prediction  correctly 
described  the  observed  i  -  2  scattering  which  is  apparently  dominate^  here  by 
the  p  =  3  virtual  source.  The  focal  circle  was  viewed  by  focusinq  the  camera  on 
the  bubble.  The£=  1  scattering  was  not  dominated  by  a  single  class  of  rays. 
Cross-polarized  scattering  from  a  polydi spersicn  of  bubbles  in  water  ha s  also 
been  seen. 

Figure  7(b)  shows  an  enhancement  of  near-forward  cross-polarized  scattering 
due  to  axial  focusing.  Focal  circles  due  to  the  p  =  ?  and  3  rays  largely 
contribute  to  the  C  ~  2  scattering.  The  forward  y  -  1  scattering  is  domi¬ 
nated  by  ordinary  diffraction  [2]  when  x  is  large.  Hie  theory  gives  I^~  10^ 
when  tj.  =  0  for  the  x  and  m  of  Fig.  7.  Forward  optical  nlory  has  also  been 
displayed  in  Min  scattering  by  water  drops  [16].  Its  description  is  complicated 
by  surface  waves. 

5.  Discussion 

Figures  3-7  demonstrate  that  for  several  angular  regions,  the  intensity 
exceeds  that  of  geometric  scattering  from  a  perfectly  reflecting  sphere  of  the 
same  a  (for  which  K  =  1).  This  information  should  be  useful  in  optical  de¬ 
vices  which  size  or  detect  bubbles  [7,  17,  18].  Due  to  diffraction  near  4>  , 
it  is  preferable  to  detect  scattering  with  .(>  <  -  u  than  to  observe  it 

with  <J)  =  90°  which  is  the  usual  practice  [l.v]  for  bubbles  in  water.  For  an 
unpolarized  source,  the  normalized  total  scattering  is  i ! -j  1  I ^ ) / ? . 
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Though  this  paper  has  emphasized  the  far-field  |  dependence,  the 
results  are  also  applicable  to  the  imaging  of  bubbles.  For  example,  if  it 
desired  to  resolve  the  virtual  sources  of  the  p  =  0  and  1  rays  in  Fig.  2, 
the  aperture  of  the  imaging  system  should  have  an  angular  width  >  Q  .  It 

*\j 

is  easier  to  resolve  the  p  -  0  and  2'  virtual  sources  since  the  angular 
width  requirements  are  reduced  to  Resolution  of  identifiable  virtual 

sources  reveals  the  size  of  a  bubble. 

Fine  structure  in  the  far-field  scattering,  such  as  that  shown  in  Fig.  5(b), 
arises  from  the  interference  of  widely  spaced  rays.  C jnsequemiy  the  positions 
of  the  maxima  will  be  quite  sensitive  to  changes  in  the  bubble's  shape.  That 
part  of  the  coarse  structure  due  to  the  interference  o\  p  -  0  and  1  rays 
will  also  depend  on  the  shape  but  more  weakly.  As  is  the  case  for  drops  [19], 
these  shape  dependences  may  be  useful  for  detecting  mechanical  resonances. 

Scattered  intensities  for  bubbles  in  water  purported  to  be  frori  Mie  theory 
have  beer,  used  in  a  study  of  cavitation  nuclei  [1/].  Comparison  of  our  Hie 
results  and  model  of  the  coarse  structure  (which  are  consistent  for  the  a  in 
question)  with  those  in  [17],  show  that  the  latter  err  significantly.  Plots 
in  [17]  predict  coarse  maxima  at  $  =  100°  for  a  =  5  pm  and  <p  -  120"'  for 
a  "  7.5  pm  which  are  not  present  in  our  Mie  computations.  We  cannot  find  any 
physical  justification  for  coaise  maxima  at  these  for  x  near  100. 

We  are  grateful  to  W.  J.  Wiscombe  for  providing  the  initial  comouter 
program  from  which  the  program  used  here  was  derived.  This  research  was 
supported  in  part  by  the  Washington  St  te  University  Research  and  Arts 
Committee  and  by  the  Office  of  Naval  Research,  P.  L.  Mars  ton  is  an  Alfred 
P.  Sloan  Research  Fellow. 
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Figure  Captions 


Fig.  1  Scattering  efficiency  from  Mie  theory  as  a  function  of  the  size 
parameter  x  .  The.  solid  curve  is  for  drops  and  the  dashed  curve  is  for 
bubbles . 

Fig.  2  Ray  paths  for  a  bubole  with  in  =  3/4  .  The  number  of  internal 
chords  p  is  given  a  prime  if  the  ray  enters  the  bubble  below  the  center¬ 
line. 

Fig.  3  Logarithm  (base  10)  of  the  normalized  scattered  intensity  pre¬ 
dicted  by  Mie  theory  (solid  curve)  for  x  =  100.  The  dashed  curve  is  the 
physical -optics  approximation  [6]  of  the  coarse  structure. 

Fig.  4  Like  Figure  3  but  with  a  linear  scale  and  x  -  10  000. 

Fig.  5  Photographs  of  far-field  scattering  for  (a)  the  bubble  radius 

a  ~  400  urn  (x  '  6360,  •  1.3°)  and  j  =  2;  (b)  330  pm  (;;  ■  4370,  3.2 

and  j  -  1 .  The  critical  scattering  angle  ■!■  is  near  the  left  edge  of  the 
photograph;  <}>  decreases  from  left  to  right.  Coarse  structure  is  manifest 
as  broad  vertical  bands.  The  rings  are  artifacts. 

Fig.  6  Normalized  Mie  backscatteri ng  from  a  spherical  air  bubble  in  water 
(m  -  3/4)  plotted  as  a  function  of  the  size  parameter  x;  for  green  light, 
a  ••  200  pm. 

Fig.  1  Normalized  Mie  neac-backscottering  (a)  and  near-forward  scattering 
(b)  for  m  =  3/4  and  x  =  3040.  For  clarity,  the  computed  r  '  was  multiplied 
by  two  in  (a)  before  plotting.  The  inserts  illustrate  soiae  of  the  significant 
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The  intensity  of  light  scattered  hv  nit  nir  buhhle  m  water  is  predicted  by  (tie  geometric-optics 
calculation  >f  M!l5fn  to  tune  a  d'voic'tit  .inn.Vr  derivative  as  the  i  ritual  suOtermK  unglu 

4>r  is  npptoiu  hed  l'.H' 'i  ts  of  dilVrartion  ri  the  angular  region  m  ar  <!>  are  described  here.  The 
Fraunhofer  diffraction  for  scatt'  ring  an.;', "a  estimated  using  a  «t*ii,«lilied  physical-optica 

approximation  A  ringing  and  itec.iv  of  mi"  iar  field  eiteivutv  i  >  nredieteU  tiiat  ;s  fainmUv  :  inula.'  to 
the  near-field  diffraction  at  a  straight  ‘  lire  Oh--  v  ition  of  nalluneter  radius  hui'hlc.  ■  n  water  with 
collimated  mnnochnmn’ir  iilumm.it ion  eimtin.i  Ue-  “xisteia  e  of  this  nitiru'.:  wl  .  h  has  a  qua  a  pot'ad 
25  mrad.  The  diflVaetinn  raleulati'ai  gives  an  approximate  description  of  tin  relative  <f>  of  the 
observed  muxima  and  minima  Fungi's  witli  a  lower  omura.st  and  spacing  ■  It  ’1  mrnd  were  also 
observed;  they  appear  to  he  caused  hv  t ! i  ii.i<  rfen  ia  e  of  lays  with  distinct  paths.  Implications  for 
the  critical  angle  muttering  of  white  light  are  discussed 


INTRODUCTION 

The  purpose  of  tilts  paper  is  to  prevail  on  a  tu  tiw  nts 
and  tin  approximate  I  heorel  ir.il  dc>i  ripi -■  n  "l  <i  1 1 1 1-. turn 
phenomena  oliservcd  when  light  i .  m  . uteri  d  near  the  cnlu  ol 
angle  Irom  a  laihlile  wit li  an  in'erior  t'Trui  -live  index  li  than 
that  of  the  outer  mi  die.  The  opt  e  a  I  wave!,  npth  is  assumed 
In  he  very  much  less  than  tin  i adiiis  id  t In-  ue.ii Iv  pin  m  . il 
inthhlr.  Approximations  given  here  i  mn  lining  tin  uehavtor 


ol  the  reflect  ion  eoetl  ii  lent  and  the  extent  of  i  p'.ei  ;'tn;;  Wive 
front  x  nil  I  an  ■  xpri  .si,m  for  I  lie  seat  tel  it  ig  in  t  .  ie  tar  held  that 
i-.  similar  in  lorm  to  the  liesiiel  (near  field  I  hit  t  ’'a  el  ion  hv  a 
t  r  light  eti.'.e  <  tli  .et  v  at  e  ms  ri  ported  here  ol  the  1 1 1  <  'me 
In.  a  i  mi  hi  dr  m  wall  r  siippm  I  some  a  .peels  ol  |  pis  .mu  ill  I  ad 
llieotv.  The  Ireatmen'  may  lie  readilv  extended  to  include 
seatlerim:  front  ot  her  eorveii  stirl.u  r  such  as  that  ol  a  cylin 
di  r.  |in  iv.i  fed  that  the  m  uu  p!e  i  ail  it  m  'em  va'ure  are  much 
I  eat  i  I  '  hull  t  lie  wav  nlelll't  h  ol  I  Igllt . 

<  i  '.Cl  ( ) | > t  ii  li  Soi  let  >  ol  Aniern  a  1.05 
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HG.  1  1  o!  il  nflorhon  *n  tho 

plant*  ol  sr.;iN< 'M(k)  nl  .1 1  ly  with  .m 
impact  pm mtiotoi  h  l.iKjor  ihmt 
tho  cntu.ni  impact  paiamotor  t\ 
Only  onr  iju.i  H.mt  of  Hu*  cross 
suction  nl  tho  tHiholo  r..  shewn 


Ha  vis1  calculated  the  utnpilur  (list  1 1  In  it  inn  nl  the  intensify 
ol  liftht  scattered  hv  an  air  bubble  in  ss  an  I  rum  a  i  nil  mint  eel 
incident  lieam.  The  rah  mint  inn  ss  as  ha  ei  I  nh  ia  mnei  rn ■  < -pi  ics 
and  the  (''rcsnel  rellect  mm an  I  In  a  nt nt  In  In  al  plane  inlei 
fact's;  fa  It'tl  1  a  t  mns  nl  I  his  1\  n  e  eon  !inin  an  a|i(>rn  •  iiii.m  mn 
that  is  nsrlnl  only  when  the  railin'-  nl  the  huhhle  a  i  very 
much  greater  than  the  ssusehn;;th  in  I  lie  miter  ineil.a  I  .el 
t  he  scat  t cm n;  all|t le  (/»  nl  a  i  a  .  denote  1 1  le  des  ni  t  mi ,  m;: le  nl 
the  scattered  fa;  t  ruin  I  lie  dire.  I  mo  nl  lie  ,  m  t  >  i  •  >u  put  a  1 1  e  I 
beam  ot  l  ays  l  see  I  II.  A  sir;ml  e  ant  tea l  ur  •  "I  I  i.n  ■  .  cal 
cutatioli  is  1ii.it  (here  should  I"  an  iloupt  din  ma  e  in  (lie 
sent  ter  ed  intensity  when  I  lie  'a  at  ter ii’".  ali;;|e  .  ■  emeil,  Hi. 
where 


ip,  "■  rr  -  L’d  .  Ill 

sinft,  -■  ri  C.!l 

and  n  =  (n„/n,  I  >  I ,  n„  and  n,  lie  me.,  i  >  .|iei  a  i\  eh. ,  the  re 
fractive  indices  ol  the  miter  and  the  nun  i  media. 

Figure  I  illustrates  the  physical  origin  nl  this  predieterl 
decrease.  Consider  how  'he  an;;le  ol  uu  idem  e  at  'he  |>uhhh 's 
surface  9  varies  a-  the  impai  t  parameter  nl  the  rav  >  un  reuse 
Kavs  with  an  impact  purainefl'i  in  eyeess  nl  a  entn  al  value  h, 
have  II  >  II,  where  d, ,  t;iven  liv  l'i|  I  ’.'I,  is  l  he  critical  a  little  tor 
a  pinin'  interlace.  At  r.iriiuij:  to  the  fteometrie  optics  ap 
proximal  inn,  rats  with  h  •  />,  n  c  ')  will  hi'  loiallv 
rejected  Hays  with  />  s  />.  will  In  partialis  lelh'clul  ac 
cn,  din;!  to  the  matt  ml  title  ol  I  lie  I  re-net  teller!  mil  coellicienls 
and  I  his  part  ial  relied  ion  leads  to  tin  uhriii’i  .  1 1  •  ■  •  e  t  ■ : « -e  in  i  la' 
scattered  intensity  for  <;i  *  •  I  i.c.  i ealculati  d  . it  t >  ruip 
intensity  also  includes  the  mult i pie  t el h  ,  i inn  and  tr  iiiMin  . 
smn  ol  the  rays  with  impact  paiametei-, /•  ■  /i,  Tin  ruvsthat 
petitd rate  t  lie  sphere  and  aie  part  i. ills  reflected  emit  ■  ihute  to 
scat  It  lit;  tor  llii  entire  ranee  n|  however,  i  hey  do  not  oh 
sett  re  I  lie  i  list  im  t  r!ec  reuse  in  the  si  at  let  "d  mien  at  v  I  nr  </>  * 
•A,  This  is  evident  'mill  I  nit  ni;;  pi  ede  1  e>n  lee;  ISntKel. 
II  lur  the  intensity  ■  it  an  an  hnl 'hie  in  \  .  i  tho  ha-  H  ~ 
HTH° ,  it  i  .  a  1  -  n  ey  nit  nt  in  the  umI1  pm  ■  1  ite.;  |  a  '  -  in  t  inie  !nr 
huhbles  in  other  llipnds  ll'e;.  ■  o|  Hi  i  'Jl  The  unerispv 
Clint  t  ihllt  ion  III  I  he  nillll  Iply  idle'  ted  and  ret  i  .III!  d  rax  s  Will 
lie  omitted  in  the  treatment  presented  here  suu  e  w-“  are  <  on 
side  r  in;;  olllv  the  el  teds  ol  dd  I  rad  mo  on  the  sea  l  lend  ill 
tensity  for  small  |<,i  -  j 

Tht*  calculations  ot  Davis  prcdn  t  t  hat  d/,  :d  .  is  di scout  in 
tious  at  Hi,  where /,  i '/>l  is  I  lie  tnl  .i  .itoied  mil  n-ity  IM 
tract  ion  ot  h;;ht  will  not  per  not  ill ,  'h;.  in  he  d 1  a  mil  imimis. 
The  method  presented  here  is  a  .nupu  i,n  )  le  d  Ini  t-.l  i.nal  nij; 
the  effects  ul  dill  rail  mu  n  the  ri';;im>  win  n  n  .  |  i .  .mall 

['  irst,  the  approximate  shape  nl  l  lie  otoal  w  ay  e  I  mnt  i ,  i  a  I 

I  Dpi  Sac.  Am.,  \  nl.  till.  Nn  '.I.  Sentemher  i'i,‘'l 


ciliated  pei am  l  rn  ally  lm  those  raVs  with  h  ">  b,  Hays  wit  h 
b  <  h,  w ill  he  mini  led  i  ee  Si  i  ID.  The  virtual  wave  trout 
is  t  hen  inodil  a  d  hy  a  ph.e  e  i  m  reel  ion  due  to  the  phase  shit i 
ot  the  Uiiliindid.o  I  n  nel  ainphlude  reflection  inellli'lent- 
tm  h  •  />,  I  he  in"  iile'd  way  e  Irmit  ih  scrihcs  the  apparent 
■  mm  i  viewed  In  m  "Ii  e'yrr  lai  Irnl'.i  Ihe  huhhle.  (lie 
Fr.i'.mlinli  i  diUi.ni  ai  nl  tlii  wave  pives  the  ih  -ired  ap 
prnxiiu.it  inn-  tm  lie  ..  oiiiim;  atiiphlude  and  intensity 

The  w.c  e  I rmo  i  i  .  .. lined  I o extend  over  ahull  plane  and 
is  not  limited  by  t  lie  limti  u/e  ot  t  lie  bubble.  The  !  real  men t 
resembles  Airy-  oilr.u  turn  correct  ion  to  I  iescartt  ■>  ’  poo 
luetrii  oplus  lre.ilmeiil  ol  the  rainbow  caustic  1  ’  Kveti 
thouph  Airs  ‘ .  approMiu.il  uni  .  <  insiders  the  virtual  wave  I  root 
to  Ih  n!  ml  .a iti  ey  ti  nt  i  anil  a  di  up  must  lie  finite),  it  is  uselul 
tm  the  do  mm  mt  polar  i,' a  I  imi  w  hen  the  radius  nl  t  he  sphere 
hi  red  I  tie  wni  h  1 1  o  !i  i i v  i  la.  tor  nt  approx i ma t el s  Mill.'1  ' 
III  mi  a  i pa l ml  t '  a  .i  m on  si  1 1. mi  criterion  max'  delineate 
tlie  ranee  ol  n  -I  tin  approx  i  mat  tons  presented  here 

tor  dill  i  art  mu  ‘e.  .,  t,  hide  in  tlm  rrpinti  nt  critical  anple 

M  ea  ni  emc  i  '  s  '  .  o  I  rn  i  nj',  i  d  ii  lull'  I  lull  I  I  iv  an  air 

lull  ii  ile  in  '\  .0  ■■■  u  n  '  c  I  I  is  t  •  I  •  ( I'l;’.  1  ol  He!  7 1  I  lies 

-III  •  I  a  '0  l  .1 1  e  1  ’  pi  c. !  I,  s  ;.  I  nl  ,  e  '  I  I'l  e  Opt  il'S  III  ,’l  rapid  de 

ere. I  in  /.  :  •!  tie  i-e  ,  0  I ,  lie  hi;  lit  and  t  he  1 1  III  It  ml 

a  1 1  i ;  u  1. 1 1  ii-  >1' n  inn  nl  t  hi  an  e.iu  i  inent ..  make  t  I'.cni  iiu-uita  hie 
tor  i  mi 1 1  e  i  .1 ,  u  d  "i  letem  t  In  orv.  Tlie  nlisersal  i  hi 

ss  0  h  la  .e|  ilh  u m  n  11  mu  upmt  ed  Ill'll'  111  St-  .  I  V  Old  ii  ale  1 1  Ml 

the  amph  lilttu  . .  :u ms  rises  a  lair  description  ol  the 

lelat  ive  pie.o  u  o  !'  •  fn.i  •  ima  ami  tiitmma  ol  t  he  mt en  at  v 
lit  1  he  re  ;imi  ■  •!  i  1 1: ..  ,.l  an;:  ie  atterinp 


I.  CUHVATUmL  of  ihe  virtual  wave  front 

I'll'  •)  j  <  '*  dm  fiat'  a  ( 'art  i -mil  eonrillllale  system  u.-ed  ill 
l  lie  ih"  '  1 1 pie  ni  id  the  s  u  ni.il  w  is  r  Irmit  The  nruy.n  is  cl  him  m 
to  in  I  he  pun  I  w  In  i  e  t  lie  i  as  w  it  1 1  h  h,  contacts  ( lie  sphere 
unit  I  lie  I  axis  I  •  pa '  ilh  i  to  the  men  lent  raV.  I  'oil  its  (  ■'  .S'  I  nil 
the  virtual  .'..ivi  t.  ail  an  obtained  by  extendi  t|;  tin  ivo 
met  r n  alls  dm  ,a  me  .ed  .  at  tered  ras  ha ek Ward  a  delate  h 
t mm  I  In  |  mint  i  !  in  I  ii  i  •  1  a  u  I .  I  •'  I  h I-  cm  ird itl.l I  e  -  v  t .  m , 

tin  ■  siirl.n  e  nl  t  In  I  hi  ill  "i  l  In  plain  al  si  at  ter  me  c.  de-i  i  ihei) 
by  1  he  point  t '  i  .s  in  a  e  h  i  u  ■  b' ) 1  -  tu'  b  l !  ,  - 

h  b,  .  ami  e  t  -  t 1 . .  i  aduis  ■  a  t  he  Imhliln.  <  '.aly  olat  mti  nl  t  he 
slope  in  nl  the  Mattered  i  n  ,;iyes  the  Inllnssauj'  par.iinelru 


iMpi.it  UMI  '-  lor  t  Ii*  vtr!  M.il 

". as c  front: 

,  -  M) 

M  *  ui")  1  |. 

lit) 

\  .  s 

"ill  (1  4  III  ')  1 

(■II 

6  br  *  /  h 

k-,/  - M 

U  '  ~  -  -  *  u 

HG  ?  flUj  •' ,*lu n »  ni  til**  lo*  ttw 

[*.Y)  of  n  potol  Onr  Mi.*  ■!  v ■  r l • :  t1  w  »v*’  Itunl  Tfio  uppiv  hflr  of 
/•  fj  pariiH*'!  U>  M'i-  »  . s v ,  ,  /i ,  .  «i t>,.»  il.-.Mtn  Outworn  tfv  p(T»>t  of  contai.t 
Wtlfi  Iho  t)uiif't*'  ■  i  I  t‘  *•  1 1. if  f  t  I  [h  •,'[  imi  Pv  point  L)  Tho  fOt.P.'  1 
iStr'iito  1*0'.  rn  s  t*<  itilruil'H.  i‘\t  to  f.M  ilit.ito  ,i  I.iylui  »»«tp.»f>si( *o  ol  tho 
W.iVO  Iff/'*1 


1  *:*>*> 


do 
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P2(Mi 


m  *=  tam/>  -  tnn|2  mis  'i/i/(i)|  (f>) 

The  Cull'lllal  ioll  ol  the  lllll  l.ll  I  lull  I  .  I.ll  dilated  llV  llci  rd ‘111;; 
virtual  wave  trouts  usi:n;  coordinates  to..")  which  have  lilt’ 
same  ori|;in  Iml  arr  rotated  ■  lockw.u'  liv  an  .mrle  i  .1  '21 
so  that  the  posit  i \  i-  r  axis  is  l In'  pal  li  wi  t  he  scull cred  1  nt  ical 
ray  predicted  l»v  rcometnc  upli.  ■  Kiptution  ill  and  ilia 
rotation  traiislnrinatiuii  rive  tin  tolluwmr  rxprcssiuii  tur  Ilia 
now  coordinates  nt  tlio  virtual  wav >■  I  runt 

u  -  x  sin'.’d,  t  y  m  '111, ,  Hi) 

r  -  —x  ons "ill,  t-  x’  sin'dl,  (7) 

i'u  obtain  a  I'aylor  oxpatisiun  of  Hu)  almiit  11  ■■■  (t,  tin-  Ini 
lowing  derivatives  were  computed  for  f)  -  0,  dritlst  -  /■  V; 
and  d  T/du  -  -  i/  d/i'i/d  t  i'ii  b/ij  1  'tin  wbi'io  i'll))  ih /tin 
and  i,  (d|  =  (hi /<l II  The  calculation  of; >  and  */  used  t  ho  r<!  1 
tiuli  (,  -  ri  sind  vvl.  10I1  is ovidoiil  troni  Ki<;  I  I  ho  ti'tin  nl  t In' 
Taylor  serifs  wliioll  is  linear  III  11  vanishes  bfinil.-.e  III  the 
coordinate  rotation,  this  is  evidc: a  since  pD',  1  11  and  1/P1,  1 

--  il  Oiisd,  ( tails',  inn  ul  oill  lie  .11  a  I  hi;:  her  ni  dor  tei  ins  in  the 
Taylor  series  (jives 

U  ■■  on/2,  (8) 

or  =  -  1/n  oush’,..  (VI) 

For  the  case  of  an  air  bubble  in  water,  numerical  calculations 
show  t  bat  Kt|.  (81  is  within  It)'1!  of  the  value  nl  1  to  I  prod  a  ted 
by  Ki|s.  tl )  l7t  provided  that  0  •  o  *' U07u  In  the  oali  iila 
tionolthedd  trail  ion  that  lulluvvs,  Ki|  t.8i  vvill  be  used  lor  the 
region  u  >  0,  (Airy's  treatment  ul  the  rainbow  also  retains 
only  the  leading  term  in  approximation  the  virtual  wave  trnnl 
but  in  that  case  it  is  a  cubic  ter  in.)  I'lie  sign  ut'  n,  mdieu'cs  tli.it 
the  reflected  wave  is  divergent 

II.  AMPLITUDE  OF  THE  VIRTUAL  WAVE  FRONT 

The  physical  optics  apprnxiiniitinn1  involves:  to  lie 
specification  of  the  polari/al mu  ul  ;l\o  nu 'it lent  wave  and.  i ill 
the  use  of  the  reflect  ion  enelfii  ion's  tor  a  pliuic  siirlace  to 
approximate  the  reflected  fields  at  each  point  on  a  curved 
surface  In  the  npplical inn  ol  the  teehmijoe  that  follows,  the 
incident  wave  has  all  elect r it  I ielil  pula  ri/; it nm  either  ent irclv 
perpendicular  to  the  jil.uie  of  -a  atlcrmc  nr  entirely  parallel 
to  it;  a  subscript  )  W  ill  he  as  .it; tied  !  he  value  t  ill  I  lie  1 11  I  ease 
and  2  in  the  second  I A  linear  com  hi  tint  mu  n!  id  looted  holds 
may  be  used  lor  other  incident  pol.iri/ul  mils  I  i'u  a  di  lint 
observer,  the  ret  lei  toil  wave  appear-  in  ,  mm  1 1  mu  \  11 1  n,d 
wave  describeil  liv  the  followiie;  1  oinplex  aeipbt ude  at  |ioints 
on  the  u  axis: 

I'll i,l)  •-  !■,(“'  “n.  Mid 

where:  1  =  \  —  1  ,  ic  is  the  'menhir  Ireipieiu  v,  /  is  the  time,  b 
=  2x/X,  X  is  the  W,1\ clcn;:l  h  in  I  lie  mil  1  1  an  dmni.  and  r .  is  I  he 
amplitude  reflection  cncllii  lent  ol  a  plane  siii'.iec  tor  tin  ancle 
of  incidetn  e  I1'': 

cos  II  --pi  -  in  Ih 1 

r,  -  .7  ,  Mi' 

Cos II  1  Ir,  -  no  ID  1 

-11 '  -  cos  il  -tin  -  -111 

ri  =  -  . .  ■■  ■  I  I ..  I 

n  ' cosll  f  (n  '  -  sin  ID 1 

Reflection  for  h  ’  b.  where  b.  -  I  ri  can  b>- describedfiv  the 
following  approximation  •  that  ::iv  e  tin-  h  ail  ui,:  dr  pi  uii  n.  1 


of  r  1  and  r  ,  on  1  ~  0,  ~  II  >  0: 

r  1  1  —  (8</n  cos II,  )I/J,  (111) 

r*i  ~  —  1  +  n-’(  Hr  In  cos 0,  )122,  (14) 

and  in  each  case  \<lr/(IO\  ••  1  '  l/-  Whim  II  >  II,  ,  the  reflection 
coelf  11  lent-  are  unimodular  ami  complex;  they  are  r,  -- 
expl  -n'p)  and  r  >  -  —  exp(—  iiV.)  wher"*: 

tan(X|/2)  =  (sin-f/  -  n  ’)i/2/ch»s 0,  (15) 

tan(X.i/2)  =  #i 2  tiin(A|/2).  Hrt) 

Rellectiun  for  b  >  h,  where  h  -  b,  «.  n  can  lie  described  by  the 
fol'nivini'  approximations  that  describe  I  lie  leading  depen 
deuces  of  7 ,  and  «/•  on  r'  -  H  —  II,  >  I): 

iS,  v.  (Hr'/ri  cos II,  t1  T  1 17) 

f>-j  ~  .’i  ‘{Hi  '/n  cosll,  ) 1  (IHi 

Approximations  for  r  in  terms  of  a  may  be  obtained  usin;; 
11  n  u/ipii,  )  For  u  '  (i,  Ki(s.  ( 17)  and  (1H)  give  ip  ~- 
,1,11'  '  where 

11,  ~  (H/llll  )|,2/(  osf),  ,  (19) 

tin  ~  nJfi  .  1  .ID) 

I'lie  specif  lent  ion  ot  '  '(11  M  is  Inrt  her  simplified  by  (a Ivilc;  /'  1 
c  -  11  lor  I*  II. ,  and  etpii  valent  I  v  for  11  <0;  the  use  ol  tin  - 
a  ppi  ix  1 11  ml  nm  is  mot  ivaleil  by  t  he  1 1 1  x  creel  ices  |  ilr/itH  \  and 
■  7, a- <  -  I  tl  and  i/i,  •  H),  respectively 


III.  FRAUNHOFER  DIFFRACTION  PATTERN 

The  scatterin'*.  ani|)liUule  at  law  distances  front  the  bubble 
is  proportional  to  the  Fourier  transform  ofr(  exp(--t/,'mi  1: 

/,(r))-  f  (21) 

where  1/  -  ./.,  -  </i  < ait  ancle  of  the  observation  point  relative 
in  )  is  a  '-aimed  to  be  small  so  thm  -ini)  ~  t/;  and  constant 
rna>;inttidc  multiplicative  factors  tii.it  precede  the  integral 
have  been  omitted.  The  approximation  described  in  the 
precedin''  ,.ei  turns  yield  the  following  exp  essiott  fur  /, 

/,.  -  i  (  e  d*i"i<!«i"i*i),,<1- r,\du,  (22) 

»/n 

wlii'ii'.  here  and  m  the  cipiations  that  follow,  the  upper  si,;n 
is  lor  polar  1 /a  1  mil  1  1  and  the  lower  one  '  ,  lor  )  —  2. 

Kvaliiatiuii  ul  I-;.,  C'.D  i-  faeiljtiitei*  liv  a  (  bailee  of' variables 
•  bat  i  c iinplcti  s  the  |i.irti.il  -ipiarc  loriued  by  litu':  )  iju  I 
I  )|  line  .  1  lo  \ ) 1  0  -  It  where 

ic  ■  >)|((i/X  I  (  os H,  | 1  1 28) 

We  will  approxiiiiale  the  integral  onlv  for  0  <  <  <8,  vvltii  It 

corresponds  to  the  region  predicted  In  ;;i  ontelric  optics  to  be 
stronr.lv  illuminated;  the.  rest rict ion  (lives  11  ■  >0.  F.ipiatiun 
(22)  may 'hen  lie  written  as  lolloWs:  /'y  -  i  t Xu  cosff,  ) 1 '  (,7 2 , 
where 


^  -  ,.d  mil  t‘* i..| 

(24) 

1,'  t.Mi  1  -  d/.’  '  jl"  ‘  ; 

1 1  -  11  1 ' '  J  ( X  u  cos 0,  )  1  1 

(2  hi 
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In  the  exponential  that  preeetles  t lit*  integral,  t/q(ie,(l)  is 
evaluated  hv  the  appropriate  variable  substitution  in  Kq.  (  JM, 
(0  replaees  ic  and  w  replat  os  .•  I;  phvsieallv.  v’ ,  be.O)  is  the 
appruxitnale  phase  sliitt  associated  with  r,  til  the  penmetric 
ray  seat  te  red  to  the  obsei  vat  ion  a  Hide  The  second  term  in 
this  exponent  nil  arises  I  roll  i  t  lie  t  mn  plot  ion  ol  t  be  si  pin  re. 

We  now  examine  the  condition  tor  omitting  the 
exp(tV  , ( )|  lat  lor  xmI  Inn  t  he  inle;:ral;  when  t  Ills  sitnpld  i 
cation  is  perniissdile,  the  inlcpral  mav  lie  estmialed  iisinp 
Fresnel  integrals.  II  1}  ,  I  .-.it  I  is  a  Mill  n  ientlv  slow  I v  xiirvio;; 
function  of the  prim  iple  eoiilriboi  ion  of  the  intcprnnd  in¬ 
curs  when  |.’|  7  heeause  the  phase  o|  exp|i  (  n  issta 

t  ionary  at  j  -•  0,  ( For  a  lucid  discussion  of  slat mnarv  phase 
approxitn.it tons,  see  Scpcl  ')  For  it  ■  0,  ■  y  it ,  and  z  ?_|y|. 
it  is  evident  t  hat  | y' ,  I  Mr  1 1  --  i/  ,!;■  .1)1,  copseqneiit  I  v  the  fat  lor 
in  ipiest  ton  mav  he  omitted  it  the  rad  tie.  nl  I  lie  In  i  hide  is  si  it 
ficiently  larpc  t1  at  i^  til  «  I.  For  example  water  has  n 
~  I  ltd! l  anti  visible  lipht  from  a  lit  \e  last  r  lias  a  wavelength 
in  water  of  n  1  IM'd'iHttm  F’nr  polarization  1  the.  etmtlo  ton 
till  i/,(-.th  |P vox  tm 'fi!l  pint1  1  >s  I  and  ha  polarizal  inn  P  it  ptves 
(n /fit M i  pm)1  1  ■>  I  where  F.t|s.  ( I'.n,  r.'tt),  amt  (Phi  have  also 
licelt  Used 

Alter  some  r ■  - . 1 1  r at u: i tut .  the  result  ol  this  approximation 
max  he  written  a-  Pillows 


f,  =,  (.i|  till"  ")ilo,  )  ..  .  .  >]  (PI)) 

where  Flu)  is  FrcsiuTs  i m t  r;tl ’ 


and  fl--")  =  S'l  ’>)  =■  ■ '/>.  For  any  stale  of  incident  pn 
lari/al  ion.  llle  billownu;  fat  Inr  pixes  the  llepentlenee  of  tlm 
intensity  on  > /  where  i;  ‘  is  the  complex  onpipate  ol  y: 

HH*  =  ((’(n  I  -I  | •'  I  |.S'(i.  I  t  'h\ (HO) 

The  normali/al  ion  used  in  Kqs.  t I )  anil  IT')  is  such  that  yy* 
•  2  as  P/>,  -  t,’>)  heetillies  lar'.e  anil  positive  As  t  .xpe.  led.  yy* 
'  0  as  </>  t/i,  het  tunes  lar::e  ami  positive,  Imwexer,  t  lie  repinn 
It1  <  1 1  is  out  o|  I  lu-  r  a  1 1 : ,  •  ■  III  lilt  x  al  h  III  V  I  esl  lor  1  lie  approx  i 
mat  ion  used  in  dertx  mi:  Kq  l.h'1 

Davis’  anahsis1  ol  I  lie  peoinet  l  ie  ..  at  tot  me  mav  be  used  to 
obtain  tbe  ma  nuli/almn  leqtiipd  to  expre  .  the  absolute 
intensity  /  associated  with  the  n-llct  let)  lav  at  a  distance  /« 
»  a  troin  tin-  t  enter  ol  the  Inilihle: 


/  ~  /„(ti/AT  /"’VS,  (Ill ) 

where  f,i  is  lli"  lilt  itlelil  mtciisitv  and  *,';,’*/8  ►  V;  for  lurpe  n 
so  that  Kq.  l.t  I  I  i  nun  it  Its  XV  It  It  Kq.  t  a  )  ol  I  tel  1  111  t  he  appro 
priate  limit  Furl  her  more,  it  tin-  Intel  It  i  tin  e  he  I  ween  the 
totally  reflet'! eil  lavs  and  the  multiply  relict  ted  ravs  are 
omit'ed,  the  total  ilitensitv  /,  Ini  nines 

/,  /ntd'/t)  '|;:p'/8  t  (,  l,,.)|.  Cl'.') 

where  (t  (<.'i)  is  the  total  emit  ri  In  it  ion  hv  the  mult  iplv  relit  •tied 
rays.  For  an  ar  huhhlc  m  water,  lilmnin  it  it  I  hv  nnpolari.ed 
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HC!  I  View  td  II ii -  appat.iliis  in  th»-  pi. me  of  scat’.. ’tint)  which  was  ht>ti- 
/ontal  m  llm  lat«  -i.iti  rv  The  hollnm  ot  the  tmtitite  was  m  .  oiilat.t  with  a 
n.ietlle  ttial  pievt'niea  Itie  h-mhli  (torn  ll'i..!iiit|  out  ot  the  sc.ittenn*.}  plane 
The  top  ol  the  hublM.-  (wtm  u  e.  the  suit,  displayed  here)  was  viewed  Oy  a 
microst ;opn  The  wavelentiai  m  ait  ol  tlm  iliominatioii  was  bJ2  8  nrn 


lip  lit .  'Calile  I II  ol  I  tel  I  shows  that  (>'  is  n  sun  ml  Ii  fund  ion 
til  t/i  lor  o  neai  '  -ind  that  (I.OOPT  and  f»" 

(7.rt°)  -  mm  ('  ( 'onsetpieiit ly  (i  lor  an  air  liuhhle  in  water 

is  approx  i  mat  t-b.  1'.  to  "I  if  I  lie  asymptotic  value  til /,'i;  "  8 

Kt  plain  m  .nil  anal. n  m  |.  a  m  to  I  he  wellknoxvn1"  Fresnel 
(rienr  /(••/'/ 1  <litlt.it  t  .  ui  hv  a  si  i  .up.  It  I  -eilpe  where  the  pa  rain 
elei  a  is,  m  1  hat  i  a  .e.  proportional  t < >  the  lateral  displacement 
ol  the  observation  1 1. . it > •  ‘.iiir*  1 1 1 1  cdpc  of  the  peoinet rit 
shadow  For  both  ease  m  I  lie  repioti  a  >iltlv.  intensity  is 
predicted  to  o- -a  ill. ate  wit  l\  an  used  la  lion  amplitude  t  hat  tie 
creases  rapidly  with  mt  re.e  inp  if.  This  behavior  is  predicted 
for  the  fut  hflii  alt  mu' v  of  a  huhhlc  heeause  of  the  abrupt 
variation  m  the  aa  phtudf  of  the  virtual  wave  front  near  the 
critic  ill  anjde  and  bet  an  .e  ol  I  lie  <  ir/  rtiTirre  of  this  wave  I  runt 
induced  by  tin  rell.s  tnm  Iro'ti  the  curved  surface. 

From  the  aloretuenl to.  ctl  hound  on  |.*|  lor  the  principle 
emit  nl  nit  ion-,  ol  the  mteci.il  in  Kq  r.M  I,  we  call  est  limit  e  lot 
th"  ic  rep  i  on  of  interest .  the  errors  introduced  in  t  lie  rah  min 
ti.ai  nl  the  I,  trom  Kq  ■  ( 10)  anti  (.')>).  In  the  interval  U  ■  n 
<  :t.  Kq  (Hi*)  pi  '‘In  i  ■  two  intensity  inaxitmt.  In  this  repmu 
the  pi  I’leiple  t  tint  i  ihii' ions  tietiir  with  0  *' u  0  -  i'  ) 

nil-  A  ‘  to)1  ’with  ir  •  :t.  With  the  ntoiemeiitium  d  A 
mill  willt  n  I  n  il.  the.  estimate  pives  ~  Mt  pm  A 
com  pa  i  istui  ol  Kq  i  i  7  i  with  A ,  ealeiil.it  ed  usinp  Kq.  1 1 ' '  i  and 
this  n  show-.  I  bat  I  In  ei  i  or  in  1 1  toll  n  ed  ill  i)\  is  less  than  TV  A 
similar  com  pa  n  iml  F  q  .  '  18)  and  (p(l|  pivc  a  In  mud  lor  t  lie 
error  in*  "I  s'  <  ’oic-cipn  ntlv  the  use  o|  a,  -  i  ip/1  is  jus 
t  diet  I  Kta  I  In  i  ti  s  .1  e  tine  a .  a,  I  lie  upper  limit  ol  ill 

tepriition  m  |  .|  i  n  lv .old  nvc  a  ir  etui  approximation.  The 
main  jn  1 1!  n  a'  i  1 1  I . .,  'In  qipt  ox  imat  n  .ns  dcst  rilietl  in  Sec-.. 
Ii  III  is  tin  .  niii!M,i  .mi  with  experiments  described  ill  llle 
set  t  ion  1 1 1 a r  n  il.  ...v 

IV.  OBSFRVtD  SCATTF.RING 

I  he.  m  i  lion  ■  1 1  ■  *  1 1  it  phot  op  i  Mpluc  nliscrvnt  lulls  nl  tin 

(T  it  teal  a  lip  1 1  X  aii.it  mils  111  I  lie  st  at  1 1  ret  I  ilit  e  Its  1 1  V  |>l  ft  I  it  ted 
in  I  lie  prt  vi on-  ■  i  i  1 1 1 1 1 1 .  I  iptn v  :t  ..hoxx  -  ,i  t|  iapraut  ol  t  lie  ap 
p  trains  I  Ii  l  illetl  wat>  r  was  placed  in  an  aluminum  eon 
tamer  The  mm  .  wall  .  ill  the  emit, liner  hail  been  lilaekt  neti 
So  that  I  It.  x  ab  .  n  In  1 1  miu  Ii  ol  the  si  rav  lip  lit .  A  steel  licet  lie 
passed  ver a  i.  ill',  t  h  roup  Ii  a  rubber  septmil  ill  t  lie  but  ton  i  ol 
t  lie  font  uln  i  I  In  o!  In  I  end  ot  I  lie  lit  t  tile  was  col  met  letl  to 
a  s\  ’  mpe  that  i  >  a  it  an  nil  air  Horizontal  boles  IP  tin  in  1 1 1 
a  met  t  r  pa-  .  .  i  I  n  i  mu.li  tin-  w.dl.nl  tin- 1  mit  aider  a  I  riphl  an 
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FIG  4  Photograph  iilltwi  colic, il  nnplo  l.ai  ilillor.  tux'  I'ea'  a  ie-  <• 
hands  hnva  a  chararteosho  separation  ot  appioxici.ilnr,  anna  'a*  v  n.  ii 
polnn/od  with  its  electrii;  fink)  parallel  to  the  si  all . hi  pl..ai 


mu,'  .rim..  ,  a.a  ic.isa,  (loci'  loft  to  right  Tlw  vertical 
imi  ,,ty  ,r,i'  i-d  hy  dillim  lion  The  incident  beani  was 


pies.  These  holes  were  sealed  wit  It  O  imp.  in  | >r*  ■■ » •  m t  w.sta 
leaks  and  were  coy  ered  by  ant  i reflect  oat  <■  >ol is!  v  uni  >  v  A 
bubble  eould  lie  rrealed  al  the  t'lal  lip  "I  'he  m  edle  b\ 
crcnsinp  the  volume  of  I  lie  syringe,  this  iiul'hle  wa  •  slahh 
provided  its  volume  was  kept  sufficiently  mall  Tin  b,ibhb 
was  approximately  spherical  in  it.  •  I'lH  r  a1  >  ■  •!’  1 1*'< >f 

diseussion  of  t  he  stability  id  bubbles  -.it  spent  lei  I  I . .  ■  d1-  ■ 

see  Kef.  II.)  The  base  nl  t  lie  bubble  wa  ■  m  h  ml  a  a  w  it  h  l  I, 
needle  durinp  the  obsery  at  ion 

The  polarized  beam  I  mm  a  h  m\V  lie  Ni  1  i  i  w.e  .  > 
pamled  to  produce  an  S  nun  diam  parallel  beam  Hie  es 
ponded  beam  passed  normally  and  horizontally  ■  hn.uph  u\u 
of  the  windows  in  such  a  wav  that  the  central  lavs  .1  i  lie  In  am 
passed  close  to  the  center  of  the  bubble.  A  enmoi >.  lewed  lie 
bubble  thrnuph  the  second  window  The  needle  .uni  the  inner 
wall  of  the  aluminum  container  were  blackened  in  ri  hu  e  tin 
intensity  ol  stray  liplit  that  would  barentl  rw  i  .e  ob  .eme.l 
t  he  scat  lerinp  ol  t  he  bubble.  I  .i  it  lit  seal  lei  ed  bv  I  lie  buhl.  It 

near  I  he  erit  ical  scat  teriuk  anp'le  had  an  alible  ol  im  idem  <  at 

the  water-inwdow  interlace  ol  approMinately  U"'' 

This  liplit  is  retracted  at  the  window  air  mlrilace  I  hi 
camera  telephoto  lens  (Nikkor  if.  manulactnred  hy  Nikon 
I nc. )  had  an  ef  fect  ive  local  lent; I  h  /  ol  '.'i IU  mm  I  lie  lens  yy  a 
placed  with  its  axis  horizontal  and  approximately  p  i'  alli  I  t. 
the  refracted  i  ritii  al  tavs  Tin  pol.u  i/.Umii  ■  In  lieidee' 
monochromatic  I  ip  lit  was  usually  i  boson  wit  h  i  lie  elect  i  a  tut' I 
parallel  to  the  seat  ter  nip  plane  l  bat  mi  respoM' I  i  ■  ■  'mill 
preceding  analysis  This  eboii  e  lor  t  lie  pi  .Ian, -a i  am  pi-e.s  t  In 
preatest  mapiiilude  lor  the  cool  I  a  lent  ol  t  1  ui  «  • )  'll' 

( 'olisiapienl  ly  it  pives  the  most  I.-Ipld  decrease  tm  |/,|  with 
inereasiiii:  •  and  it  should  la  the  polari/.at  ion  ti  u  w  lu  eh  lie 
omission  of  t  he  wave  fri  mt  wit  Ii  u  ’'Din  Kip  i I  will  pindie  e 
the  smallest  error  in  the  piedi<  i”d  seat lerinp  I  lie  diameiei 
nl  the  bubble  in  tin  ■■rnttfrim:  /i/am  yy.i  n  "  med  with  a 
mieroseope  located  m  the  air  directly  ahove  the  buhlili 

I’llntopraphs  were  taken  with  I  lie  camera  leu.  loeuwd  l" 


I, in'  'h  . I  'h.  1 1 i 1 1 1  plane  was  at  the  optical  transform 
i  i  i  a  i  io  h  ie.  I  .el  if  ■  .a  ...  where  and  </>, '  denote 
no  nr  .  ,  i  lie  i. om  o  n  .ide  the  cell  i li.it  correspond  lo./i  and 

ui  •  la  w ,c. a  . . i.|"i'iii|y  iiieremenls  in  the  distance  / 

n  I.  I  1 1  I  i  y.i.vh  ,  sep.iliye  are  related  to  the  increments 
■  I.  \  A  /' im  \l  ■  ■  I  Numerical  calculnt  ion  ol  t  be 
,u  a.  i  :;l a  ami  M  e  water  air 'iiterfaces  show 
a  a..  .  in  in  a  m  c.  Im  related  to  the  eoirespiiiidinp  Aij' 

ie.  'ie  o.'iiiv.nh  a  nepiipible  error  The  success  of  this 

i,„  i  h... I  n  I . .  the  choice  ot  the  plane  of  the  viewed  window 

liih.il  ii  o,  m  at  Iv  perpendicular  to  ( lie  critical  rny.  These 
it  .uli  .  i.a.ay  hr  combined  topive  Ai;  *  Al/nl'  [C)f  course,  if 
i  'a  i  ot  it.  ti  ei"e  ext ci  mil  to  the  sphere  had  been  uniform  so 
■  ha  I  io  ■  x\  "oh  iw  v.  a  needed,  it  woiih I  not  lie  necessary  to  ill 
\  im  ,  i  a  i"  ,v  m  mi  a  ai .  Then  A/  'nh'  could  lie  used  as  the 
1 1  to  i  ."a  nsi  iiiieier  in  Ki | .  r.’li  hectiuse  of  the  Fourier 

'tail  I .  n  oil  1 1  l  if  Opel  t  Ie  .  1 1|  a  lells. 1  | 

i  i  e ,  I,  .,  |  a  "it  iy  e  i  e| 'todm  l  loti  ol  a  [ihotopraidi  ol  t  lie 
III  I,  al  .inch  -'  U  lei  ii  ip  yvil  Ii  a  bubble  diameter  of  TI  5  -fc  ().(), ri 
a  a  a  I  he  •  ..I,  i  odii  ales  a  Aij  i  el, at  ive  to  an  nrbit  rarv  refer- 

,,,  e  n  I  im  him  used  was  Tri  X  with  an  exposure  time  of 
I  i ■.  *  t  he  hi o  el  \  ei  t ii  al  dark  hands  xvit h  a  characteristic 
pa.  to:  l  ipi  i  Miualelx  'r,  mrad  (  —  I  1°  I  are  examples  of 

ije  mi . .  ii. mini. i  predicted  hv  K<|.  (Ml!)  The eoncent ric 

,  o .  ol  i  ■  1 1  on  a-  tbit  at  e  c  entered  close  to  the  most  cent  ml  dark 
tire'  ,  i  •  appat  i  nl  l\  ,m  ,u  1  it  act  ,ass.  a  anted  xvit  ll  spurious  re- 
t  ie, a  a  -  n  mt  hei  I  tom  i  he  exit  rvuidnxv  ot  from  the  telephoto 
a  o  I  ai  il  oh  tar" at  i in  the  lelt  ol  the  lip. me  is  indicative  ol 
lit"  i e. 1 1 a  i ion  in  intensity  that  oeeii rs  when  le  <  0. 

Mi  t . „ t*  tisit omet .a  cans  were  obtained  from  two  negatives 
i  ,  1. 1,  do  . ii.  a  ipi.ioi  Mat  iy  e  comparison  yvil  Ii  the  theory. 
I  i,  ■ .  oinpaia  s  i  in  relat  iye  intensity  inferred  from  a  scan 
:  i  m  pal  i  "  lor  a  buhlili  diameter  ol  I  .lib  1  (I. Of)  Him.  Two 
liei  p.u  iiiii  h  i'.  yye.e  used  lor  each  nepative  to  make  the 
"  1 1  o  ■  .  ai  w  1 1  h  t  lie  I  hem  y  possible:  til  The  local  ion  of  the 
,  .  i  .no  ,  1 1  al  1  ",  cmild  not  lie  determined  experimentally 


.1  (lot.  See.  Ami  ,  Vol  lit).  No  ft,  Sept.  ciIht  i'l.'i 
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Ifomn  microdensi'oiiwter  sc, in  ol  Hie  photi  implm  neu.ilivo  IP,.  iiiiUMo 
diameter  A ,  IS  I  00  mm  Dll'  solid  I  wve  I  .  nn«n  I,v  t  q  I  no  which  r.  tin- 
result  ot  ttie  ptiysii  ,il  '■■ptii.s  .tpprnom.,11011  I  ai,  i  ai.iuieltirs  wen1  .kJiusIiki 
to  facilitate  this  comp  insnn 


so  that  an  arbitral  v  constant  angle  w.i-  added  to  I  be  abscissa 
of  the  measured  intensity  profile  in  optimize  t be  agreement 
with  the  theory.  I  h )  Tile  ah-.oliile  intensit  x  could  oot  lie  tie 
tormined  experimentally  so  that  t h  experttia  -iit.it  mieositv 
/,.  was  obtained  I  is  mull  i  pi  v  in;’,  mrasiin  men  is  by  a  *•<  nistant 
factor  K  that  seas  selected  to  optimize  aitteement  with  Ivp 
(HO).  'I'be  mat;  n  it  title  of  (» '  was  sill  ficiettl  Iv  small  t  bat  it  was 
neglected  in  the  cnmpulalion  of  the  theoretical  prediction  hi 
that  gg*  could  be  used  as  the  prediction.  The  relationship 
between  the  measured  transmittance  of  the  negative  >  and  I, 
is'-'': 

/„  =  K(r  -  r,,1'1)  (Tt) 

where  fo  was  measured  at  an  unexposed  region  ol  the  negative 
and  the  contrast  index  >  was  estimated  I  mm  manufacturer's 
literature  to  he  0.70  for  the  conditions  of  film  development. 
The  r ( r ;)  measurements  were  smoothed  Ity  hand  to  reduce  the 
magnitude  of  the  circular  fringe  artifact  in  the  itilcmd  /, 
The  scale  of  the  smoothing  only  s  1 1 • ; 1 1 1 1 v  exceeded  I  lie  spa-  up; 
of  those  fringes.  The  validity  ol  I  he  pmci  dme  was  su list  an 
lialed  by  eomparinj;  parallel  but  dislnn  !  scans  front  the  same 
negative. 

Figure  ft  illustrates  ..  close  .igrcomciil  ol  t  be  relative  pou 
t  ions  of  the  maxima  and  minima  ol  I  lie  mini* it  \  band-,  w  il  b 
those  predicted  by  !•’,< | .  (alii,  l  ie  a  values  wi-ie  calculated 
using  Ki.|s.  0..0  and  with  n  I  d.tb.i  and  'u  I  i->  mm; 
l'(lc)  and  .S’prl  were  obtained  limn  table-  The  values 
and  locations  lor  the  maxim,  i  and  inimni.i  weir  obtained  I  root 
Fresnel's  original  l  a  bill,  it  ion  tm  i  be  de  i  i  ipt  ion  ol  edge  dll 
fraction  (as  iptoted  by  N •  1 1 . i » 1 1  1 )  There  is  an  apparent 
discrepancy  ol  I  lit'  I  hem  >  w  il  li  t  b.  I.  .11.  ,w  in-  I  eat  ores  m  i  be 
measurement:  III  there  appears  to  be  a  diilerem  e  ot  mini 
sities  ill  the  region  n  '  I  Mini  eauiiol  I--  explained  by  l  be 
omission  ol  ( ’  l  rum  the  ealeiil.il  nm;  I .  i  •  lie  a  in  pi  it  ode  ol  the 
intensity  oseill.il  iolis  cm  eed  -  I  lie  pred  i.  I  lull.  I  ol  n  •  .mil, 
( iiil  I  here  is  an  unexplained  drop  o:  !  be  i utensil  y  near  ic  :t 
Allot  her  neg.it  i\  e  .limvs  ipialil  at  lx  ely  ii.iil.il  ili-ngl  eeliiciit 
lull  also  a  similar  .e/rei-meul  lor  i  lie  In.  a  I  ion  ,  1 1  I  be  maxim.i 
and  minima  with  the  theory.  I’be  >pi.uiiiiv  end  ipi.ililv  ot 
measurements  that  ran  be  obi  allied  w  il  li  1 1 ,,  a  all. able  appa 
rat  us  are  not  adoi  pi.ile  lor  del  er -uiiiiii  m  l  In  ■ .  b  I  eaei  eu  . .  a  ibis 
ealeiilal  ion. 


A  comparison  ol  photographs  taken  with  bubbles  of  dif¬ 
fering  radii  shows  that  the  angular  spacing  of  the  intensity 
oscillation  decreases  with  increasing  a,  which  is  ex|*ccted  front 
the  torm  ol  K<p  C.M).  This  decrease  is  also  evident  to  the  eye 
if  the  radius  e.  increased  as  the  bubble  is  viewed  throu„h  a 
teles* ope  locusi  d  mi  infinity  Intensity  variations  similar  to 
thoge  displayed  in  Figs.  I  and  b  were  also  observed  with  the 
incident  polarization  u (  j  -  1 . 

V.  FINE  STRUCTURE  IN  THE  SCATTERING 

Card  ill  examin.il  am  ol  the  photographic  negatives  revealed, 
in  some  ium-s,  an  om  illation  in  the  intensity  periodic  in  ij  wine  It 
was  superimposed  on  the  diffraction  stmclurc  described  in 
See.  IV  1  '  For  example,  osi  illations  in  intensity  with  a  period 
ol  appio-  -il.  i  I  'i  :  mi. id  are  visible  in  the  original  photo 
graphs  j  •  i  mi  ■  ■  ha  e,l  hi  l-’ig.  I.  (It  is  unlikely  that  this 
stun  lull  X  !!  I"'  •'  "  in  the  printed  version  appearing  ill  this 

lourii.il  lie i  n  ..  *  .'i.  o lx c  ,  I  in  Fig.  b.)  Kxaniinat  ion  ot  t be 

neg.ii.e  I,  -  ill -a  1  la  1 1  Mige  spacing  r;  decreases  wit  li  in 
ere.e  in,  •  i  .  •  m;  tin.  loseupe  w  as  used  to  measure  >| 

fr.. in  lb.  leg  -a.  e.i  I  mlibb".  with  1 .7  mm  2«  $  -.2  mtn. 

I  be  i,  'il:  u.  i.u  i  *  mud  be  fitted  by  the  equation,  i;  = 
\i.\  .ii.  \  i  .*  I  i  o  i,  however,  it  was  difticiilt  to  oh 
lam  *  mi  a  :•  . '  i...  ••  lueineiils  from  each  negative  because  o! 
am  a  .'ii.l  l-  '  .  1. 1  i.i  I  ol  the  oscillations  Knergy  asso 

da i eil  will;  •  •  i  mi  i  i\  .  incident  with  h  <  h,  will  eventually 

lei.xi  if.  '  d'1  '  ,'ae,  i .  1 1  with  iji  From  the  linearity  ol 

I  lie  i  ip  i.i  l  e  i"  oid  i  lie  magnitude  of  A,  i  I  is  likely  that  the 

oseillaiae.  m  *.m  ,  d  h\  t  he  interference  of  (hut  light  wit  h 
the  ei  1 1  a-. all ,  i .  :  le  * .  •  i  light ,  however,  the  nu  in  her  of  rel'lee 
t icin'  and  x  no,  '  a  1 1  ,e  path  of  this  interfering  light  has  imt 
lieen  .!•  •■I'o".  I  l  i."  i  p  nation  ol  tli<>  maim  virtual  source: 
foi  I  .  I  \ 

Si  mi  I.o  .  i  I  ,  ’ ,  i .  a  ed  m.i  illations  can  he  observed  in  the 
rainbow  '  t ■  " in'  imm  liquid  drops  with  il,  >  zt, ,  and 
monoebt "I..  e  i 1  am n  a  loii.  The  author  has  found'"  that 
l licit  ■  pa i  it.:;  1  , -  .  et-i'lv  propoit ional  ton  and  for  drops 
of  xyleiii  '  "u  to  water:  A  ~  t  Mi.  These  fringes  were 

found  loh'  o'!  i  deteei  mg  uiicroii  amplitHile  shape  o.-i' 
dilation  -  a  o  I'm  ■  i<  i  -!i mieti  drops  since  changes  in  t  lie 
drop-.  ■  i  1 1  ■*  i  ni'i  ui  mgiil.ir  shift  in  tin  position  of  these 
Irin, a  ■"  '  on  . -.jin  ol  1 ,  wbei!  the  shape  of  the  drop  os*  il 

late, I.,-  i  i  j  a.  i.  •  a  ,i  o|  a  lew  l  Iz.  the  fringes  appeared 
to  bln i  i"  1  i  ,  a  il  I.  •  illations  ot  the  shape  ol  I  lie  bubble 
should  .  •  the  iiie  si rud ur(‘  lo  blur,  luileid, 

bill)!,!'  1  in-  .  ii,.'  i .  n  e  ito  I  by  room  viliralions  may  b.ixe 

bio  1 1 1  - 1  i  i"  '  ■  a  tin-  e  negalives  wliel'i  none  were  i  i, 

•*et vi  -I  t*  »p.  -  .'"a  a*  re  typically  I /HO  si.  The stniilaiitn  « 

ami  dilli .  '-:n  d  angle  scattering  xxilli  rainbow 

si  at*  i  mi"  •  -  ■  1  n  I-.  '  mi  pa  ring  Fig.  -I  with  photographs 

ol  t  la  1 : 1 1 1 " i  1 "  i  ■  i  !  i  .  i  u 1 1  I 


VI.  DISCUSSION 

The  ei  ilu.il  n  do  seal  t.  ring  of  eight  bubbles  was  photo 

graphed,  the  'mlido  . . 'lei  .  were  m  the  range  l.ti  '.Ml  mm 

K.ii  1 1  pin.'  1  i pi-,  i  a i ml.ir  in  appearance  to  Fig.  I  in  that  die 
ring  m  ;  "i  ml  ej.  ,u  x  <  me  lo  Ibe  or  it  u  at  angle  is  1  |Uiie  visible  I.u 
eithei  tale  I  I  nil  id'  ol  I  ."la  I  IZ.U  ion.  Consequent  l\  tllcle  !  - 
i  ,  n  1 1 1,  , 1  - 1  ■ ' '  l  e  a  i .  i  I  be  pi  iul  id  ion  ol  geoniel  l  ie  opl  n  - 
m  ,ir  1 1  i  .  i  ;  ’ :  ' .  nd*  I  or  ai  i  I  a  ibl  lies  ol  t  In-  ,i/e  il  i  xx  a  1 1  a 
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Furthermore,  the  measurements  of  the  positions  of  th(>  in¬ 
tensity  maxima  and  imnima,  measured  relative  lo  an  nude 
lermined  angle,  are  in  reasonable  agreement  with  K(|.  (dll) 
An  experimental  lest  ot  several  aspects  i il  |  his  ealeulat ion  mav 
re<|uire  (he  use  ot  photoeleet  lie  del  et  t .  ns  heeause  ot  limit  a 
tiolis  in  the  present  appoint  u  r.  These  i  lie  hide  t  lie  following 
features:  (i)  the  amplitude  nf  the  ringing  ol  the  intensity,  hi) 
the  absolute  positions  in  the  maxima  and  minima,  (ini  the  rale 
of  decay  ot  the  intensity  for  te  5  1,  (iv)  the  rang  •  ot  radii  for 
which  K(|.  CUD  may  he  used,  and  (v)  qnanhl.i"  i.mire 
ments  of  the  interference  phenomena  describe  .  i  see.  \ 
Unfortunately,  the  largeness  of  the  huhlile  radii  makes  it 
impractical  to  calculate  the  Mie  scattering. 

A  brief  discussion  of  this  use  of  the  physical  optics  ap¬ 
proximation  is  in  order.  Surface  wave  phenomena1  do  not 
appear  to  he  included  in  this  calculation;  Fir;.  f>  is  evidence  that 
they  are  not  needed  tor  an  approximate  description  ot  thesi 
observations.  Surface  waves  can  he  included  in  the  Watson 
transformat  ion  method  which  has  been  used  to  I  muni  late  the 
scattering  from  a  cavity. Is  l"  A  comparison  with  the  results 
in  Fitf.  ft  has  nut  been  performed  heeause  of  I  he  mm  plica  led 
nature  of  predicted  sealtenni;  amplitude  j<w>  Kef  1H,  Kps. 
(4.92)  and  (  1. 101  )|;  mde"d.  the  present  author  h  as  been  unable 
to  determine  il  the  ringing  is  ei insistent  with  those  amplitudes. 
A  secnrtd  shortcoming  of  the  approximations  presented  m 
Secs.  I  III  is  that  the  curvature  of  the  huhlile  orthogonal  to 
the  plane  of  sraltei  ing  is  omitted  until  Kq.  Ctl)  where  its 
geometric  effect  is  included  in  the  «•'  factor. 

Diffruction  may  have  a  major  effect  on  the  far  field  up 
pearancc  of  the  erit ic.il  angle  scattering  for  hiilihles  in  water 
illuminated  hv  collimated  white  light.  I.et  </>, '  mid  </■, '  denote 
respectively  1 1>,  tor  violet  and  for  red  light.  The  dispersion  ot 
water  and  Kps.  (I)and  (2)  give:  </>,'  —  0,  "  ~  17  mrad;  con¬ 
sequently  the  geometrical  scattering  (Fig.  1 8  of  Kef.  I),  as  seen 
through  a  telescope,  should  appear  bluish  in  the  region  </»,- " 
<  0  <  0,.'.  It  is  evident,  however,  from  Fig.  5  thill  diffraction 
smearing  of  the  otherwise  abrupt  change  in  intensity  tor 
monochromatic  illumination  will  tend  l>.  reduce  significantly 
th'  coloration  predicted  by  geometric  opt  ies 

This  work  was  motivated  hv  an  investigation  of  optical 
methods  for  measuring  the  refractive  index  of  superheated 
liquid  drops  suspended  in  immiscible  host  liquids.  At  suffi 
cicntly  high  temperatures,  such  drops  frequently  have  n,  < 
n„  and  they  should  exhibit  critical  angle  scattering  phenomena 
similar  to  that  described  for  hubbies.  Kvon  though  the  dis¬ 
tinctive  critical  angle  scattering  predicted  by  geometric  optics 
is  altered  by  diffraction,  ft  may  he  interred  from  the  absolute 


angles  of  the  intensity  maxima  and  minima,  provided  the 
present  (henry  is  found  to  give  I  he  ubsolutr  angles  when  ti  is 
known. 
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Scattering  by  bubble  in  water  near  the  critical  angle:  inter¬ 
ference  effects  (P.  L.  Marston  and  D.  L.  Kingsbury)  0.  Opt. 
Soc.  Am.  71,  192-196  (1981)  [An  Erratum,  reproduced  below, 
was  published  in  0.  Opt.  Soc..  Am.  71,  917  (1981).  This 
notes  errors  of  transcription  whicFTare  not  corrected  in 
the  reprint  given  here.] 


errata 


In  the  caption  to  Fig.  3  of  our  paper,1  i)  m  <t>  -  </>,  should  be 
written  *)  =  </><  -  <l>,  which  is  the  definition  of  >j  given  in  the  text 
of  our  paper.  In  the  line  below  Kcj.  (lf>),  *■  (S/nnl'^cosfl, 

should  be  written  ns /l  i  =  (8/nn)l/l!/cos«,.  Also,  the  top  line 
of  the  right-hand  column  on  page  194  should  be  written  as 
", . .  error  from  using  Bq.  (22) . ,  . 
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A  physical -optics  ii 1 1 1 in >x i inn t  h >n  is  derived  fur  light  scattering  by  dielectric  spheres  wit  li  refract  in-  indices  less  l linn 
their  surroundings,  and  it  is  applied  In  air  liulihlcs  in  water.  The  approximation  gives  the  coarse  structure  ill  the 
scattering  when  the  scattering  angle  0  is  near  the  critical  scattering  an)  le  t/>, .  where  0,.  -»  H.'t°  lor  huhliles  in  water. 
Diffraction  has  been  observed  to  he  important  in  the  critical  region  because  of  an  abrupt  change  in  the  amplitude 
of  the  relied ed  wave  |1*.  I,.  Marston.  .1.  Opt.  Sin-.  Am  1*9,  I 'JOS  1211  I  l!Wtll|.  Interference  that  is  due  to  a  refracted 
wiivi  produces  oscillations  in  the  intensity  with  an  angular  quasi  period  of  magnitude  I  A/ri  I1  ’  rad  near  0, .  where 
,\  is  the  wavelength  and  11  is  the  radius.  Unlike  diffraction  related  os.  illations,  the  interference  oscillations  in¬ 
crease  in  magnitude  as  forward  scattering  is  approached.  < fptical  tunneling  through  bubbles  is  also  discussed. 


INTRODUCTION 

An  approximate  description  of  the  scattering  of  light  by  large 
dielectric  spheres  may  he  obtained  from  ray  optics,  provided 
that  the  observation  angle  is  not  in  a  region  in  which  wave¬ 
front  diffraction  is  dominant.’  In  the  ordinary  case,  in  which 
the  refractive  index  of  the  sphere  exceeds  that  of  the  sur¬ 
roundings,  forward  scattering,  the  rainbow,  and  the  glory  are 
the  regions  in  which  diffraction  is  important.1  :*  In  the  less 
explored  case  of  a  sphere  whose  refractive  index  is  less  than 
the  surroundings,  a  new  region  appears  in  which  diffraction 
is  important.’*  This  region  may  he  called  the  critical  scat¬ 
tering  region  because  it  is  associated  with  light  reflected  from 
the  sphere's  surface  when  its  angle  of  incidence  is  close  to  the 
critical  angle  normally  associated  with  total  reflection  from 
plane  surfaces.  Diffraction  effects  in  this  region  have  been 
observed  and  described  with  a  physical -optics  approximation. 1 
The  purpose  of  this  paper  is  to  discuss  the  effect  of  the  in¬ 
terference  of  the  reflected  wave  with  certain  rays  that  pene¬ 
trate  the  sphere  and  are  refracted  to  angles  in  the  vicinity  tit 
the  critical  scattering  region.  This  interference  leads  to  a 
modification  of  the  previous  model*  that  becomes  significant 
in  the  critical  region  when  an  air  bubble  in  water  has  a  radius 
<0.4  mm. 

I,et  the  scattering  angle  </>  denote  the  angle  of  the  distant 
observation  point  measured  with  respect  to  the  direction  of 
propagation  of  the  incident  plane  wave  and  the  sphere's 
center.  The  critical  scattering  angle  is  </i,  =  v  -  2 II, ,  when 
0,  =  arcsin  (n  ')  is  the  critical  angle  of  incidence  and  n  = 
(n„//ii)  >  l,/i„  and  n,  being,  respectively,  the  refractive  indices 
of  t  lie  outer  and  the  inner  media.  Air  hubbies  in  water  have 
n  ~  4/11  and  </>,  ~  82.82°.  Figure  1  illustrates  several  rays 
that  have  a  deviation  angle  0  of  50°,  Knelt  ray  is  characterized 
by  a  parameter  p,  where  p  -  1  is  the  number  of  reflections 
from  the  internal  surface  and  p  =  0  has  only  an  external  re¬ 
flection.  Uavs  that  enter  the  bubble  below  the  center  line  are 


indicated  with  a  prime.  The  angle  of  incidence  measured 
from  the  surface  normnl  is  It,,,  and  the  corresponding  refracted 
angle  is  Rays  with  p  ?  0  have  Up  <  II, . 

The  previous  physical-optics  approximation*  described  the 
diffraction  of  the  virtual  wave  front  associated  with  t  he  p  = 
0  ray  only.  It  was  shown  that  diffraction  is  important  in  the 
region  1 0  —  </>,.  |  (A/a  l1'-’ rads,  where  n  is  the  sphere  radius 
and  A  is  the  wavelength  in  the  outer  media;  consequently, 
diffraction  is  important  near</i,  even  if  A/a  is  quite  small.  The 
present  paper  shows  that  the  interference  that  is  due  to  the 
p  =  I  ray  leads  to  a  modulation  of  the  scattered  intensity  as 


0  ?  1 


Kig.  I  ltav  paths  In  the  scattering  plane  with  a  scattering  angle  0 
=  .ri0°,  The  number  gives  the  rav  parameter  />  The  flashed  ray  is 
forward  scattering  that  is  due  to  tunneling  Infill  the  reflected  ray  and 
is  discussed  in  Sec  I. 
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a  function  «l  «/>.  The  p  i  2  rays  will  lie  omitted  since  they  lead 
to  virtual  wave  fronts  much  smaller  in  amplitude  than  the  p 
~  d  >'nd  p  —  1  wave  fronts.  The  p  =  2'  ray  leads  to  a  super¬ 
posed  low-amplitude  modulation  of  the  scattered  intensity 
with  an  angular  period  ^A(-Vn),  where  a/ A  is  the  lateral 
separation  ot  rays  0  and  2’  after  emergence  from  the  sphere. 
For  an  air  bubble  in  water,  ray  tracing  gives  A  -v  0.82  with  (/> 
=  «/>, ,  whereas  observations  close  to  0 (Sec.  V  of  Ref.  4)  gave 
A  =  0.70  ±  0.04.  The  angular  period  of  this  modulation  is 
when  a  »  A.  The  emphasis  of  the  present  paper 
is  on  the  coarse  structure  in  the  scattering;  consequently, 
primed  ravs  are  omitted.  Mie- theory  results'1  (to  lie  published 
separately)  substantiate  these  approximations  in  the  region 
(i/i,  —  40°)  <  0  <  0,  for  n  —  4/0. 

1.  DESCRIPTION  OF  THE  VIRTUAL  WAVE 
FRONTS 

Diffraction  and  interference  effects  in  the  scattering  are 
modeled  by  first  using  ray  ootics  to  describe  virtual  wave 
fronts  positioned  within  the  si  atterer.  These  wave  fronts  are 
then  allowed  to  diffract  to  th*  fur  field.  Fortunately,  it  is  not 
necessary  to  express  the  exact  amplitude  and  phase  of  the 
wave  fronts  in  terms  of  Cartesian  coordinates.  This  is  because 
approximate  dependences,  dciived  below,  justify  the  use  of 
the  stationary-phase  upproxin  a. ion  in  the  diffraction  inte¬ 
grals. 

Figure  2  illustrates  Cartesian  coordinates  used  in  the  de¬ 
scription  of  virtual  wave  fronts  in  the  scattering  plane.  The 
origin  is  chosen  to  be  the  point  where  the  reflected  ray  with 
On  =  ft ,  contacts  the  sphere.  The  positive  e  axis  is  the  path 
of  the  scattered  critical  ray  predicted  by  ray  optics.  The 
positions  of  the  virtual  wave  fronts  associated  with  rays  0  and 
I  are  denoted  by  the  functions  l\,{u)  and  iq (u),  where  the  u 
axis  is  perpendicular  to  the  t1  axis  und  lies  in  the  scattering 
plane.  Phase  shifts  of  ray  0  that  are  due  to  reflection  beyond 
the  critical  angle  are  not  included  in  r„  but  are  accounted  for 
with  a  separate  phase  factor  in  the  description  of  the  ampli¬ 
tude.  The  reference  phase  is  selected  such  that  )>»(())  =  0. 
Ray  tracing  gives  in  =  iq  at  u  =  t),  and  the  leading  term  in  a 
Taylor  series  lor  in  is1 

I'n  =  <ui  (D 

where  o  =  -to  cos  0,  )  '.  The  Taylor  series  for  iq  may  be 
obtained  by  including  effects  of  refract  ion  and  t  he  change  in 
velocity  within  the  sphere;  the  leading  term  is 

i'i  =  'hiui'.  (2) 

» 

The  resulting  wave  fronts  are  shown  in  Fig.  2.  The  calculation 
ol  the  interference  of  waves  0  and  I  uses  the  exac*  path-length 
difference  that  is  not  easily  expressed  as  a  function  of  u. 
Effects  ot  wave-front  curvature  orthogonal  to  the  scattering 
plane  will  be  approximated  with  normalization  factors  to  be 
introduced  subsequently. 

In  the  description  of  the  scattering,  it  is  only  necessary  to 
consider  two  cases  with  orthogonal  incident  polarizations, 
since  the  general  case  may  lie  obtained  by  a  linear  combina¬ 
tion.1  The  subscript  /  is  assigned  the  value  I  when  the  inci¬ 
dent  electric  vector  is  entirely  perpendicular  to  the  scattering 
plane  and  the  value  2  when  entirely  parallel.  Waves  scattered 
to  the  tar  field  associated  with  p  =  0  and  p  =  I  rays  may  be 
simulated  by  a  wave  on  the  u  axis  with  the  following  ampli¬ 
tude: 


V 


Fig.  2.  Positions  of  the  virtual  wave  fronts  for  p  =  0  and  p  =  1 
si  uttering.  The  dashed  ray  is  the  incident  critical  ray  that  is  lellected 
to  become  the  e  axis.  A  portion  of  the  bubble's  surface  is  also 
shown. 


Up./Ut.t )  =  h,,j  exp|t(-  hvr  -  u)( )],  id) 

where  i  -  v  — I  ,  u>  is  the  frequency,  t  is  the  time,  und  h  = 
2 ir/\.  The  exp(-  iwt )  factor  will  he  omitted  subsequently. 
Ray  tracing  and  the  physical -optics  approximation  give  the 
following  amplitude  factors1-': 


h,yj  =■  r, ((),)),  u  <  0,  (4) 

h,).j  =  it  >  0,  (ft) 

/ii.,  =  0,  it  <  0,  (ti) 

h\j  =  [  1  —  rq((t|)”|  \  kl)  .  it  >  li,  (7) 


where  the  rf  are  Fresnel’s  amplitude  reflection  coefficients 
for  a  plane  surface,'1 


r  i(W,,)  =  : 


_  sin (/),,  -  Op) 
sintp,,  +  II,, ) 

r,(0p )  =  'JM&LZJhd  , 
tantp,,  +  (Ip) 

and  ftp  is  the  refraction  angle  predicted  by  Snell’s  law. 


18) 

(«) 


ftp  =  urcsintn  sin  ()„).  (Ill) 

The  Stokes  symmetry  relations  have  been  used  in  the  deri¬ 
vation  of  Kq.  (7),  where  I)  is  Van  de  Hoist's  divergence  factor,1 
which  may  be  written' 


li  _  _ sin  0 1  cos  0 1 

2|1  -  (n  cos  ff|/cos  P()|sin  0 

0  is  the  scattering  angle  oft  hep  =  1  ray  predicted  by  ray  <>p- 
t  ics, 


0  =  2(pi  -  t/|),  (12) 

and  iv  is  a  positive  constant  normalization  factor  determined 

subsequently. 

I  he  phase  shif tsM  ft)  are  associated  with  the  (nearly)  total 
reflection  of  ray  0  as  ft,,  exceeds  0,  (tunneling  through  large 
bubbles  can  be  neglected;  see  Sec.  4): 

tan((i,/2)  =  (sin- -  it '-I1  -/cost/,,.  lilt) 

tan(iVi/2)  =  n  -  tan(i),/2).  ,  pp 

I  he  ft,  vanish  when  0,\  =  ft, ,  but  they  become  positive  as  0„ 
exceeds  0, ,  which  represents  an  tuliaturmvtil  of  the  p  =  t) 
virtual  wave  front.  (The  reader  is  cautioned  that  some  an 
I  hors  write  ft,  and  A.  with  an  incorrect  sign  for  their  choice  of 
the  sign  ol  w  in  the  time-dependence  factor.11) 

Approximate  dependence  of  hpj  on  small  values  of  u  is 
obtained  by  first  relating  ft,,  ton  torn  ~  (land  ft,,  ~  ft, .  This 
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leads  to  the  following  proportionalities:  1  —  ht)J  a  (— u/h)i/- 
for  u  <  0,  ,'t  i  j  a-  (u/o)  ’  -  for  u  >  0,  and  1 

^  f)jU  u‘.  (If.) 

where  f1\  =  (8/nri)1-’  cos  f),  and  fi,  =  n-fi |. 

2.  FAR-FIELD  DIFFRACTION  AND 
INTERFERENCE  OF  THE  SCATTERED  WAVES 

Let  R  denote  the  distance  to  the  observation  point  from  the 
center  of  the  bubble.  When  R  »  kn-,  the  scattering  ampli¬ 
tude  is  proportional  to  the  Fourier  transform  of  U()j  +  f/| 
which  may  be  written  as  /,  =  fa,  4-  /|  ,,  where 

/,,,>(sin  ij)  =  J'  )  expt-  iku  sin  rjklu  (16) 

and  1)  =  0,.  —  0.  An  approximation  for  /o  ,  with  rj  >  0,  which 
was  introduced  in  lief.  4,  is  to  let  (/„ ,  =  0  for  u  <  0  because 
( hyj  decreases  abruptly  as  u  becomes  negative.  It  is  not  clear 
if  ray  optics  can  be  used  to  model  the  virtual  wave  in  the  region 
in  which  \dhuj/Au  |  is  large;  however,  its  use  leads  to  the  dif¬ 
fraction  of  reflected  energy  away  from  the  critical  region  and 
supports  the  approximation  U»,j[u  <  0)  -  '•  Most  of  the 
diffracted  p  =  0  wave  in  the  critical  region  comes  from  the 
virtual  wave,  where  u  is  small  and  positive,  and  an  approxi¬ 
mate  result  is  obtained  by  using  Eqs.  (1)  and  (15)  in  the  entire 
positive  u  domain.  The  integral  may  be  approximated  by 
using  the  principle  of  stationary  phase1"  with  the  result  that 
fnj  =  (A a  cos  0,.) 1/2 go, ;/2,  where4-1 1 

go  j  ^  expti'Yo., )  |F(u>)  -  F{-  «)|,  (17) 

w  =  [(a/A)cos  9, 1I/,J  sin  17,  (18) 

7n.>  -  - ft){w/2)l/1(\a  cmO,VM  +  — - •  (19) 

4(v 

wheic  the  first  term  of  Eq.  (19)  is  the  phase  shift  of  the  geo¬ 
metric  ray  reflected  to  <fi  and  l'(u  )  is  Fresnel's  integral, 

F(w)  =  exp[i(ff/2)z-|dr.  (20) 

When  the  observation  point  is  in  the  critical  region,  Kq  (2) 
may  be  used  in  the  evaluation  of  [\ y. 


wave.  Consequently,  the  error  from  using  Kq.  (2.'1)  in  the  total 
scattering  should  be  acceptable  even  if  rj  is  small. 

The  scattered  intensity  normalized  to  that  of  the  geometric 
result  for  a  perfectly  reflecting  sphere  of  the  same  size  is 

I]  =  ‘/z  U’oj  +  A’i,;|t  (24) 

where  the  normalization  was  obtained  by  noting  that  |a'oj  |  ► 

\  2  as  te  ■■»<».  'l'he  actual  intensity  is  t  he  incident  intensity 
multiplied  hv7  l,(a/R)'- 74.  The  value  of  s  in  Kq.  (7)  is  deter¬ 
mined  to  be  2  by  noting  that  only  then  does  at.;  become  the 
value  predicted  by  ray  optics.7  Because  of  the  interference 
of  the  ttpj,  the  intensity  is  sensitive  to  the  computation  of  (7oj 
-  71,).  Comparisons  with  Mie  theory  have  shown  that  the 
phase  difference  given  by  Hlqs.  (19)  and  (28)  is  limited  in  its 
usefulness  because  of  errors  introduced  hv  Eqs.  (1 ),  (2),  and 
(15).  The  agreement  with  Mie  theory  is  significantly  im¬ 
proved  when  (y„j  —  71 ,)  is  taken  to  he  the  exact  value  pre¬ 
dicted  by  ray  optics17: 

7o  ,j  ~  7i  j-  2/m(cos()|  -  cosflo  -  n  ~ 1  cos  /q )  -  ?i,  (()„), 

(25) 

where  On  is  predicted  bv  ray  optics  to  be 

<>»  =  (ir-  0)/2  *  0,  +  (i?/2).  (26) 

To  evaluate  Eqs.  (22)  and  (25)  for  a  given  value  of  0,  it  is 
necessary  to  solv'  the  transcendental  Eqs.  (10)  and  (12)  for 
0 1  and  pi.  The  evaluation  of  Kq.  (24)  requires  the  evaluation 
of  Eqs.  (7)- (14),  (17),  (18),  (22),  (25),  and  (26). 


3.  DISCUSSION 

Figures  8-6  illustrate  the  result  of  three  approximations  for 
the  scattered  intensity  for  an  air  bubble  in  water  normalized 
to  the  ray-optics  value  for  a  perfectly  reflecting  sphere.  In 
each  figure,  n  =  4/8.  The  ray-optics  result  for  a  bubble  is 
equal  to  1  +  2.*i  i,/-’ when  tj  >  Oand  ry  (On)- when  17  <  0.  It  ex¬ 
hibits  a  cusp  at  the  critical  angle,  which  is  also  evident  in  the 
rav-optics  predictions  plotted  by  Davis7  and  by  Welford.1-’ 
'I’he  predictions  for  the  diffraction-only  model  developed  in 
Kef.  4  are  given  by  |g,> ,  1 1/2.  The  present  model  must  even¬ 
tually  break  down  in  the  region  where  17  <  0  because  of  the 


fij  ~  h  1,,  exp|-  i\ku  sin  17  +  (Auru<‘72)j|du.  (21) 

The  phase  of  the  complex  exponential  is  stationary  when  u 
=  ii,  where  17  =  a  cos  0,  sin  r).  Since  h  \ ,  is  a  slowly  varying 
function  of  u,  the  stationary-phase  approximation  may  he 
used  to  evaluate  Eq.  (21)  provided  that  U  is  well  within  the 
domain  of  integration.1"  The  result  is  f\ ,  =  (Aa  cos 
gl  .,/8,  where 


AT.;  =  2/i  i.,(f),)  exp 

1(7  +  71;)  +0(A/n)'6-\ 

1  \4  ) 

(22) 

Ti.;  = 

A  (sin  »/)- 

2« 

(28) 

In  Eq.  (22)  the  stationary-phase  point  that  is  used  in  the 
evaluation  of  h\ s  corresponds  to  the  f)  1  given  by  Eqs.  (10)  and 
( 1 2).  The  use  of  the  stationary -phase  approximation  breaks 
down  as  r)  -►  0,  since  Ti  -*  0;  however,  this  is  just  where  h  * 
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Fig.  8.  Three  approxiinntions  for  scattered  intensity  is  a  (unction 
ol  1/  —  '/'  ~  1 />,.  with  11  =  4/,'!,  the  electric  field  perpendicular  to  the 
scattering  plane,  and  An  =  10,000.  The  dashed  curve  is  the  result  of 
rav  optics,  which  adds  the  intensities  of/;  =  0  and  /;  =  1  scattering, 
'l'he  dotted  curve  is  the  diffracted  /;  =  0  wave  only  and  is  from  Kq. 
I  IT).  The  solid  curve  is  given  by  Kq.  (241,  which  includes  the  dll 
traction  of  the  p  =  0  wave  and  the  interference  of  the  /)  =  1  wave. 
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Kit;.  4.  Same  ns  Kit;,  4  lull  with  parallel  polarization. 


‘  ,!o  to 

r,  OF  GUFFS 

Kit;.  •>.  Same  as  Kit;.  ■'  hot  .  it  h  parallel  polarization, 

assumption  that  ( (u  <  (5  =  0.  In  the  figures,  the  morlel 
hits  been  extended  into  that  region  by  evaluating  Kip  (24)  with 
g,  ,(t]  <  0)  =  0.  Fresnel's  integral  F(icl  was  evaluated  with 
it  nunierieal  approximation.1 1 

The  bubble  sizes  in  these  figures  are  manifest  by  consul 
ering  illumination  by  a  He  Ne  laser  with  a  wavelength  in  air 
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of  (>52.8  nm  and  a  wavelength  in  water  A  ~  474.(1  tint.  Then 
ha  =  10  000  gives  a  =  0.755  mm.  When  measuring  the  scat¬ 
tering,  refraction  corrections  may  be  required  to  relate  the 
scattering  angle  in  the  water  to  that  seen  by  an  external  in 
st  rument.4 

The  diffraction-only  theory4  predicted  oscillations  in  the 
intensity  that  are  homologous  with  respect  to  bubble  size, 
polarization,  and  refractive  index.  They  decrease  mono- 
tonically  in  magnitude  as  r/  increases  and  have  an  angular 
quasi-period  £(A/n)1/2  rad.  Figure  5  shows  that,  for  j  =  1  and 
ha  =  10,000,  interference  alters  the  scattering  significantly 
for  >/  2  10°.  After  the  fourth  maximum,  the  oscillations  in¬ 
crease  in  magnitude.  It  is  noteworthy  that  the  quasi-period 
of  the  newly  predicted  oscillation  is  also  roughly  (A/n  )l/i!  rad 
Jits  may  he  shown  by  Kqs.  (1)  and  (2)|.  Figures  4  (1  show  that 
the  effect  of  interference  is  even  larger  for  j  =  ‘2  and  lor  ha  ~ 
1000.  'The  intensity  averaged  over  a  quasi-period  is  roughly 
the  ray -optics  result.  Consequently,  for  large  spherical 
bubbles  in  white  light,  the  ray-optics  results  are  useful  except 
a  lien  ic  <  2,f>,  where  diffraction  corrections  become  essential 
(ic  =  2.5  corresponds  to  tj  =  4.4°  at  ha  =  10, (KM)  and  r)  =  14.0° 
at  ha  =  1000.) 

The  experimental  data  plotted  in  Ref.  4  were  for  ha  *  10 
950  and  j  =  2.  A  comparison  of  Fig.  4  with  those  data  indi¬ 
cates  that  it  is  not  surprising  that  the  diffraction-only  model 
correctly  described  the  relative  positions  of  the  first  four 
maxima  but  that  it  underestimated  the  magnitude  of  the  in¬ 
tensity  oscillations  for  tj  >  5°.  A  new  apparatus  has  t)een  built 
that  permits  observations  in  the  region  n  >  7°,  which  was 
obscured  in  the  original4  apparatus.  Visual  observations  for 
air  bubbles  with  a  -v  1  mm  in  monochromatic  light  suggest 
that  the  magnitude  of  the  intensity  oscillations  tends  to  in¬ 
crease  with  increasing  rj  but,  as  of  the  time  of  this  writing, 
photographs  of  a  quality  to  permit  quantitative  measurements 
have  not  been  obtained.  These  observations  preceded  and 
motivated  the  present  theoretical  effort. 

The  intensities  given  by  Kq.  (24)  have  been  compared  with 
Mie  theory1  14  with  ti  =4/11  and  ha  =  25,  100,  1000,  5000,  and 
1 0,000.  The  computations  were  performed  using  Wiscombe’s 
Mie-scattering  algorithms1'’ with  minor  modifications.  The 
comparision  shows'’  that,  for  0  <  rj  <  110°,  the  course  structure 
is  approximately  described  by  Kq.  (24).  Krrors  in  the  pre¬ 
dicted  locations  of  intensity  oscillations  are  typically  less  than 
one  tenth  of  a  quasi-period.  Furthermore,  the  phase  differ¬ 
ence  of  the  scattering  amplitudes  for  the  two  polarizations 
given  by  Kq.  (28)  (below)  reproduces  the  coarse  structure  in 
the  Mie  results.  Kquation  (24)  differs  most  from  the  coarse 
structure  in  the  Mie  results  when  j  =  2.  This  may  be  due  in 
part  to  its  omission  of  the  p  =  2  wave,  which  has  a  smaller 
amplitude  for  j  =  1  than  for  j  =  2.  F'or  example,  at  >/  =  50° 
the  intensities  of/)  =  2,  scattering  in  the  units  of  Figs.  5  (5,  arc 
given  by  ray  optics'  to  be  0.012  and  0.015,  respectively. 

Kquat  ion  (24)  does  not  describe  fine  struct  are  in  the  Mie 
intensities  which  is  most  prominent  when  j  -  1.  This  si  rue 
lure  has  a  quasi -period  ~A(A/n)rad  (0.5°  for  ha  ~  1000), 
which  indicates  that  it  could  be  included  in  the  model  with  the 
addition  of/)  =  2’  scattering.  Ray  optics  verifies  that  the 
amplitude  of  2'  scattering  in  the  critical  region  is  largest  when 
/  =  1.  Kxcept  for  the  line  structure,  thd1  Mie  intensities  dc 
crease  gradually  as  r;  decreases  through  zero.  The  error  in  Kq. 
(21),  with  )/  negative,  generally  decreases  with  ha  as  diffract  ion 
becomes  more  important . 1 
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Equation  (24)  may  In-  converted  to  the  units  commonly  used 
in  Mie  algorithms  through  the  following  transformation: 

|.S-, !  ■  =  (A-nt-7,/1,  (27) 

where  .S’,  are  the  scattering  amplitudes  in  the  usual  units.11''111 
The  phase  difference  of  .S'  |  and  .S',,  is  needed  if  the  scat  ter  inn 
for  iirhitrary  states  of  the  incident  polarization  is  to  he  cal¬ 
culated. 1111  The  phase  difference  is  predicted  to  he1 ' 

arg(.N\i)  -  urg(N|)  =  argtgo..,  +  am..')  -  argU’iu  +  am  t). 

(28) 

where  arg  denotes  the  complex  argument  of  the  quantity  in¬ 
dicated.  The  present  calculation  uses  an  exp(—  iu)t )  time 
factor,  whereas  the  .S',  are  often,11’1  hut  not  always,11’ specified 
by  using  exp(  +  iu>( ),  so  the  right-hand  side  of  Kip  (28)  should 
he  multiplied  hy  —  1  for  some  applicat ions. 

The  coarse  structure  in  the  scattering,  which  arises  because 
of  the  interference  of  the  reflected  and  refracted  light,  may 
he  a  useful  experimental  tool.  The  angular  position  of  the 
maxima  should  shift  if  there  are  changes  in  the  size  of  a  bubble 
or,  for  the  case  of  tt  drop  with  n,  <  n„,  if  the  refractive-index 
ratio  changes.  Shifts  in  the  structure  of  rainbow  scattering 
have  been  useful  for  detecting  shape  oscillations  in  drops, 1H 
and  a  similar  technique  could  he  used  for  detecting  small 
changes  in  the  scattering  from  hubbies. 

Diffraction  associated  with  the  critical  angle  and  structured 
scattering  that  is  due  to  interference  should  also  he  present 
in  the  scattering  hy  cylindrical  bubbles  or  dielectric  cylinders 
with  n,  <  n„.  for  a  circular  cylinder  with  a  symmetry  axis 
perpendicular  to  the  incident  beam,  the  critical  scattering 
angle  and  ray  tracing  in  the  sentterinp  plane  are  the  same  as 
those  described  here  for  spheres.  The  divergence  factors  used 
in  the  description  of  the  virtual  wave  fronts  differ  from  those 
of  a  sphere.  The  geometric  optics  of  a  large,  perfectly  re¬ 
flecting  cylinder11'  show  that  (neglecting  tunneling)  |hn  ,(n  > 
t))|  is  no  longer  constant  but  is  "  |sin(«/»/‘2)| 1  The  normalized 
magnitude  and  positions  of  the  intensity  oscillations  should 
differ  from  those  of  a  sphere.  Critical  angle  diffraction  and 
interference  structure  should  also  he  present  in  acoustic 
scattering  from  fluid  spheres  and  cylinders. 

4.  ELECTROMAGNETIC  TUNNELING 
THROUGH  SPHERICAL  BUBBLES 

The  present  model  hits  assumed  that  the  reflection  is  total 
when  II  n  >  II, .  The  existence  of  evanescent  waves  within  t  he 
bubble  wiil.  however,  frustrate  the  total  reflection,  as  in  the 
case  of  plane  dielectric  slabs  separated  by  an  air  gap.  M  The 
effect  of  tunneling  on  nearly  total  reflection  at  curved  di¬ 
elect  ric  interfaces  has  been  studied  to  determine  losses  from 
curved  optical  waveguides”;  however,  the  sense  of  the  cur 
suture  differs  from  that  in  the  present  case.  The  effect  of 
tunneling  through  bubbles  can  be  estimated  by  using  the 
plane  slab  results.  1  in  which  the  gap  <1  is  equal  to  I  lie  bubble 
diameter  and  H„  becomes  the  angle  ol  incidence  at  t  lie  plane 
surface.  The  result  of  this  approximation  is  that  the  wave 
that  has  tunneled  through  the  bubble  will  always  lie  trans¬ 
mitted  in  the  forward  direction.  The  dashed  ray  in  h’ig.  I 
shows  the  rav  transmitted  because  of  the  trust  rated  reflection 
ol  I  lie  / i  =  ll  rav.  For  each  /)  =  II  ray  with  II,,  >  H, .  t here  is  a 
forward  transmitted  ray  with  r  =  II,,.  The  intensity  trails 
mil  lance  I  h rough  a  plane  gap  is  found  to  lie  of  order  (/.•</! 
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when  II, i  =  II,,  and  it  decreases  rapidly  as  II,,  exceeds  II, 
Tunneling  should  have  a  negligible  effect  on  the  reflection 
from  a  bubble  unless  ka  is  small.  This  conclusion  is  also 
supported  by  the  success  ol  the  present  model  in  reproducing 
Mio  results'’  when  ku  =  2,r>. 
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Mie  scattering  algorithms  were  used  In  Compute  scattered  intensities  niui  phase  differences  lur  air  iiutitiles  in 
water.  Results  are  plotted  as  a  I'mictinii  ol  tile  srattennit  angle  </>  in  the  general  range  <>l  do  UUU  lur  si/e  parameters 
l,o  nl  do.  Hill,  l.ooo.  and  10.000  leurrespondini:  lo  radii  ii  -•  1  .d  pm  to  o.h  non I.  As  v  di'ereases  lielow  tite  eritieal 
seat  ter  ini;  angle  at  H'J.S",  I  lie  intensity  increases  and  undergoes  broadly  spaced  nsci  Hat  ions  that  are  described  by 
a  physical -optics  approximation  developed  liv  the  authors  in  a  separate  publication.  Mie  scattering  also  exhibits 
finely  spaced  oscillat ions 


INTRODUCTION 


( Ihserval  ions  of  light  scat  I  eretl  by  air  I  at  I  tides  in  water  in  I  lie 
vicinity  of  tile  crit  inti  scattering  angle1  indicate  t  tint  dilTrat 
t  ion  is  significant,  even  if  I  lie  lull  title  rad  ills  a  is  as  large  as  I 
tit  in.  The  coarse  st  met  tire  in  I  lie  critical  scattering  region  lias 
been  described  with  a  physical-optics  approximation  that 
accounts  for  I  lie  (lit traction  of  light  reflected  front  t  lie  bubble's 
surface  and  t  lie  interference  of  certain  refracted  light/’  In  t  bis 
Letter  we  use  Mie  theory  1  •’  to  calculate  far  field  scattered 
intensities  and  phase  differences  in  the  critical  region  of 
spherical  bubbles  in  water  and  compare  these  exact  results 
wit  It  t  hose  of  t  be  physical  optics  approximation. 1  Previous 
applications  of  Mie  theory  to  bubbles  emphasized  t  be  ealcti 
lat  ion  of  the  scattering  efficiency  factor'1  *  and  t  lie  radial  ion 
pressure. M  In  contrast  with  the  scattering  by  spherical 
drops,1"111  the  exact  angular  structure  ir  the  scattering  by 
bubbles  is  not  wei!  explored. 

Consideration  of  the  approximate  scattering'’  introduces 
the  features  to  be  expected  in  I  he  exact  Mie  results  and  la 
cilitutcs  the  choice  of  the  range  and  frequency  of  the  angles 
needed  for  a  descript  ion  of  t  he  angular  st  rticlure.  figure  I 
i  I  lust  rales  several  rays  with  a  scat  tering  angle  i/m  tl  50°.  In  t  his 
I  igure  and  the  Mie  eonipulat  ions  I  hat  follow,  the  relative  re 

I  ract  ivc  index  it,  In,,  is  taken  to  be  li/l.  where  a,,  is  I  lie  refi  ne 
live  index  of  the  water  and  n,  is  the  refractive  index  of  the 
bubble  contents.  Rays  tire  characterized  by  a  parameter  />. 
where  />  —  I  is  t  he  number  of  reflections  within  t  he  buhl  tie  and 
/>  =  It  lias  only  an  external  reflection;  I is  the  angle  ol  inci 
deuce  at  t  he  /it  li  ray  when  it  I  irst  touches  t  he  surface  of  I  he 
bubble.  When  sin  II  »  //,  hi,,,  t  he  /)  =  0  ray  is  totally  reflected 

I I  In  effects  of  tunneling"  and  surface  curvature  are  neglected  I 

The  corresponding  condition  on  the  scattering  angle  ol  the 
reflected  ray  is  i  ■  -  HIS'T.  where  i,«, .  the  critical  scat 

ter  mg  angle,  is  related  to  n,/n„  through  Kqs.  I  I )  and  CJI  ol  Ret. 
I .  I  for  air  bubbles  in  seawater.  •;»,  may  be  as  large  as  M  I". I 
I  lie  physical  opt  ics  approx i mat  ions1  lead  to  t  lie  predict  ions 
ol  st  met  me  in  t  he  angular  scat  I  ering  wit  h  an  angular -spread 


of  t lie  order  of  < Mi) 1  rad  |  ~  I  TKA'tt  I  1  ■’ deg |.  where  A  is  the 
wavelength  t/t  the  miter  medium  and  I;  =  'Jir/A.  Structure 
with  a  quasi -period  nl  that  magnitude  will  be  referred  tons 
roar  si’  si  meltin',  in  coni  rast  to  /ine  strut  lure  in  the  scattering 
that  has  a  quasi  period  in  the  crit  ical  region  of  less  limn  A/ri 
rad  (or  IM50 °/hn  I.  The  f  ine  st  rttcl  u re  original es  in  part  from 
the  interference  of  tlie/i  =•  2  ray  with  the  rest  nl  the  scattering. 
To  describe  the  general  effects  of  the  critical  angle  on  I  lie 
scattering.  Mie  com  put  at  ions  are  needed  wit  h  a  spread  ol  |  •/<, 
-  </>|  somewhat  larger  than  (Ml  I1  with  angular  steps  some 
what  smaller  than  A  Mi  to  prevent  sampling  errors 


COMPUTER-PROGRAM  CONSIDERATIONS 


The  Mie  solutions  were  computed  by  using  the  I-'OHTHAN 
MIKVO  subroutine  developed  by  Wiscninbc,  ‘  1 1  chosen  in 


l‘  ij:  I  Kay  pat  hs  in  tin-  >cn1 1 1'lini;  pin  lie  wil  h  n  scjilteriii^  alible 
■><»'.  I  lie  number  adjacent  tn each  ray  uivr'-  !he  rav  parameter  /> 

I  lie  appmxi  iii.il  uni  j  *  i  \  e  n  in  Kef.  :f  UK.  -hide's  the  ini  erf  etenre  and 
Hit  t  rat  !  i«  m  "I  wau  ^  ;t  ssnrial  «•«  |  with  ravs  o  and  I 
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preference  In  his  Ins) cr  MIKV  I  mill  inc  because  ol  I  he  menu  irv 
limiliit inns  of  i ho  HI’  1000  computer  system  used  in  the  work. 
Wo  found  it  necessary  lo  convert  nil  VVisootnho's  inlornnl 
vnrinlilos  lo  III1  doiihlo  prooision  1 1  I  significant  digits).  nod 
wo  wore  I  lion  aide  In  duplicate  his  'utiles1 1  for  seal  I  on  ng  liv 
a  drop  with  n, In,.  =  I  d.  iWiscomhe’s  siiliruulines  assume 
I 'I -digit  single  precision.  I  Must  ci implex  arithmetic  was  re 
written  as  real  arithinet  to  I  mm  also  of  t  ho  absence  of  con  ;>lo\ 
dotthlo  precision.  The  logarithmic  derivative  of  the  III 
cut t i  Bessel  function,  ootnnionlv  denoted  liv  A„.  wascahn 
luted  hy  downward  recursion,  which  is  known  to  ho  stable. 
The  program  was  tested  against  Wisoomhe's  tables  for  /i,/n,. 
=  l.a  of  the  complex  scattering  amplitudes  and  effiibucv 
factors  and  against  the  scattering  elTiciem  y  factors  tor  sma.'1 
size  parameters  when  u,/n„  =  It.T.d  as  ;*  von  in  Refs.  7  and  K. 

Lot  N  denote  the  mini  her  of  terms  included  in  the  cotnpu 
tat  ion  of  the  Mie  series  and  /Vimii  denote  the  number  of  angles 
for  which  the  Mie  scattering  is  computed.  Since  MIKVtt 
uses'’’'1 1  N  =  f<u  +  1. Dot  An )'  -'  +  2  and  N.,„„  must  also  lie  in 
creased  in  |iroportion  to  An  to  resolve  the  line  strut  o,  exe¬ 
cution  time  is  roughly  proportional  to  I  An ) '  for  a  fixed  spread 
of  angles,  Whereas  the  results  for  An  =  2d  required  a  fraction 
ot  a  minute,  those  for  An  =  ID, HDD.  whole  the  angle  step  size 
was  D.DDd0,  required  2.  a  h  for  HDDD  angles.  I  These  times  were 
realized  after  veetori/.ation  of  two  loops,  which  speeded 
u))  execution  hy  a  factor  of  approximately  I.) 

Higher  values  of  ,V  than  that  given  in  Wiseonihe's  etpialion 
above  were  tested  for  the  bubble  ease  with  An  =  1  DO  and  IDDD 
and  found  to  yield  identical  results  to  six  decimal  places. 

In  contrast,  the  physical -optics approximation  program  was 
executed  in  seconds,  regardless  of  size  parameter.  K<|iiatious 
used  in  that  program  are  listed  at  the  end  of  Section  II  of  Hof. 


KKSUr.'IS 

( 'nictitated  intensities  /,  as  a  fund  ion  of  t  tie  scattering  angle 

are  shown  in  Figs.  2  I.  The  normalization  was  chosen  so 
that  / , t >;> f  =  I  represents  perfect  reflection  according  to  geo 
metric  optics,  '1 '  At  a  distance  A’  »  Au  '  from  the  center  of 
the  sphere,  the  actual  /  polarized  ittlcnsilv  is  the  incident 
I  polarized  intensity  nniltiplied  liv  /, tu//i  I  7  1,  where  /  =  I 
denotes  polarization  of  the  electric  vector  perpendicular  to 
t  he  scat  tori ng  plane  and  /  =  2  the  parallel  case.  The  van  tie 
1 1 u  1st  normalization1  used  in  Wiseonihe's  routines  retptiretl 
that  the  squared  modulus  of  the  scattering  amplitude  com 
pilled  with  MIKVO  be  multiplied  hy  (2/Aril  '  to  be  expressed 
in  I  he  units  of  Figs.  2  1. 1 

As  is  shown  in  the  figures,  t he  exact  solution  suggests  that 
a  fine  ripple  si  met  tire  is  superposed  upon  the  coarse  structure 
that  is  described  by  the  physical  optics  approximation.  The 
ripples  have  a  quasi  period  of  about  D.H2  l  >\/<;  I  rad  near  i/i,  . 
decreasing  slightly  with  To  avoid  severe  distortion  re¬ 
sulting  from  sampling  errors,  it  was  necessary  to  use  an  angle 
step  size  less  than  ,\/7u  rad.  This  fine  structure  is  dti"  to  t he 
interference  of  ravs  not  included  in  the  approximation,  pri 
ina r i I x  the  p  =  2’  ray.  Consequently,  the  ripples  should  he 
damped  out  for  bubbles  containing  a  gas  that  absorbs  light. 
Fine  structure  from  air  bubbles  has  been  observed.1 

The  plots  in  Fig.  a  show  !  lie  phase  difference  of  the  scat 
I  ered  waves 


(b) 


log  '.!.  Calculated  normalized  scattering  intensities  lor  An  =  ID, (Mm. 
The  electric  vector  is  parallel  to  tile  seatiering  plane.  The  solid  curve 
is  Ironi  Mie  lla-oiA  The  dash  curve  is  the  phvsieal  optics  ap 
proximal  ion  go en  in  Ret  2. 

a  =  nrgiN  d  -  argt.S' 1 1,  ill 

win  re  the  S',  ate  Mie  scattering  amplitudes  for  a  time  do 
pendence  of  expt  ~i u :l  I.  where  w  is  t he  frequency  and  I  t he 
t  imc.  This  phase  difference  can  he  used  along  with  / 1  and  / 
to  compute  the  scattered  intensity  and  polarization  for  arbi 
trarv  states  of  incident  polarization, 4  The  physical  optics 
approximation  lor  o  is  given  by  Kq.  I2NI  of  Ref.  2.  As  is  de¬ 
scribed  in  Ref.  2.  it  was  necessary  to  reverse  the  sign  of  a  from 
the  output  of  t he  MIKVO  program  because  Wiscombc  assumes 
a  time  dependence  of  expt  t  iw't  I. 


DISCUSSION 

For  ■;>  >  <;>, .  the  agreement  with  the  approx  hunt  ion  is  best  at 
low  size  parameters  where  diffraction  accounts  for  the  seat 
taring  in  this  region  in  the  approximation. 

The  coarse  structure  oscillations  at  angles  >/>  <  «/•,  are  similar 
to  the  predictions  of  diffraction  theory,  which  gives  a  rellecled 
intensity  proportional  to  t he  squared  modulus  ol  Fresnel  in 
legrals  |K,q.  DID)  of  Ref.  l|.  The  first  maximum  of  the  re 
fleeted  intensity  occurs  at 1 

ip,  -  y.  -  1 .21  \/u  cos  II,  l 1  1  t rad t  (21 

The  Cornu  spiral  behavior  of  the  Fresnel  integrals  predicts 
a  decrease  in  the  amplitude  of  I  la  se  intensity  oscillations  with 
decreasing  <:>  and  also  a  decrease  in  their  quasi  period,  The 
amplitude  decrease  is  most  evident  at  huge  size  parameters 
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Ki(i,  It.  I, ike  Kitt.  2  Imt  with  (n)  A a  ■  HMK)  and  the  electric  vector 
perpendictilar  to  scattering  plane  ami  tl>)  An  ■»  1(M)  and  the  electric 
vector  parallel  to  the  scattering  plane. 

(Fig.  2)  where  the  diffraction  pattern  is  less  spread  out;  at 
small  size  parameters  the  diffraction  decrease  is  lost  because 
of  the  increasing  amplitude  of  the  intensity  oscillations  as¬ 
sociated  with  the  interference  of  the  p  =  1  ray.  These  oscil¬ 
lations  are  also  of  quasi-period  <(A/n)l/2  rad,2  diminishing 
with  </>•  In  addition  to  the  j  polarizations  of  Figs.  2  4,  the 
physical-optics  approximation  has  been  verified  lor  the  j  not 
shown  and  for  both  j  at  An  =  5000.  Both  the  approximation 
and  the  Mie  theory  predict  that  the  positions  of  the  coarse 
maxima  depend  only  weakly  on  j;  the  coarse  maxima  lie 
slightly  closer  to  </y  when  j  =  I.  The  approximation  predicts 
that  the  intensity  variations  of  the  course-structure  oscilla¬ 
tions  are  greatest  when  j  -  2,  which  agrees  with  Mie 
theory. 

Mie  theory  shows  that  the  finc-ulructurr  intensity  varia¬ 
tions  are  larger  for  j  =  I  than  for;  =  2;  this  is  probably  because 
p  *  2'  scattering  is  greatest  when  j  =  I  For  lit)0  <  •/>  <  90°, 
‘.hi  positions  of  the  fine-structure  intensity  oscillations  for 
j  =  2  are  typically  shifted  from  those  for  j  =  1  by  one  half  of 
a  fine -struct ure  quasi-period,  except  when  </>  is  slightly  less 
than  ■/>,.. 

It  is  evident  from  Fig,  5  that  the  coarse  structure  of  the 
phase  difference  for  Mie  scattering  roughly  follows  that  of  the 
physical-optics  approximation.  This  has  also  been  verified 
at  ka  =  25.  The  gently  sloping  curve  of  the  approximate  re¬ 
sult  near  </>,-  may  he  understood  in  terms  of  the  phase  advances 
bj  of  the  (totally  reflected)  p  =  0  ray;  &■>  is  approximately 
(r»„/ri, )2^!  near  </>,  . 1,2  Comparison  of  Figs.  5  and  9  shows  that 
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the  maxima  in  the  coarse  structure  ol  |A|  are  associated  with 
minima  in  the  coarse  structure  of  the  /,  An  anomaly  is 
present  in  Fig.  5(h)  for  •/>  ~  97°.  Mie  computations  of  /,  in 
dicate  that  l\  is  quite  small  at  this  anomaly;  for  Am  =  10(1,  lt 
has  a  minimum  value  of  H  X  to  4  in  the  units  of  Fig.  9(h)  at  <l> 
=  911.705°.  The  | dot  of  A  for  ka  ~  25  reveals  a  similar  anomaly 
at  i />  ~  29°. 

The  range  of  bubble  sizes  associated  wit  h  Figs.  2  5  will  be 
illustrated  by  considering  green  illumination  with  a  wave 
length  in  air  of  550  mu  and  a  wavelength  in  water  A  of  4 19  mu; 
then  Aa  -  10,000  gi\  --s  n  =  (552  pm  and  An  =  25  gives  n  =  l.ll 
pm. 

It  is  remarkable  that  the  physical-optics  approximation 
describes  the  general  features  of  the  broad  intensity  maxima 
at  </>  ~  95°  for  An  =  25.  A  stationary  phase  approximation 
was  used  in  Bel's.  1  and  2  toeliminate  a  phase  term  |i /,(■'.  u  ) 
in  F,q.  <24)  of  Kef.  !|  from  the  diffraction  integral.  The 
arguments  given  in  Kef.  I  that  justify  that  approximation 
would  appear  to  break  down  for  An  -  25.  The  rough  success 
of  the  model  indicates  that  the  maximum  can  lie  attributed 
in  part  to  t he  first  diffraction  maximum  and  in  part  to  the 
interference  of  the  refracted  wav  associated  with  the />  =  1 
ray;  however,  25  may  be  close  ,o  the  lowest  value  of  An  for 
which  this  model  is  useful. 

The  principal  purpose  for  developing  and  testing  the 
physical-optics  approximation2  was  to  see  if  the  coarse 
structure  in  the  scattering  could  be  related  to  the  dilTraetion 
and  interference  of  reflected  and  refracted  waves,  The  model 
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scattering  ANGLE 


Kit:  I  Like  Fig.  2  Imt  with  Ail  =■  29  nail  the  eleelrii  vector  (a)  per 
penilii  ul.-ir  nail  (til  parallel. 
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.’>  I  ‘I  lase  difference  tor  the  siiil  tcriiii:  amplitudes  in  (lie  Hvo 
pulariznllmi  slates  Ini  In)  Ini  -  1 1  MM>  mid  ill)  An  =  |ihi.  Thesnlid  curve 
i'  I  ri  >ni  Mic  theory  mu'  Ki|  i  1 1.  The  dashed  nine  is  I  lie  plivsn  ul 
.iptii's  n|>|>r< >x i mm m >11  hv  I'li |  i'.IHi  ul  Hi'i,  :: 


mnv  Imvi'  praelienl  applications,  however.  even  though  il  omils 
till'  line  structure  in  I  lie  scnltcriug.  For  «,\imi|*li'.  in  the 
scattering  hv  a  splu'l'ical  polvdispcrsion  id'  liulililcs.  the  line 
structure  would  be  lust  anil  llic  course  structure  retained  il  the 
spread  ul  liulilile  sizes  nr  wavelengths  was  mil  inn  large, 
Furthermore,  I  lie  iiindel  limy  lie  ilselul  in  I  lie  d'-sign  of  in 
st ruments  I Imt  ih  pend  on  the  seat tcring  ol  light  by  I  in  Id  lies. 
Kxainples  ineliide  lal  deviees  Inr  in  ■■ iln  optical  ineasiiremeiil 
nl  inicrnlnililile  populations  al  sea;  dil  laser  Doppler  ane 
ninnielry  nl'  I  mill  lies  in  I  ii  pi  ids;  and  lei  I  lie  del  eel  inn  nl  liulilile 
illipi'i'lect  inns  in  idass.  .Inlnisnn  and  Cooke1”  have  used  (al 
In  i  ill!  ain  evidence  nl  a  lack  nl  I  mill  lies  at  sea  wit  Ii  radii  nl  less 
Ilian  approximately  .III  pin.  in  apparent  disagreement  with 
similar  aenuslic  nieasiireuients.  Their  instrument  made  use 
nl  scattering  wil  h  */.  ~  !Hi°,  which  is  shown  by  l'’i|;s.  H  and  1  in 
lie  much  smaller  Inr  mielnlnilililes  I  hall  Inr  a  geometric  re 
Heel  ion  front  a  perlecl  rel'leclor.  Moth  the  physical  optics 
approximation  and  the  Mic  I  henry  suggest  that  liulilile  <le 
I  eel  ion  insl  ruineiils  reiving  on  liislal  ic  seal  lering  should  ac 
cept  light  with  less  I  han  I  hal  of  I  he  first  dill racl  ion  maxi 
mum  as  predicted  by  Kip  I'd).  ( ienniel  l  ie  npl  ies1 "  gives  /,  > 


I  Inr  all  i/i, ,  conseipiently,  il  significantly  overeslimaies 
I  he  scattered  intensity  close  In  •/»,  when  hn  <  1000. 

Keller’ 111  has  attempted  In  use  light  scattering  with  </»  ~ 

I  III1,  In  measure  I  he  size  spec!  rum  nl  niiernluihliles  in  water  for 
the  purpose  nl  testing  models  of  cavitation.  Deferences  11 
and  In  give  angular  scattering  patterns  (reported  to  he  ob¬ 
tained  from  Mie  theory!  with  A  ~  0,(i:i2H  pm/l  .Till  and  a  = 

I .  a,  7.a,  and  10  pm  (which  correspond  In  l;u  ~~  1 0.  (ill,  100,  and 
Iddl.  Comparison  of  the  results  nl  our  Mie  program  and 
madel  show  that  Keller's  patterns  err  significantly;  they  do 
not  include  line  structure  and  they  incorrectly  describe  the 
coarse  si  ruclure. 

We  are  grateful  to  W.  .1.  W’iseombe  for  providing  the  iuitita! 
compuler  program  from  which  the  program  used  here  was 
derived  and  to  II.  Medwin  for  bringing  I  be  work  of  Johnson 
and  Cooke  loour  attention.  This  research  was  supported  in 
part  by  I  be  Washington  Stale  l  'Diversity  Research  and  Arts 
Cum  mil  tee  and  by  I  be  ( )|T  ice  of  Naval  Research.  1 1  1  .Mar 
ston  is  all  Alfred  I*.  Sloan  Research  Follow. 

‘  A u I  hor  to  whom  correspondence  should  be  addressed. 

KF.FF.RF.NCFS 

I .  P.  I ..  Mar-Ion.  "I  'lit  ion  I  angle  seal  I  ering  by  a  liulilile:  |>li  v-u-nl 
optics  approximat ion  and  observations."  •!.  dpi.  See.  Am.  (ill, 

I  do;,  i-.ti  |  1 1 uvn );  7(1,  :id;iiKi  iliisin. 

;!.  P.  I.  Xtarslon  and  I).  I.  Kmgshury,  "Scattering  b\  a  hobble  in 
water  near  the  cril  ical  angle;  mlcrlcrcncc  el  I  eels,"  J.  I  ipl .  Sin 
Am.  71,  lir.'  1  in*  1 1 !is 1 1 , 

:l.  I  i  Mie.  "Itei  I  rage  /ur  Opt  ik  trailer  Medina,  spozicll  kulluidnlcr 
MotalloMingon."  Ann.  I  Mi  vs.  ( l.eipzigi  25. 1177  t  la  i  IlliiKi 
I  1 1.  ( '  van  do  I  InUl .  I. u;h  I  Srttl  I  I’l'int!  In  Small  I 'hi  In  /is  I  Wilev, 
New  York.  lll.'iVl 

a.  W  .1.  W  i-ecanbe.  "Improved  Mie  scattering  algorithms,"  Appl 
1  >|ll  111,  latla  1:1(111  I  1 1IHIII 

ii  II  1 1  Roll.  R  (1(  iuinprechl .  and  C.  M.  Sliepi  eviih,  '"I  heorel  a  al 
light  si  altering eoelln  ienls  lor  relative  retractive  indexes  less  than 
unitv  and  lor  totally  reflecting  spheres."  .1  (Ipl  Soe.  Am.  It, 
IH  |  I  Ilia  1 1. 

7.  R.  11.  Roll  el  ill  ,  Tahir s  o /  l.inlil  Snil Irnim  Initul  inn*  IF 
Mulligan  Press,  Ann  Arbor.  Mich.,  I'.lhXi 
s.  1,1,  /.ehnanovieh  and  K .  S.  Shifrin.  Tahir*  o/  l.iuhl  Sm  Hr  ruin. 
I'url  III  Coi7/ieieals  o/  h,\ limliaii  .S',  alien  in; ,  mill  l.whl 
/Vissnri'  1 1 1  vilromeleorologieal  Publishing  House.  I.einngrad, 

IllliSI. 

P  M  Kerker.  77e  St  nttri'iiu:  n/  I, mill  am/  Ul/u  i  /■.’/ , a  Irniiiuiiiirln 
liinlnii  inn  I  Academic.  New  York ,  1 1  till  1 1 . 

In  I  V.  1  lave.  "Seal lering  ol  visible  light  by  large  water  spheres," 
Appl  (Ipl  S,  I. ‘.a  llil  HllliXI. 

II  W  .1  Wiscoedie.  "Mu  Sealleliug  t ‘alculal ions;  Advances  in 
'  1  '< -<  Il  1 1  i<  |  to  and  Fast.  Vector  Speed  <  oiupulrr  <  'odes,"  NI'AR 
I'ecl  a  mat  Note  Nl’AII  TN  1  III  I  STR  I  Nat  ional  Ceiiler  lor  A I 
mosplierie  Researeh,  Rouliler.  Colo,.  III71II. 

1"  (  ,  K  I  lavis,  "Sealiering  ol  light  hv  an  air  liulilile  ill  water."  ■).  I  >|i! 
Soe,  Aim  la.  aV.t  ...HI  I  I li:.  .  1 

1:1  R,  1 1.  Johnson  ami  l\,  C  I  'nuke,  "  Rubble  pupil  lal  inns  and  speet  rn 
ill  roaslal  waters:  a  pliol ograp'.ie  appl'uaeh,"  J,  1  ieupliys.  lies, 

Nl,  :i7lil  aVlili  I  1117(11. 

II,  A  Keller.  “Kill  Slreulieht  /.alilverlalireh,  aiigcwandt  rur 
Reslinunung  ties  Kavilalionskeiinspeklrunis.''  Ilplik  .'17,  Ilia  1 Y < i 
( I'lYlll. 

I  a,  A .  Keller,  "The  inl'luenee  oil  lie  ca  vital  ion  nucleus  spec!  rum  nil 
cavil  al  ion  iucepl  ion,  invest  igaled  with  a  seal  lered  hglil  ci  nulling 
md hod. "  J,  Rasic  Kng  III,  PC  ll'.’a  I  HIV:.’ I. 


SUPPLEMENT  TO  PAPER  NO.  4 


This  supplement  compares  unpublished  Mie  computations  (solid  curve) 
with  the  physical -optics  model  from  Paper  3  (dashed  curve).  The  refrac¬ 
tive  index  ratio  is  n./n  =  0.75.  The  notation  is  as  in  Papers  3  and  4 

i  o 

except  in  Figures  SI  and  S2  we  give  the  total  normalized  scattering 
I  =  (1^  +  1^)  / 2  for  the  case  of  unpolarized  incident  light.  For  Figures 
S3  -  S5,  scattering  for  the  same  ka  but  for  the  orthogonal  polarization 
is  given  in  Paper  No.  4;  the  associated  figure  from  that  paper  is  noted 
in  parenthesis. 

FIGURE  CAPTIONS  FOR  SUPPLEMENT 
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Log 
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I  for  ka 
u 

=  25. 

Fig. 

S2 

Log 
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Iu  for  ka 

=  100. 

Fig. 

S3 
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for  ka 

=  100  ( 

see  F ig. 

3b). 

Fig. 

S4 

>2 

for  ka 

-  1000 

(see  F ig. 

3a ) . 

Fig. 

S5 

5 

for  ka 

=  10000  (see  Fig 

,2a).  See  a  1  so 

Fig. 

4  of  Paper 

1  for  (J) 

from  75“ 

’  to  85°. 

Fig. 

S6 

Pha 

se  difference 

for  scattering  amplitudes 

(as 

in  Fig.  5a, 5b) 

for  ka  =  25.  Note  that  the  model  does  not  give  the  fine 
structure  present  in  the  Mie  result. 
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One  of  the  causes  of  light  scattering  in  glass  is  the  pres¬ 
ence  of  small  gas  hubbies.  In  critical  applications  such  as 
optical  fibers  this  scattering  is  not  negligible;  a  microbubble 
trapped  at  the  core-cladding  interface  may  scatter  a  signifi¬ 
cant  fraction  of  the  light  in  the  fiber.1  Here  we  use  Mic  theory 
and  physical  optics  to  describe  the  main  features  of  bubble 
scattering  in  glass.  We  anticipate  that  these  results  will  be 
useful  for  the  detection  of  bubbles  in  quality  control  appli¬ 
cations. 

We  are  unaware  of  prior  publication  of  the  Mie  intensity 
for  0  180°  scattering  from  bubbles  in  any  dielectric.  Struc¬ 
ture  in  the  Mie  scattering  from  bubbles  can  he  related  to  dif¬ 
fraction  and  interference  near  the  critical  scattering  angle  and 
to  the  Brewster  null  in  the  reflectivity  of  parallel  polarized 
light.  For  scattering  from  a  sphere  with  a  refractive  index 
which  exceeds  that  of  the  surroundings,  the  critical  scattering 
angle  is  not  present;  also,  the  effect  of  the  Brewster  angle  is 
less  pronounced. 

We  approximate  the  bubble  shape  as  a  perfect  sphere  and 


assume  the  glass  has  a  refractive  index  of  1.46,  as  does  the  glass 
used  in  the  cladding  and  outermost  core  layer  of  a  typical 
optical  fiber  preform.  (In  the  resultant  fiber  a  bubble  trapped 
at  the  core  -cladding  interface  is  broken  up  into  a  sequence  of 
much  smaller  bubbles  elongated  in  the  direction  the  fiber  is 
drawn. 1 )  We  take  the  bubble  gas  to  be  air,  so  that  the  relative 
refractive  index  n,/n„  cv  0.685.  A  bubble  is  described  by  its 
size  parameter  ka,  where  k  =*  2ir/A,  A  is  the  wavelength  of  the 
light  in  the  glass,  and  a  is  the  bubble  radius.  Here  we  con¬ 
sider  bubbles  with  size  parameters  5,  25,  and  100.  For  A  = 
0.43  pm  (which  corresponds  to  632.8  nm  in  air  for  visual  in¬ 
spection)  ka  =  100  gives  a  =  6,9  pm;  ka  =  5  gives  a  =  0.35 
pm.  For  A  =  0.56  pm  (0.82  pm  in  air,  a  typical  operating 
condition  of  fibers)  ka  =  100  gives  a  =  8.9  pm. 

If  the  bubble  is  illuminated  by  a  collimated  light  beam,  the 
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Kig.  1.  I.ognrithm  (bsse  10)  of  the  normalized  scattered  intensity 
predicted  hy  Mie  theory  (solid  curve)  for  ka  -  100  and  the  electric 
field  perpendicular  to  the  scattering  plane.  The  dashed  curve  is  the 
physical  optics  approximation  given  hy  Kq.  (24)  of  Ref.  It  which  is 
useful  when  scattering  angle./)  is  close  to  the  critical  scattering  angle 
./.,.  ~  94°.  The  approximation  fails  to  describe  forward  region  (./•  < 
20°)  and  the  backward  regt  >n  0/i  ~  180°). 


Fig.  2.  hike  Fig.  1  but  with  electric  field  parallel  to  the  scattering 
pintle  (y  =  2  scattering). 
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Fig.  It.  Scattering  for  ka  =  25.  Mie  theory  for  j  =  I  (dotted  curve 
and  7=2  (solid  curvet.  Dashed  curve  is  /  j  from  physical  optic 
approximation. 
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lingular  distribution  of  the  scattered  intensity  is  described  by 
Mie  theory. -  At  a  distance  R  »  ha-  from  the  bubble  center 
and  a  scattering  angle  </>,  the  ;  -polarized  scattered  intensity 
is  given  by  the  incident  ./-polarized  intensity  multiplied  by 
/;  (</>)( where  j  =  1  denotes  polarization  of  the  electric 
vector  perpendicular  to  the  scattering  plane,;'  =  2  the  parallel 
case.  The  normalizat  ion  of  the  funct  ions  Ij('t>)  is  chosen  so 
that  /,((/>)  =  1  represents  perfect  reflection  according  to  geo¬ 
metric  optics.2  '1  Using  a  program  described  in  Ref.  5  incor¬ 
porating  a  modified  version  of  the  MIKVO  subroutine  devel¬ 
oped  by  Wiscombe, 11  we  obtain  the  />((/>)  functions  shown  in 
Figs.  1  4. 

The  coarse  structure  in  the  vicinity  of  the  critical  scattering 
angle  </>,  may  be  understood  with  reference  to  a  physical  optics 
approximation  developed  by  Marston.2'7  Most  of  the  scat¬ 
tered  intensity  in  this  region  is  due  to  three  types  of  rays  which 
we  denote  by  the  parameters  p  =  0,  1 , 2,  where  p  is  the  number 
of  chords  within  the  bubble  and  the  p  =  0  ray  has  only  an  ex¬ 
ternal  reflection.4  When  the  local  angle  of  incidence  0  satisfies 
sinfl  >  sinW,  =  n,/n„  the  p  »  0  ray  is  nearly2  totally  reflected; 
the  abrupt  boundary  of  total  reflection  at  0,.  results  in  a  far- 
field  diffraction  pattern  of  this  ray.  Rays  with  II  <  II,.  are 
mostly  transmitted  inside  the  bubble.  Rays  with  p  —  1  are 
refracted  twice  and  interfere  with  p  =  0  rays  in  the  region  </» 
<  (/>,.,  where  </>, .  =  ?r  —  2 II,.  is  the  scattering  angle  of  rays  inci¬ 
dent  at  II,..  Mnrston’s  physical  optics  approximation  con¬ 
siders  the  diffraction  and  interference  of  p  =  0  and  1  rays. 
Small  quasi-periodic  fine-structure  variations  ride  on  the  large 
quusi -periodic  coarse-structure  variations  and  are  due  mainly 
to  the  interference  of  certain  neglected  p  =  2  rays  Ip  =  2'  rays 
in  the  notation  of  Refs.  It  and  5).  If  the  average  intensity  over 
an  angular  range  containing  several  fine-structure  variations 
is  desired,  the  approximation  may  be  used  to  predict  /,  in  the 
range  20°  S  </>  S  »/»,.  =  911.54°.  This  model  breaks  down  for 
small  bubbles,  as  can  be  seen  in  Fig.  4,  due  to  its  use  of  the 
Fresnel  reflection  and  transmission  coefficients  of  planar  in¬ 
terfaces  and  the  stationary  phase  approximation  of  diffraction 
integrals. 

The  physical  optics  approximation  predicts  a  quasi-period 
5  (A/«),/2  rad  for  the  coarse  structure,  and  a  fine-structure 
quasi-period  of  0,95(A/«)  at  </*,• ,  decreasing2  with  </>.  These 
features  are  clearly  evident  for  ha  =  25  and  100,  as  is  the 
predicted  broad  decline  in  intensity  for  -/>  >  </>,.  - 
arcsin|1.2(A/u  eosfl,)l/2|  (identified  in  Ref.  5  ns  the  approxi¬ 
mate  location  of  the  last  coarse-structure  maximum). 

The  broad  minimum  in  l  >  near  1 10°  corresponds  to  light 
scattered  from  the  Brewster  angle:  parallel  polarized  light 
incident  on  a  plane  glass  air  interface  is  totally  transmitted 
when  II  =  114.41°.  Thus,  when  ha  is  large  enough  to  use  the 
Fresnel  coefficients,  we  expect  a  negligible  contribution  to  /■.> 
at  i/i  =  180°  -  2014.4°)  =  1 1 1.2°  from  the  externally  reflected' 
(/)  =  0)  ray.  This  >l>  can  also  be  written  as  2  arctan(u„/n, ). 

For  ha  =  5  (Fig.  4),  the  line-structure  quasi-period  is  ap¬ 
proximately  equal  to  the  coarse-structure  quasi-period.  The 
minimum  in  /•..  at  102°  and  the  slight  dip  around  45°  are  t  he 
last  remnants  of  this  structure.  As  ha  is  reduced  below  5,  we 
find  that  the  scattering  pattern  approached  that  of  a  dipole 
radiator,  as  predicted  by  Rayleigh  scattering  theory:  / ](</>) 
becomes  a  flat  line,  and  /•_>( </>)  is  symmetrical  about  a  minimum 
at  1)0°. 

We  have  compared  these  plots  with  ot  hers  made  for  n,/n,, 
=  0.75  (air  bubbles  in  wafer)  for  the  same  values  of  ha  and 
found  them,  as  expected,  quite  similar.  The  most  noticeable 
difference  is  the  slightly  larger  quasi -periods  for  bubbles  in 
glass.  In  both  cases,  the  best  angles  to  look  for  scattering  will 
be  those  angles  where2  </>  5-  </>,•  —  (A/<.<)1/2.  Corrections  are 
sometimes  needed7  when  computing  the  effect  ive  scattering 
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Fig.  4.  Like  Fig.  3  but  with  ha  =  5. 


angle  observed  by  a  detector  due  to  refraction  of  the  scattered 
light  as  it  leaves  the  glass  and  enters  the  detector.9  The  coarse 
and  fine  structures  present  when  ha  Z  25  may  be  used  to 
obtain  independent  estimates  of  a.  That  diffraction  is  im¬ 
portant  at  i /),.  can  be  seen  bv  comparing  the  Mie  results  to  the 
geometric  prediction,4  including  only  thep  =  0  ray,  that  /,(</>,.) 
=  I  and  with  the  physical  optics  prediction  that  /;(</>,.)  =  0.25; 
the  Mie  results  lie  much  closer  to  the  latter  than  the 
former. 

The  result  that  scattering  is  best  observed  when  i/>  5  0,.  - 
(A/a)1/2  rad  can  also  be  used  in  certain  cases  for  scattering  by 
nonspherical  bubbles  if  their  approximate  size  is  known.  For 
example,  if  the  bubble  is  spheroidal  and  either  the  major  or 
the  minor  axis  lies  in  the  scattering  plane,  the  critical  scat¬ 
tering  angle  predicted  by  ray  optics  is  the  same  as  that  pre¬ 
dicted  for  a  sphere.  The  diffraction  reluted  shift  in  the  first 
maximum  is  roughly  ( A/u),/2  rad,  where  a  becomes  an  average 
(which  depends  on  the  orientation  of  the  spheroid)  of  the 
semimajor  and  semiminor  axes.  The  details  of  the  coarse  and 
fine  structures  will  depend  on  the  bubble's  shape. 

The  approach  to  the  scattering  described  here  is  to  allow 
the  light  to  enter  and  leave  the  glass  via  plane  surfaces  and  to 
neglect  surface  reflections.51  Scattering  near  90°  has  been 
used  to  detect  bubbles  in  glass.111  For  glasses  wit  h  «„//»,  =  1.5, 
</>,  =  90.4°  and  the  Brewster  scattering  angle  is  112.6°.  Pre¬ 
vious  applications  of  Mie  theory  to  bubbles  are  reviewed  in 
Refs,  5  and  I  I  Approximations  for  scattering  by  inhomo¬ 
geneities  within  fiber  waveguides  have  recently  been  pub¬ 
lished. 1212  but  these  do  not  describe  the  angular  structure  lor 
scatterers  with  ha  >  I . 

This  work  was  supported  by  t  he  Office  ol  Naval  Research. 
1’.  1,.  Marston  is  an  Allred  1*.  Sloan  Research  Fellow. 
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Append i x 
COMPUTER  CODES 


Three  programs  were  used  to  generate  our  plots.  MIESC, 
incorporating  a  modified  version  of  Wiscotnbe's  MIEVO,  produced 
the  Mie  and  CRIT  the  approximation  results.  The  versions  given 
here  store  results  in  disk  data  files  where  they  can  be  accessed 
by  the  plotting  routin  MIEPL  for  display  on  the  HP  2847A  monitor 
or  for  paper  plots  on  the  HP  9872A  plotter.  All  programs  are 
written  in  Hewlett-Packard  RT-IVB  Fortran.  CHIT  evaluates  t  lie 
approximation  developed  by  Mnrston  ami  Kingsbury  (1981). 

MIESC 

MIESC  asks  for  the  following  inputs: 

1.  REF  INDX  is  the  complex  relative  refractive  index  n. /n0 . 
The  imaginary  part  must  be  non-positive. 

2.  BETA  is  the  size  parameter  ka  ~  2rra/A,  where  A.  is  the 
wavelength  in  the  outer  medium. 

3.  AST  is  the  smallest  scattering  angle  for  which  scattering 
angle  for  which  scattering  functions  are  computed,  DPHI 
the  angle  step  size,  and  NUMAMG  ~  Na„g  is  the  total  num¬ 
ber  of  angles. 

If  NUMANG  is  even,  results  are  calculated  for  the  angles 

<Pj  ~  AST  +  i(DPHI),  i  =  0  to  NUMANG /2 ,  and  their  complementary 
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angles  (180°-  <pf  )  .  (If  NUMANG  =  1,  90°  must  be  the  only  angle. 

If  NUMANG  =  0,  only  the  extinction  efficiency  factor  Qtc#  is 
computed . ) 

For  example,  if  n,  =  1.33  -  i(.08),  n  =  1.0,  and  scat¬ 
tering  functions  are  desired  for  <p  -  40*  to  90°  with  increments 
of  .5°,  one  should  type  in 

/. 7508, -.0282, 40.  ,  .5,101 

MIESC  then  asks  for  the  name  and  cartridge  number  of  an  existing 
datafile  in  which  to  store  the  output.  What  data  is  stored 
depends  on  what  the  output  buffer  RBUF  lias  been  equivalenced  to. 
In  the  version  included  here,  each  record  v/ill  consist  of  three 
real  numbers  giving  <P ,  I,  (<|>),  ?nd  Ij(<£),  and  the  records  are 
entered  in  order  of  increasing  4> ,  except  that  the  record  for 
(180*-  <j> )  immediately  fe  .lows  that  for'  <j>  . 

If  the  phase  difference  DIFA  is  not  desired,  the  code 
that  computes  FH1,  PH2 ,  and  DIFA  is  skipped. 

To  bring  subroutine  M1EV0  down  to  more  manageable  size, 
the  original  comments  were  largely  omitted.  The  commented  ver¬ 
sion  is  given  in  Wiscombe  (1979).  Most  complex  variables  have 
been  converted  to  2 -dimensional  arrays  or  to  separate  variables 
prefixed  with  "R"  or  "I."  For  example,  NUM  became  NUM(l)  ana 
NUM(2),  denoting  the  real  and  imaginary  parts,  respectively,  and 
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ZET  became  RZET  and  IZET.  (This  was  done  because  of  the  absence 
of  double  precision  Fortran  in  RT-1VB  Fortran.)  Most  complex 
arrays  were  similary  converted  to  separate  arrays  for  the  real 
and  imaginary  components. 

Subroutine  names  prefixed  with  "DV"  refer  to  functions 
in  the  HP  Vector  Instruction  Set,  which  considerably  speeded  up 
the  execution  of  two  former  DO- loops. 

To  keep  the  loaded  program  within  a  32K  size  limit,  the 
original  MIEVO's  complex  BIGA  array  was  discarded.  (BIGA,  which 
stored  the  A„ ,  contained  as  many  elements  as  the  number  of 
terms  used  in  the  Mie  series.  Thus  for  ka  =  10  000,  a  total  of 
2(10  088)  =  20  176  double  precision  must  be  stored,  requiring 
161  408  bytes.)  Required  An  values  are  instead  stored  on  and 
retrieved  from  a  system  disk  using  EXEC  calls. 


The  use  of  the  Hewlett-Packard  minicomputer  instead  of  a 
regular  "pay"  computer  system  was  absolutely  necessary.  Extra¬ 
polation  from  Table  8  of  Wiscombe  (1980)  indicates  a  MIEV0  exe¬ 
cution  time  of  over  100  seconds  on  a  CRAY-1  for  ka  =  10  000  and 
NUMANG  =  8000.  The  several  size  parameters,  repeated  testing, 
and  dozens  of  plots  required  made  our  choice  of  computer  systems 
inevitable. 
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CRIT 

Required  CRIT  inputs  are  the  relative  refractive  index 
RI  (now  assumed  real);  the  size  parameter  BETA;  PTS  and  STRT, 
which  determine  the  angles  for  scattering  computations;  and  LAST 
ANGLE,  the  largest  PH  (c p)  .  N  =  u0/ni  -  1/RI  is  used  to  calcu¬ 
late  the  critical  incident  angle  0C  : 

THC  =  arcsi n ( 1/N ) . 

(Arctangents  are  used  in  the  program  due  to  the  unavailability 
of  the  arcsin  function.)  The  critical  scattering  angle  is 

then 

PHC  =  Tr  -  2  ( THC  )  . 

For  a  given  incident  angle  TH-THC,  the  angle  PH  at  which 
the  refracted  ray  is  scattered  is  given  by  an  equation  derived 
using  Snell's  Law  and  elementary  geometry 

IH  -  2 ( arc  sin S  sin(TH )/N )  -  TH). 

For  a  desired  smallest  PH,  one  chooses  STRT  such  that  TH 
=  STRT(THC)  substituted  in  the  above  equation  gives  this  PH. 

As  CRIT  executes,  TH  i s  stepped  up  in  increments  of 
E  =  (THC  -  THillit  ) /PTS .  As  TH  approaches  THC  the  corresponding 


step  sizes  for  PH  grow  steadily.  Consequently,  at  selected 
angles  TH2 ,  TH3 ,  and  TH4 ,  the  step  size  in  TH  is  reduced. 


The  function  G1  is  the  geometrical  divergence  factor 


sm  ( TH)  cos(TH) 

G  I  rr  — i — - - - 

Zf  1  -LrtcM(TH)/cos(R)J|si»,<rPH) 

in  the  expression  for  the  intensity  of  refracted  radiation 
(Davis,  19S5;  van  de  Hulst,  1957),  where  R  =  arcsin(Nsin(TH) ) . 

IP  and  IS  are  the  intensities  of  the  parallel  and  per¬ 
pendicular  refracted  radiation,  respectively,  computed  using 
Fresnel's  intensity  transmission  coefficients  (1  -  FP )  and 
(1  -  FS)  for  plane  surfaces.  The  coefficients  are  squared 
because  this  radiation  is  refracted  twice  --  at  entrance  and  at 
exit . 


The  variable  W  is  the  upper  limit  of  integration  on 
the  Fresnel  integral  used  in  the  reflected  radiation  calculation 
(Marston,  1979).  Subroutine  FRES  evaluates  this  integral  using 
a  numerical  approximation  (Abramowitz  and  Stegun,  1965). 

The  phase  differences  B1  and  B2  between  the  reflected 
and  refracted  rays  for  perpendicular  (Bl)  and  paral  lel  ( B2 ) 
polarization  have  three  component:'.:  the  phase  shift  the 
reflected  ray  undergoes  at  the  bubble  surface  (DELI  and  DEL2 ) , 
the  phase  path  difference  E’J’AC,  and  a  Tt/4  phase  difference  aris¬ 
ing  from  the  curvature  of  the  scattered  wavefront.  (For  the 
reflected  ray  this  phase  shift  is  already  included  in  the  Fres¬ 
nel  integral . ) 
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CRIT  is  now  ready  to  compute  the  scattering  functions  at 
angle  PH: 


s,  -  tflS  +  (l/,{2){  FC  +  iFS ) exp ( - iBl ) 
S»  =  VI?  +  (i/f2)(FC  +  iFS )  exp  ( - iB2 ) 

where  FC  and  FS  are  the  output  of  subroutine  FRES. 


MIEPL 

Output  was  plotted  by  program  n'.TEPL  using  the  HP  1000 
Graphics  package.  MIEPL  is  not  included  in  the  program  list¬ 
ing  which  follows. 
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^  ^  >  ^  .  r  Ht.MiT  L  tHi"  '‘liUOV  1  M!|  f:  0  f  ¥  1  *  >  >  f  >>  f  f  V  •<■  *  *  i*  '■¥ '+  f  •*  ♦  ;' 

I  t  f  f  f  t  f  ♦  ^  f  f  ►  ♦  f  t-[il.l  I  UM  T  r-  1  i  jij  ,'t>UK  i  t  i  ••  L.  *-  * 

BLOC!'  C'hTh 

couple  ::  s&mok  ,  s t ,  S2 

0  0  »j  U  L  i.  p  p.  e  c  1  S  1 0  h  e  e  r  h  ,  x  n  u ,  R  1 

COMMON.  ■  l  NOU  r.'  PI  ■  2  > .  '.■•Bwft-  ,  SK  *51  2  > .  S2-:  S  I  2  '  .  OLTw  ;.MIK  512  • .  H 20 J  ' 

♦  I  lUMHNG ,  i  IMP  I  .MM.  OO-;  f  ,  L'JOUT  ,  tCNWD,  U-LM  ,  I  '^El  T  .MP0,  I  UP  1 U,  1 Jp  tU 

END 

ppogphm  m t  f sc  • 

c  HHGULwR  C'  I  MEM  '.  I  ON'S 

couple;'  sbnci  ,  si.  02 

DOUBLE  PI*  El  ] 1  DM  BL  Co,  HMU,  F'l ,  t*  •* .  R  I  ,  UJM  .  HuT 


PEmL  M2CUT  .  I  1  ,  I.',  IN  '.'CM 

DIMENSION  IPCe.  !  M  '  ppup.'  ?  >  NmMEC  ?  >  _ 

COMMON/  IMOMT/M--  2  »  '  K  •  1  •  *f<12  ’ .  S2'-  PM  2  < ,  BE  TO  ,  /MU'.  5  ’  2  >,  M^L-  . 

+  NUMmHG  ,  MNP  1  ■  MM  .  OEHT  ,  I.UOUT  ICNWD, ILEH  IGCCT  MB  u  ,  1  UP  U',I  UP*  O 
EQUIVALENCE  '  POUF  '•  1  !  ANGLE  PBUFv 2  .It  •• ,  •  PLOP  3  >  •  1 2  .» 

PNTm  i:  ir.'l  1.1.  • .  I  TYPE/-?/,  i  rc/li.1,  I  OF'TN/0. IL.  0/ 
c  IL  =  t>IM'  PGUP  .'*2 

LU  I  U“LuGLU'  .ID.) 

WR  II E<  LU  t  M .  2  0  0  0  »  .  ..  .  _  .. 

2000  FORMaTt  "PP  I'l  rnt-PI  •  rOPCOUT  ,  COMPUTE  PHhgEu  ’  '■  1  r  *  jP  »■ 

READ  '•  LU  I M  •  2 '.'0  1  >  IP,  t  ■=.  r  ,  t  PH 

2001  FORMAT  t 3h1  ' 

IF  <  I  R  .  EO  ,  1  MY  ■>  LUOUT  - 
IF  (IP.ME.lHYi  LUC'i  IT  -  l.UIM 
M2CUT  =  0. 

PI  a  3  ,  1  -1  1  G  I'YC/'  YSODo 
COM  -  R  I  /  .  1  0D3 
■SCOt I  -  SMGl.'  COM  * 

•  1, 02  ^  1 HDY ,  BETA.  DBF  T  mST.  ''PHI,  NUMHNG"  > 

*“  '  PFul')'  I  u  { 1 1 ,  •  r-t'  t  ).RD  BETH. DDET,mST, DPMI  .m.iriMHG 

c  I  i  CMS  FOR  NUMERIC  Ml.  INTEGRATION  0'»EF'  SOLID  ANGLE 

ShST  -  S  I  H'  C.MGL'  AS  1  1  «='  C  ON  ’> 

OoTEF:  =  ■  F'f" H  1  t-'.-CC'M  ' <■  *2 
1FC  I  ‘,T  ,  1  If:  .  1  ■•'<  <  GOTO  2 
WRITE'.  L  M  I  M  .  2003  ■ 

2003  FQF'Mm  1  ■.  “FtlT  1 1  .*  nif.".  CHS-.  CR  MO,';:  CM,'"' 

PE.!<  l,  i  j  I ! : .  .  04  )  ‘l-.oir.  j  •  .  t  ihI'IE’  2  >  MhITID  J’-t'-R 
;  nil  i  F'.il-  Mu  f-  ?u2  i  .  ; 

I  AU.  OPEN'  K"'B  ,  IEP  MhME  ,  1 0P  TM  .  I  SC  ,  I  l  ■' 

It  ■  f  !-'Kp  ,  l_  V  0  •  GO  ro  B,v? 

Pp.MF'  I  ■  -  BE  Tu 

rbi. ip ■  ::■ .»  =  r,  i'  i  .' 
ppijr  s  ■  -  r  i  _  ■ 

Cull  wi-  i  n  ■  ir-  i..  t. . Rour ,  sl  > 

mi i pi 1 1 1 j B!  -  •  ir  i-:,ii n  ., '5  I I  .......  , 

i-  i  r  Mi.ir.MuL  ,  >1  i  T  I.t  1  DESIRED  ,  M'IlT  PM  _•  r-’  i  ■•’ ■  *  -  B  •  ••  *  *  *■ 


Hi  OTu  =  I'li  l"  .  i'.h-Im'  I',  .  j  '•  2 
C-O  14  IT- ■  =  1  .  GUM'!:  O't 


NMR1  -  NO  MhNiS  +  I 
NH  =  HMP  1 1-  2 
DO  4  I  *1,  nil 

<4  XMLK  I  >  =  OCO  h,5T  +  i  1-1  '  .  D  P  H I  >+CL'N 


*  CALL  MIEVO-  n;  , 

IF  0  NUMANG  ,  EG  .  0  .•  GOTO  12 
TINT  «  U. 

DO  10  I  =  1  ,  (HI 
RK  =  0  . 

UNCLE  =  m  it  T  + 1  1-1  ’+PPHI 

PHI,  PH2  A  PC  THE  PH, .OF.  » ‘MOLES  IH  pCCPECS  ;  PHD  IF  IS  DIFFERENCE 
IF  f  I  PH  .  HE  i  HV  >  COT 1. 1  ;• 

PH  1  - -57  .  OOF:;;  f,  h  i  hM1  hIHaG-  ■=■  1  •  I  ■  '/REAL'  1 I  '  >  ’ 

PH0--57  LOSS  »•«  qToili  wtHAG*  i-~>'  I  1  VR£+.L>  S2*.  I  "< 

IF\ REAL'  S  1  ■  I  ■  LT  .  u  .  >  PH  '  *PH  H  i  PO  . 

IF.'  PH  1  ,  QT  ,  \  \{  0  ,  i  PH  1  -PHI  -  2>2  u  . 

IFC  REwL<  S2'  T  ■  > ,  i_  T  .  0  ,  >  PHE=PHO  M  00  . 

IF>'  PM2 ,GT  .  100.  >  PH.l  -PH2~~6o. 

PHD I F®PH2-PH l 

IF<  PHD  IF  ,  OT  ,  15«C».  «  PHDIF=PHD]F-:n'.D, 

IF< PHDIF . LT , -  1  SO .  *  PHD  IF -PHOIF  +  3G0 . 

?  li  «  <  REAL’’ Si'.  I  »  ■  > +2  +  •  A  IHaC>.  1 '  I  »  )  2 

11  a  .'li  +0:1; :« 

12  =  <  REaL<  -:.2«  I  >  '  )**£■*<  hIMaGv  ‘*.2<  1  >  5  1  Z 

12  =>  <12  +  02  .' 

I H T E 1 1  a  ALOOT'.  ,  Oh.;  I1+I2  )  > 

DEGPOL  *  ■  12-11  12  +  11  '< 

IF<  AMO  DC  hUi'iLE,  2.  .»,  EO  ,  0  .  >UR'i  TE'.  LL'OUT  ,  1  0  0  1  >  ANGLE , I  1  ,  I  2 
tFC  l:j f  .£«•!.  !HV'»  CALL  WKITFC  IDCB,  lERR.R&UF  H  ' 

DO  COMPLEMENTARY  ANCLE 

IF  CO  HP  .  ANGLE  n.liiC- *  I  Otic  DEC  IP.  ED  ,  i  HE  FINAL  DO-LOOP  IN 
SUBPOUT  I  He.  MIEVO  CmHNOT  EXCLUDE  S'  MNP  1  -  J  *  FH'S 
ANGLE  a  ISO  -ANGLE 
IF<  HUHhNG  .  El.' ,  1  I  GO  IC  1  2 

I  1  C  -  '•  PF.uL  ■  G  1  v  HMP  1-1  ■  i  '»  +2  +'.  h  J  HAG".  C  1  *  .IMF  1-1  •  •  >**2 

II  =  -11 +02  • 

I2C  =  '.PEhL  •••2>  MtTF  1  -  1  .>*+2+'  h  I  WAG'.  S2f.  Nl'lf  1-1  ■  >  '+*2 

12  *  <  12+02) 

INTEN  =  ALOGT.  ,  0 *  '•  11+12  *  • 

THESE  LINES  PERFORM  SOLID  ANGLE  IOTEGP+  ION  wT  PM I «P 1/4 
3  2INT  a  <  COS'-  RL  + -.CON  *  »'•+•♦  2  +••  11  +  11'-  *  +  '  OH*  1  \  CN  *  m  +2+>:  I  L  +  I  21-  > 
TINT  a  TINT+l  »EIMT  +  f"-;.Tf.r' 

RK  a  RF.+DPHI 

IF  >  S  I N'.  hNGLL +SC '.'H  '+C0'";.'.  RT'  +SC OH  '.  C.E  r  '  GO  r0  8 

IF-'  HNC.r"’  .‘HOLE  /  2  '.CO  h.  1  WRITE  ’  LUOL<r  I1'1’!  wHC,:.  E  ,  i  I  ,  T  2 
til  CONTINUE 

RUNT  mLOGT'  T  I  hj  , 

ANGLE  =  EETh 

11  =  FLINT 

12  a  PL INT 

IF*;  l  1  H  ,'  >  CmI.L  I.If  1TF<  ICCC  .  ILFf  .  'F  i  L. 
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IE 

r  a 
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2  007 
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1  000 
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C 


IF<  1ST  .  GO  .  I  HV  >  ChLL  CLOSE'.  I  DCS  ' 

WRITE'-  LU  IN,  2  007  »  I  ERR 
FORMA  I'  "  I E  R  R  =  "  ,  15  > 

USE  THIS  FORMaT  FOR  MORE  DETAILED  OUTRi.IT 
FORMAT'  1H5..5X,  HIE  SI2E  PaPaMGTEP  =  FS  .  2 , i  5X ,  REFRACTIVE 
k  'IHoE:;  =  '  .  r>1  4  ,  6  ,  Et2 .3//  ANGLE  S-SUb-1  ',2ix,  'S-SUB-2', 

♦  15:-',  I NTENS  I  TV  ,  2X  ,  OEGPOL  .8.  ,  '  I  1  ‘  .  1  2/  12'.' 

FORMhT <  FT  ,  2  ,  5  E  i  4  ,  6  ,  FS  .  4 , 2E  I  4 . 6  ,  E6  .  3 
FORMAT"  F7 , 2, 1-  1  2  ,  ~  > 

END 


SUBROUTINE  NIEVCk  IU  > 

COMPLEX  SBACK,  St,  S2 
DOUBLE  PRECISION  BET H , NMU , P I 
REAL  N2CUT 

COMMON/ 1  NGUT/P I  >:  2  • .  SBACK  .  SK  Si  2  >,  S2-;  51  2  >,  BETA  ,  XMU"  512  >,N2CUT  , 

*  HUM  AUG ,  NMP  1  ,  Hi  i .  GE/T  ,  lUOUT  ,  KNUD,  I  LEM,  I  SBC  T,  MR  0,  IOP1  0,  IGF  20 
LOGICAL  NO  I  MAG ,  NOhNGS 

COMPLEX  ANI'li.BNMI 

D I  MENS  I  ON  RSP<  256  > ,  I  SP'  256  > ,  PS  M<  256  > ,  I  SN-  256  ■  .  PSPSC  256  ,  I  SPSC  25' 

*  PSMS'i  256  > ,  I  SMS'  256  '.PIN'.  256  .'.PI  NN 1  <  256  • .  1  HP'-  256  >  ,  TAUN<  256  > 
DOUBLE  PRECISION  RBICA.  IBI  C.A  ,  P  IN  ,  P  INN  1  ,  NP  1DN ,  COEFF  ,  VTEM'1  256  !> , 

PSP  ,  I  SP  ,  PSM  ,  T  Si-1 ,  PSPs  ,  I SPS ,  PSMS  ,  I SM  .?• ,  AN'1'  2  .  BN*.'  2  >  ,  RANP  ,  I ANP  , 

*  PBNP  ,  I  BNP  ,  RhNPM  ,  IhNPM  .  RBI  1PM  ,  IBM  Art .  PZET  ,  I  CRT,  P2ETN,  12ETN, 

*  RZTNP1  ,  IZTHP  I  ,  TMP  .  TaUN  .  RE 2,  REZINV  ,  FN, RN, FNP  1  ..  TWONP  1  ,  RNP  1  , 

+  2 1  NV<  2  >  ,  FFC  2  ' ,  AC'.  2  >  ,  DEN"  2  .» ,  MUM';  2  >  .EPS!  ,  EPS2  ,  TEM<  2  > ,  TMAC  2  >  , 

*  NTN<  2  >,  DTD'.:  2  .- .  TK  2  >  .  PI  OR  I  V  .  PS  IN,  PS  IMP  t  ,  CHIN.CHINP1  ,  BUFFI  32  > 
EQUIVALENCE  <  S  1  •  1  '..PIN',  1  >  >  ,  •.  S2n  1  >  .  TAUN"  1  -  ». 

*  S 1  >'  257  .> ,  P I NM 1  <  1  ;»  > ,  <  S2<  25  7  > ,  TMP'  1  >  ’ 

DATA  EPS  1  /  1  .  0-2/  ,  EPS2/1  .  0-8/  .  MAMIT/  10000/ 


IFC  BETA  .  LT  ,  0  .  .OR.  R  I<  I  >  .  IE  .  0  .  .OR.  RK2),GT,0.)  STOP  10  00 
HOANGS  =  NUMAi-G.EO." 

NO  I  Mag  =  PI'  2  i,  gE  ,  II2CUT 


CALCULATE  NUMBER  OF  TERMS  IN  MIE  SERIES  '  h  LEAST  UPPER  BOUND 
USING  EMPIRICAL  FORMULAE  FITTED  FOR  SIZE  PwPAMETEPS  UP  TO 
2  0, 0  0 0 


IF0BETA.LT.  4200.  >  NT  =  BETA  M  05  +  BFTa  *  ".  1 
IFC  BETA  .GE  .  4200  .  )  NT  =  BE  Th  M  •*  BET  A  *  ♦ 1  1  .  , 
NTP  1  =  NTH 


COMPUTE  BIGh 


IF 

IC 

.- 1 

ARRAY  DIGA  HAS  ALREADY  BEEN  GEN 

IF  '■  I 
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■  GO  TO  2  i"i  Ci 
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* 

1  4 

P  F 
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-!C  C  ON  S  .‘STEM  fT'-.C  FOR  aRRA  6  1 

C  A  L  L 
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C  •  4  . 

i  T  Phi  Is  (hi  ,  I  H  I S  L  ,  I  E 1  T  > 

ICNWD 
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ILFN 
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Tfib  1 
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2INV'  1  .)  =  TEt-K  I  VBETn 
2IHVC2.')  =  TEMT2VBETH 

PREPARE  FOR  DOOM-RECURRENCE - 

COMPUTE  INITIAL  H I GH-ORDER  BICA-M  USING  LENTZ  METHOD 

15  FF<  I  >  =  NTPt  -2INVT 1  > 

F F <  2  >  =  NTPt  *ZINVC  2  > 

MM  -  -1 
KK  =  2* NT  +  3 

AKc  1  1  =  T  MMt-KK  .>*2  IN  VC  1  .> 

AK<  2  >  =  <  MM-rKK  :»*21NVT  2  .< 

BEN1:'  1  >  =  A  K  <  1  > 

DENT  2  -  AK<  2  ) 

CALL  DC  DVT  (i,  IDO..  0,  .  FT,  TEN  * 

NUNC  t  >  =  DENT  1  ')  +  TEH'  1  > 

NUNC  2  >  -  DENT  2  >  +  T  t’i'K  2  ’> 

KOUNT  =  I 


2  0  kviUNT  =  KOUNT  + 1 

IF< KOUNT , GT , MAX  1 T  >  GO  TO  40 

CALL  DC  DVT  -  1  ,  NUNC  1  > ,  NUMC  2  '• ,  hII  ,  TEM  > 

CALL  DCDV<  -  1  ,  DENT  1  » .  DEN*.  2  >  ,.iK  ,  TMa  > 

IF  <  TEM<  1  >  ,  GT  .  EPS  I  ,  and.  TlTAt  1  i.GT.EPS  1  >  GO  TO  30 

ILL-CONDITIONED  CASE --STRIDE  TWO  TERMS  INSTEAD  OF  ONE 
MM  =  -MM 

KK  =  KK+2 

AKT  1  >  =  «■'  MM  *KK  '»  *Z  I  NV>.  I  ') 

AK<  2  :>  =  <  MM+K.K  >*2I NV<  2  > 

NTH',  1  >  =  AN,  i  >  ♦: M U N '  1  >  -  mK<  2  n"NUMC  2  '  +  IDO 

NTNT  2  -  -  AK*'  I  >*NUM<.  2  >  +  hK>  2  >*MUMT  I  > 

DTDT  1  >  =  Ml  t  >*DEN<  1  >  -  hK(  2  DENT  2  >  ♦  IDO 

DTD<  2  T  =  A*  .  1  >  t'DEiK  2  >  +  AK  >  2  INDENT  !  > 

CALL  DCDVT  U  .  NTN1  I  > .  NTH'  2  >  ,  DTD  .  TEM  > 

TMA>:  t  >  =  TEM'  1  .  >i-F.;  1  -  -  TEN'.  2  «>  FFT  2  .) 

FFT  2  >  =  f  EM'.  1  » »'•  r-  •  2  >  +  TEM<  2  >‘i  FT  1  ) 

FFT 1  >  =  TMAC t  i 
MM  =  -MM 

KK  =  KKt-2 

AK'<  !,)-■:  MMH  V  '■  I  Z I M V 1  t  .* 

AK-:'  2  )  -  •.  MiNKK  '  *-21 HV*  2  ■ 

C  ALL  DC  DVT  0  HI  I  IT  I  ‘  I IUM'  2  .» ,  IITN  .  TEM  .> 

NUIT'  t  =  Ml  '  1  >  *•  TEM'  1  ) 

MUM1'  2  •  =  h>  ■  2  ’  *■  TEM'  2  • 

CALL  DCDVX  0  .  DEN'-.  I  > ,  IDLCM-  2  > .  DTD,  TEM  > 

Df  N'  I  i  -  Hi,':  I  ■  +  TEM  i  > 

t'L  N<  2  ’  =  Hi,  2  >  +  n~M'  2  * 

KOUNT  -  KOUNT+t 
GO  T0  20 

3  0  CALL  DC  DVT  0  MUM1.  I  '  .  Hi  :!•!■'  2  >  .  DEI-1  TT  * 

TEN'.  t.:»  =  TT’  j  )  *  F  F  •.  1  .»  -  TT*  i  >*FK\2  > 

FF'  2  '  =  TT'  !  1 +  FF-  2  +  TT'  2  1  *FF<:  1  > 

FFT  I  '  -  TFi:  I  ! 

'  'if  c  I '  f-  0  ‘  CONVLF  Ci  f t 

IF',  DhD  T  ID  II  •  ,  L  I :  .  EFC2  .  hND  .  DhC;  T  r  ’  2  .*  1  l.  T  EPS  2  <  GO  TO  SO 

MM  =  -MM 
K  K  =  I  K  4  2 

rif  \  I  »  -  1  ITM  ill  i ».  Z  I  tr.  >  !  I 
M K 1  2  >  -  ■  MM . ;  (  > I z 1 N -  ■  2  • 
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C  hLL  C> C C> V ■:  0  1  [>  0  ,  0  .  ,  HUM  ,  TEH  > 

HUM'  1  >  =  nt  •  1  >  ►  TEM'  1 

MUM1'  2  '  =  mi  <  E  ■  +  TEM'ii  .' 

CULL  C>  C  C>  V  0,  IDO,  0.  ,  DEM,  TEH) 

DEN>.  1  :>  =  nt  ■,  1  >  +  TEI-K  1  .> 

DEM<  2  >  =  At<  2  >  +  TEM1'  2  > 

GO  TO  20 

40  UPITE'  LUOUT , P.O01  1  HT,  BETA , PI , AK , HUM , DEH, TT , FF 

OOI  FORMA  T  v  /,  V  com  INUED  FF.t-iC  T I  Or!  FOR  M-iUB-HT  BAILED  fO  COHVERCC 
.*  HT=  ,  I.;/'  •  »  '  ,E20.8/  PEER  IHD£Ii  =  *  ,  2E20 , 3/  MK» * , 2E20 . 8/ 

*  HUM-  ,  r  Et'O  ■  8.  ’  '  DEH-  ,  2E2  0  .  :?/  '  TT-'.2E20  8/  FF  =  '  ,  2E2  0  .  ft  ) 

STOP  1 002 

5  0  RBIGh  =  FFv  1  .) 

IB  IGA  =  FFi  2 
BUFF-:  31  :>  =  PBIGh 
BUFF'  32.)  =  IBIGh 
MB  ■-  31 

DOWMWwPD  RECURPEMC'E  for  BICa-H 

DO  7  0  H  K  HT  ,2,-1 
TEM<  1  1  =  M » ZINV-:  1  ’  +  PBIGh 
i EM*.  2  >  -  H  »ZIHV'  2.)  ♦  IBIGA 
CALL  DCDV-:  0.  1D0, 0.  ,TEM,TMA  > 

RBIGm  -  H  +  ZIHm  !  >  -  r Mh>  1  > 

IBIGH  =  rmZIMV'  2  >  -  TMh'2> 

MB  =  MB -2 
BUFF':  MB  )  =  PPIGm 
BUFFCMB+1  )  =  IBIGh 
IF  >:  MB  .  HE  .  1  '  GO  TO  7  0 

CALL.  Z‘  'EC  ■  2  ICMI.ID,  BUFF  ,  JLEH,  IOF'1  ,  I0P2  > 

IF  1  HE!  O  .  2  .'  GO  TO  ?  O 
MB  -  7  7 


i  Ml  '  =  l  OP  1  -  1 

I  OP 2  [SECT -2 
7  0  CONTI  DUE 

Call  e: ;(•:<:*  2, 1 r m ( ; c  .buff,  ill*  iofi  ,  iopz  > 

SO  MB"  *  MO 

I  OP  10  •=  I  CPI 

i  np'  1  =  i  op.: 


;0o  Cum  'HUE 
HO  '  MB  O 
I  DPI  -  TOPIC 
I  OP 2  a  l  IIP 2  0 

C  wLL  EXEC  '1,1  0  HOC' ,  LUFF  .  ILEH,  I  CP  1  0  ,  IOP2  0  • 

PI  OP  IV  -  ID'!. 'PI'  l  ' 

inn  iwl.  T  ~  i."  ihOT  I  T  !  £':■  I  D  FOP  EFFfCIr.HT  C  hLCULaT  I  OH  OF 
MMIFRICmL  CO"  FF  I  C  I  EH  T'O  IN  Hit  'SEprEO- 
FN  =  *  C •  0 
PH  -  ID" 

MM  ~  1 

I H I  T  I  uL  t  OTi  P  !  i.  n  l  T  I -BESSEL  FUHCTIOH  .'ETh  r.ip  UPWARD  PECUFPEUC 
pc  ih  =  u  ;  i ••  lL.tm  ' 

CHIU  =  -  ITT,  > 


r-riNF-i  ■■■  psih.-bf.th  -  emu 

CHI  HP!  =  CHIH/Bt  Th  «■  PS  IN 
P  Z  E  T  r  <  -  P  1 1 1 
I ZCTN  3  CHl.-J 
PZTNP)  ---  PS  I  HP  1 
IZTNPt  =  CHI  HP! 

INITIALIZE  PP£V I UUS  CuC'FIC IEMT3  -  m-SUC-N- 
FOR  USE  It)  mS  .TIME  TPS  FhCI  OP  SERIES- 
hNI-I  !  =  'C'.O  i.i.u.i 

BNH I  =  1  0 . 0 ,  u .  0  1 
QEXT  =  0 . 

IN  IT  1  Al,  t.TE  F’li  PJM  hUD  SUMS-  fu-  + 


C  ALL 

D7M0V' 

o . 

D  o . 

n  . 

r"'  P  ,  ! 

.  MM  > 

CALL. 

DVHOV'. 

o , 

[-•ii  . 

0  / 

I  SP,  ! 

/MM  » 

CALL. 

f'VM  O'-1' 

0  , 

,  D  0 . 

h  , 

PSM.  1 

,  MM  • 

CALL 

[>  Vi'lO  7’ 

-  [•■  0 , 

0  , 

I  ■;  m  .  i 

,  mm  > 

CALL 

ir/MOV' 

n 

[>  0 

M  , 

.  PSi  2- , 

)  MM 

On.L 

DvriOV' 

i'i 

.  [I'I 

o , 

.  *  '  >  '  . 

)  MM 

CALL 

[>7!tC'.V< 

;  m 

DO 

■  1 

P- -1-1 - 

1  MM 

CALL 

DVKC.'V-: 

i. 

.  DO  . 

H  J 

,  i 

!  *‘M 

CALL 

DVMC7- 

.  D  O  . 

'  0 , 

,  C  I  Ml-!  ! 

,  !  M , 

call 

07  MOV' 

1 

.  D  O  . 

1.1  , 

,  P  I  1  !  ,  1 

.  MM  • 

WRITE1' LUOUT  ..20'  S  ■  ill'  ,1'OUHT  .  1 1.  ,  tlMf-lntlC 

forma r>  ,e  i z .  s  ' 


DO  5  0  0  )i  =  I  ,  NT 
PE'- 1  Oh  =  EUR  I'  ms 
I B I  On  =  run--'  Mb'  I  ' 

COMPUTE  T-'F  V -.PIOUS  N  imEPKhL  COF.FF  U  IFNI 

run  "  r-in 

i  ii.ii.'! !  (•••  i  -  (-  ti  ■ up  » 

;  (IP*,  -  !  ['I..--  m 

i.  r r  -  P  D  t  ;  i  ; 
nr  1  OH  =  i  D  •'  *i:  M 


ChLCUI.hTF  THE  Mil-  SERIES  C  OEFF  I C  I  ENTS  lITT 


IF',  MOIMhC  ■  GO 

TO  “  1 1  0 

GENEPhL  :  h 

:'-E 

Call  r.--c i.i  .  r 

d,h.  imo-c 

i  t  r  .  ■ 

TEH*  1  -  -  TCt"-' 

1  )  i-  SM.'lETm 

TNAC  1  :>  =  TEN- 

i  - 

n  !■  ;■  >  ip 

ITiA-  2  >  -  TEA' 

1  ii  I  Z  TUP  1  *- 

T  E'!'  2  '-PC 

.  TE'.-i 

i  •*!■;- imp  ’ 

'  i  run 

1  TMhv  i  ■ 

It'-  t-  1  I-'U' 

H.- 

Mil.  =  •  ’  TED 

|2'IU.  I  N't  '* 

n-i'-i  I  >  - 

',  Ti  lw  1  • 

T  V  t  '  V-i  • 

111  : 

I E M '  '■  '  =  Pin 

'irr'DH  -  Pi 

■  .j  .-  v  I  R  i  CiH 

TEI'U  2  ‘  =  PI'  t 

• «-  i  i  C.  h  i  P  I 

'  ■•:[■'!  c.h 

T|-!m'  i  -  -  i'Ct-1- 

I  ■  i  S 1  T : :  i-  !  - 

I  L  '■  ’  .  ‘  i  C 

Tf-tw  2  .  =  Tf 

;  '  * 

TC-  1  r  " 

LU'  1 


i  r '  1  -i  It.  -  •  i  il  l',  i  '  < 


J  I  P  ;.  i 


DC  Th 


1,  B-SOB-N- I  .') 

.  AT  uL.L  hHGI.  E  :■ 


NEEDED 


lE-h  hHD  LITTlC-B 

t : 

u 

PS- 1 tiP  »  ••  *  iMw  2  -  ' 
I  ir  .  -PS  I  tl  ■•♦-TI-Ih'  2  1  > 


r  i-  '  '  . 


ii 

iF-'lllr!  2 


300  CONTINUE 

NQ-mES 0 R P T  I  0 H  C h E 

TECK  1  .)  =  ■  RIOPI  /  ♦■PF: !  i,h  <  I'll. 'BETA  i*RZTNPi  -PZETM 
TEN''  2  '  =  R 1  OP  1 V « i  Gh  f  FM/BL  I  H  LMCrTtlP  1  - IZtTN 
TMA<  I  i  a  \  R  1  OP  I  V  t  RL<  I  Gh-*  Ftl/EE  T A  >*PS  I NP  1  -PS  I N 
ChLL  dcc»v  ■;  o .  rt n :■  t  ■  u ,  ,  tch  mi  ■ 

TEM1  1  >  a  •>  rj'.  1  >  *  RL»  !Sj»i*FI-r‘8tf  A  >-*RZTNP  1  -PZETN 

TEI'K  2  >  *  v  PI'..  1  '»+P.BI  Gh+FM/BETA  )*  I  ZTNP  1  -  I  ZE  TN 

TMAC  1  >  ~  •  P 1  •  1  >*RBIC,m*-FN/BETA  >*PSINP1  -PC  IN 

CALL  DC  0’v\  U  ,  T i  Ih1  1  i,0,  ,  TEH  ,  BII  < 

350  OEXT  a  01:.  r ■  t-TMi.irlP'  l  CNGL  Hlk  1  1  '♦SNGL'  Bit  t  >  >  » 

I F <  HOANG'S  •’  00 Tu  40  0 

PUT  CUE  SEP IEO  C OEFF I C I  CUTS  IN  FORM  NEEDED  FOR  COMPUTING  S+, 

PhNP  =  CORFU  hH'  1  '  t-BN'  1  >» 

I  MNP  s  C'TErr  ■*/  hi.''  L'  > ♦EN>.  2  >  > 

RC  HP  *  COE  FT  +•  Hit  1  >~PIK  I  i  < 

I3NP  =  COEFFf  hi 0  2  >•311'  2  •  > 

RANPM  -  MLHRhNP 
ImMPM  =  MM+IhIIP 
RBHPM  *  MI-H-rplIP 
I  BN  Pit  =  Mi'iflONP 

MOD  UP  SUMG  I'MUI  L  UPWARD  PEC  UPS  I MG  mN'JLmP  FUNCTIONS  LITTLE  PI 
MHD  LITTlE  ThU 

CALL  DVMPV*  ,;MU,  i  .PIN  1  .  VTEIT,  t  .  HN  > 

CALL  DVSUB> V  TEN .  I  . P I  MM t  ,  t  , THP ,  1  . UN  ) 

CALL  DVSM'f- FN,  TlIP,  I  .  T.-OJN,  1  MU  • 

CALL  DVSI.iD'.  TutIM  I  ,  p  I  Ml  t  .  !  TmUH,  1  ,  NN  • 

call  l'-vmo  /■  >■  i ii ,  i  . ri mm i  ,  i . hm  > 

ChI.L  DVR  I  V'  M  I  DN  ,  TiTP  t  VTEIT,  1  PIN  I  .  NN  ' 

ChLL  1 1  v ,  iC'l"  C  HIM  I  ,  I  ,  ThUH  .  !  .  /  !L,T  .  1  ,Ntl 

ChLL  DVPIV  <Mi '  .  7  TL,  i  t  P  p  1  ,  !  ' -P  ,  1  .  Nil  • 

CALI.  C'VP  I"'  I  hNP  ,  V  TEN  ,1  1  :  T-  .  t  ,  ISP  ,  1  .Nil.' 

C  ALL  OVP  I  V>.  C't  It'  IT  .  /  TEH  ,  I  ,  F  SMC  .  1  ,  l-CMS  ,  I  .  NN  > 

ChLL  DVP  IV'  IlMir  IT  »  VTC.iT ,  1  ,  l  SITS  ,  1  ,  I  CMS.  I  ,  HN  ' 

CALL  C' VC l  IR1  PIHHI  I  TmIJII,  I  , VTEM,  1  ,  NN  > 

call  dvi  r/'  ppm- ,  v n.iT ,  i  ,rvit,  t  ,pom,  i  ,nii  > 

Cm1,  i  uvr-tV'  M.oir  vimt.i  i;U.i  iom,  i  tin* 

ChLL  C'Vrr,'  Pn!!!-M  V  I  tC  i-l  .  t  .  RCC'C  ,  I  POP-.  ,  1  ,  IIH  • 

call  D'-t-  IV'  i,.NpM .  if  u  t  i .  ,  *  isps,  t  ,nh  - 

UPI'hTL  RElEVhHT  OUmNT  l  TIES  Pi.R  NF.KT  Ph'S'S  THFOUGH  loop 
45 U  ITM  ■=  -  MM 

FN  •=  FHPI 
RN  Pm- I 

HlllT  I  -  i  MCI  ;.||i, l  nl!'  I  ■  '  SHGL.i  hII>  c  i  ■  i 

BN  I'll  -  I  I  III,!  ■  Mil  I  I  I  VNO.I  '  IT  Ml.  ,?  »  I  ' 

|'|||.|T||  III.  :  i'.Til  CL',1.  nun  T  1 01.  ■  BY  ui  !'■  r.'C  EF.t.URREHCE 
rzlt  :  ■.  r  i 1 1  ,i  i ;  r  i .  ;  l-  ’  h  '  1 1.:  nip’  ~p;x  ru 
IZTT  -  ■  Tl.iOUr  1/1  ETh  -MhTNP! -izetn 
Rzr  hi  -  mi  "'  t 

I  Z L  III  -  ICTIIPI 
F  '  '  :•  I-  '  i 

I ;;  r  i  ■  i  •  i .  r 

PC  III  =  I  '  I III  I 
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PSINR1  *=  P2TNP1 
MB  =  MB  *2 

IF  <  ME-.LT.  33  >  GO  TO  500 
MB  =  1 

I 0P2  =  I0P2+2 

IF  < I0P2.LT. ISECT >  GO  TO  490 
I OP  I  =  IOP1+1 
I0P2  *  0 

490  CmLL  EXECC 1 , ICNMOj BUFF, I L  E  N , I OP  I , I0P2  > 
500  COM T I HUE 


QEXT  *  2  .  *i?Ei<T/BETH**2 
IFC  NOHNGS  )  RETURN 

RECOVER  SI  rtND  32  FROM  S+,  S- 

00  800  ,1  =  1  .MM 

SJ<  J  >  =>  0 . 3>*>CMPL»  3NGL  <.  PSP*;  J >+P3N<  J  >  >,SMGL<  ISP1,  J  >♦!$(*<  JO) 

$*<  ■i  i  *  0 , 5*CNPL/*.  $HGL‘  PSP'.  J  >-R$M<  0  )  .» ,  3NGL*  1  SP'  J  I  3MC  J  >  >  ) 

S  1  <  NMP  I  -  J  i*0 . 5*CMPL/C  3NGLC  R3PS-,  J  mRSNSC  J  >  >  ,  3NGLC  ISPS><  .)  >■*  I3MS<  J  ) 
800  32CNMP1-J  .>*  0,5  m"  MPLX\  3MGLC  PSPS'  J  > -ROMS',  J  >  > ,  '.NOLC  ISPS-:  J>-I3MS><  J  ) 
C 

RETURN 

END 

SUBROUTINE  DCDVCN.P1 , R2.S.Q) 

C  PERFORMS  COMPLEi;  mult  I  PL  I  CRT  ion  8,  01  VISION  IN  DOUBLE  PREC 

DOUBLE  PRECISION  R I  ,  R2  .  S-:  2  > ,  0-'  2  > ,  0 1 V 
IF  <  P  2  .  E  0 . 0  .  .PHD.  S'  2  > .  EO  0.  *  GO  TO  10 
DIV  =  SC  1  >t *2+3C 2 >  *  +  2 
IF  C  N  t  1,2 

1  C'C  1  >  ~  D30*T<  'Rif  ♦  2  *P2  *  *2  ».  '0 1  '• 

RETURN 

2  O'  1  ■>  »  CPI  +  S'  1  >  *r:  J.  ♦ S',  2  >  >.  0  I V 
i”1*'  2  *  *  C  P 2  *  i-1  1  i  ■  F  1  *S*  2  »  D  .  V 

RETURN 

10  OK  1  >  «  R 1  /S'.  1  ,i 
QC  2  >  =  0  . 

RETURN 

END 
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T  =  0  0 0  0 4  1C  OM  CP  0  0 023  *.(•*•  INC  00013  BLKS  P*<?0O0 


C 


’TH4  ♦>»C»W  2  GHT ,  23  JULV 

PROGEnM  CPU  w 

C 1MEM  -  1 11,1  PBUF1'  1  >,  I  C>CB<  144  ''H*TtE<.  - 
rLl  rt,HO,2P.I3  M3  1MTEK 
DOUBLE  PPEC!  3I0tl  N.'-r.THC 

SSS,1  1  ‘e«  *  - ;.««» > .  <  3  -■  t  -p 

„  ,  pgupv  4  » ,  -.«:*>  * .  -  peuf  ‘ .5  * .  -'2  3  > 

ri-C/O.’  iOPTrl.  m. '  .  1 L.  1 
LU  IN  =  LOGUJ’  I  DO  Mil  •’ 

WRITE-.  LU  IN,  . . *  ..  ,  T  HtinLE-, 

F,lRMu  l  •  "ENTER  R  2  .  t.t  •  h  ,  P  .  . --.RT,Lt4 

PEADvLUIN,  «•  '  PI  .BETh,PTS,::.TRT,  :TmNu 

WRITE*-  LUIN.iOvI  •  .....  ,  -  ,-H  :. 

FORMAT’S  3h:i.  12  > 

PI  s  ?,  1  41932 6 5 3 5 3 C> 0 
p  4  *  p  1  4  . 

$02  =»  EQRTt.  2.  ' 

h  m  I  ,  C'O/R  I  ,  . 

THC  3  PATrtM’-:  1  .  C*0.T'30RTi  H  +  +2-I  .Du.'  > 

TH2  ~  THC-U.oE  2 
TH~  «  THC-O.fE'2 
TH4  =»  THE-  I  .  ?E-3 
PHC  3  P  1  -  2  .  *-THC 
NS  -  M  Ml 

A  =  BE Trt/G  2  C»0»PI.» 
hK  s  $uPT'  h  "'03’  SN'jL’-  THC  1  1  1 
IH  a  GTR'ftTHC 
p  •=  ■  tmc-th  j/pt  :• 


2000 

2  0  0 1 
2  002 


E4 


EE 


b^OPEUeC’HTH  FILE  OEilMGATED  B. 

rnLL  OPEN’.  ICiCD  .  Ibl;P  .  MhIIE,  lot-  TH  .  ..-,L 

IF’  TEPP  ■  LT  .  0  ■>  GO  T.;i 
PBI.IF'.  I  >  *  BEth 
Ppi.iFO  .3  •’  -  ’  •  •'  M 

K<U.to'o:’  p.jm'.  r  ioh '  i  f : ' gig  ted  ;n  pbof  eeoiv, 

cat  L  MPtTf"  I  f » C  E*  tt.l;P  ,  PP'  IE  ,  I  L  ■ 

i  a  ';.;T  DO:  IIP  Cn-.LE*  „ 

•  '  r  T  ,,  .,  T  ■■■••■■  r-~-  I  [*«'•—  H  +  '"'T  >'*'*>£  '  •’ 

R  =1  t>H  T  *-01 '  H'lE.-  ■  • 

0  =  P-TH 
.?  .  K> 

_  i;.o$'  TH  • 

PC>  =  PH  IN  DEGREES 
r.|;(P7  : 

J  compute""  v»h"de  hhlgt  divergence  fhl  :op 

T . ' ■  2  f  •  !,.•*■"■  I  (’■  PH  '  ' 

=  p  *  TH 

TmM'  D  '.••TmM-  - 

■  --IN’  i'  ■  •  3  '  >»*2  , .  ttc1-.  , 


PH 

ET 

r 

pr> 


31 


FP 


0  EHCE  3T wTEMCNT 
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1 


J 

i 


1 


IP  =  G  1  M  .  ♦ '  I  ,-fp  M.  f  2 
I  -I-  =  Gl  *4  .  -FS  >4  +  2 
ETA  =  PHC-PH 
W  =  S  I  N<  ETt-t  /*Ak 
TH8  -  THC+ETA/2 

DELI  »  ATaNC  SORT*  S I  MC  THE  ':•*  *2  - 1  .  /NS  >/C  uS<  THB  >  ) 

DELS  =  2  .  +hThN<:  NS  >1 aNC  DELI  >1 
DEL  1  =  2  .  +DEL1 

ETAC  =  2  ,  ■+E-ETm-*<>  CO;-'  THB  >-C  T  +COS';  R  >/N  > 

B  I  =E  T nC  +  P4  K'EL  1 
B2»ETaC+P4+DFL2 
SI1-SORT*  I-*-  • 
s  1 2* SORT <  ir  :• 

CB  I  =  COS'  B 1  - 
C  62= COS'  B2  ■ 

'5  B  )  -  S  I M  C  &  1  .> 

SB2»SIN< B2  • 

GOTO  40 

THIS  CODE  USED  WHEN  PHSPHC 
30  ETAD  =  ETAD-CETh 
ETA  =  ETAD/5? , 2S53 
PH  =  PHC-ETh 
PD  =  PH+S7 , 2S5S 
W  *  SILK  ETA  .>  ♦  AK 
S  I  1  *  0  . 

SI  2=0, 

CB  1  =  1. 

CB2= 1  . 

set  =  o . 

SB2=  ii , 

4  0  call  freS'.  m,  fs,  n;  ■> 

FS=FS/S02 
FC-F  C/SQ2 

S- 1  P=FG  *C8  I  +  F S  ♦  S P  i  t  :  I  t 
S  2  R  -  FC  +  C  B2-*  F  S  *  S-  C , .  +  S  I  2 
S 1  I =  FS +C b 1 -FC f  Sb i 
S 2 1 =  F S  +  C B2-F  C +  SE  2 

COMPUTE  PHwSE  DIFFERENCE 
PHI  *5?  -  I'ASS+h  1  All-.  S!  I /SI ft  * 

PH2=57 . 2  SCO  ♦  hTmNi'  S-2 1  /S2P  > 

IPX  S1P..LT,  0.  )  FHI  =  1  so  ,  4  PH  I 
IF'  FHI  .  G  T  ,  too,  ,  FHI  =  PH  I  -'It-  0  . 

I  F<  S2P  .  L.  T  ,  0  ,  >  P H 2  =  I  8  0  -  +PH2 
IF<  FH2  ,GT  .  I  C  O  -  .*  PH2-=PH2--7aO  , 

DIFh-PH2-PH I 

IFF DIFA , GT ,  I  SO .  >  D  I Fw-0 I Fm-3o 0 . 

IF1!  DIFk.LT.  -  I  0  .  i  D1Fh  =  DIFm  *30  0  , 

COMPUTE  TOTAL  INTENSITIES 
T  I  S  = '  S  I  ft  *  *  2  +  S  1  I  ♦  *  2  > 

T  IP=i  S-2F  *  42  +  S2  1*42* 

I HTEH=ALOG  T  c  5*\ TIS  +  TIP >  > 

I  FC  MMOOC  PD  ,  I  E  -4  < .  EO  ,  n  ,  .  up  j  TEC  LU  IN  .  I  0  0  n  >  F  D  r  I S  ,  T I P .  D I  Fa 
0  FORMhT'  cS  ,4  2E  t  4 . 0  F  I  ft .  4  > 

CALL  UP I TFC  I  DC B ,  I  ERR , PBOi  ,  I L  > 

IFOH.LT  THC  >  GOTO  SO 
IFC  PD.LT.  S  r.iHG  >  GOTO  3..1 
GOTO  OS'? 

5  0  TFC  IH  ,  GE  ,  Til  7  ■  E  -  C  2 

if-  in  ,  ge  ,  ru  e  •  L - L"  3 

I  FC  I H  ,  GE  .  TH-I  '  E  F4 


tm-th+e 

IFOH.LT.THO  GOTO  10 
ETAD=0 , 

C'£TA=  ,  05 
GOTO  30 

900  CALL  CLOSE*:  JDCB  > 

WRITE*'  LU1N,  £003  >  IERR 
£003  FORMAT* ’'IERR3,\  14  > 

STOP 

END 

SUBROUT  It  it  FREG<  WE ,  FS ,  FC  ) 

THIS  SUBR  CALCULATES  THE  FRESNEL  INTEGRAL  FROM  0  TO  WE 
W  3  H6S*  WE  > 

F  a  (  |  ,  ♦  1 9£E*W  £  .  +  1  .  7‘3£^W+3 . 1  04  *W*  +•«*  > 

G  *  1  ,  A  £  ,  +4 , 1  42*W+3 . 492*W+*£+6  .  <»7  0*W**3‘> 

A  =  1 ,57079S*W**£ 

C  =  COS*A> 

S  *  S  2  N<  A  > 

FC  =  , 5+F^S-C^C 
FS  =  ,5-F*C-G*S 
IFCWE.LT.0.  >  GOTO  3 
FC  *  FC+ . 5 
F  S  3  F  3  +  ,5 
RETURN 

3  FC  3  ,5-FC 
FS  3  ,5-FS 
RETURN 
END 
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Paper  No.  7 


Glory  in  the  optical  backscatteriru)  from  air  bubbles  (D.  S.  Lanyley 
and  P.  L.  Marston)  Accepted  for  Publication  in  Physical  Review 
Letters 


Abstract 

Observations  of  light  backseat torod  from  air  bubbles  In  a  viscous 
liquid  demonstrate  an  enhancement  due  to  axial  focus  in,;.  A  physical-apt  i.es 
approximation  fo’*  the  cross-polarized  scattering,  correct  J  y  describes  the 
spacing  of  regular  features  observed.  The  non-cross-polarizod  scattering  is 
not  adequately  described  hv  a  single  class  of  rays. 


7  G 
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The  Mle  solution  for  electromagnetic  scattering  by  a  sphere  frequently 
does  not  lead  to  direct  interpretation  of  the  angular  scattering  pattern.  Con¬ 
sequently,  models  have  been  developed  to  facilitate  an  understanding  of  the 
structure  in  the  scattered  intensity  present  where  intensity  is  plotted  as  a 
function  of  the  scattering  angle  <j>  or  the  size  parameter  x  =  k(l  (k  =  wave- 

number;  a  =  sphere  radius).  These  models  have  emphasized  the  angular  regions 

2  3 

where  diffraction  is  important  for  a  drop  of  water  in  air:  the  rainbow,  ’ 

<fi  =  180°,^  ^  and  <p  ~  0°.^’^  In  the  scattering  of  light  by  a  spherical  air 

bubble  in  a  liquid  or  in  glass,  the  real  part  of  the  refractive  index  of  the 

sphere  is  less  than  that  of  the  surroundings  and  the  models  must  be  signifi- 

7  8 

cantly  modified.  New  phenomena  appear,  such  as  diffraction  ’  in  the 

region  of  the  critical  scattering  angle  ■]>  .  Here  we  report  the  first  detailed 

observations  of  backsc.attering  by  air  bubbles  in  liquids  and  give  a  model 

which  describes  some  of  the  observed  features.  We  refer  to  tills  as  glory 

3-5 

because,  as  in  the  case  of  drops,  the  i|>  180°  scattering  is  enhanced  when 

x  is  large. 

3  4 

Van  de  Hulst  ’  gave  a  partial  explanation  of  the  enhancement  for  crops 
by  noting  the  axial  focussing  of  those  backsca ttered  rays  which  have  a  non¬ 
zero  impact  parameter.  When  modeling  this  focussing  in  the  far  field,  dif¬ 
fraction  provides  an  essential  correction  to  ray  optics  because  the  facLor  in 
the  scattered  intensity  which  accounts  for  geometrical  divergence  of  the  rays 

goes  to  00  ns  <)>  -►  180°.  Examination  of  this  factor  in  ray-optics  models  of 

9 

scattering  by  bubbles  shows  that  this  «>  is  not  restricted  to  drops.  We  have 
modeled  the  backscattering  with  a  physical-optics  approximation.  The  proce¬ 
dure  is  to  (a)  compute  amplitudes  in  an  exit  plane  in  contact  with  the  bubble 
via  ray  optics,  and  (b)  allow  this  wave  to  diffract  to  the  far  field  where  the 
distance  from  the  bubble’s  center  R  >>  kci‘\ 
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Fig.  1  illustrates  several  rays  which  lead  to  backscatter ing .  The  paths 

are  determined  by  the  number  of  chords  p  and  m  ~  m./m  where  the  refractive 

i  o 

indices  of  the  inner  and  outer  media,  m.  and  m  ,  are  taken  to  be  real.  Fig.  .1 

i  ( i 

is  drawn  with  m  ^  =  1.403  which  corresponds  to  an  air  bubble  in  the  dimethyl- 
siloxanc-polymor  liquid  used  in  the  experiment.  All  rays  satisfy  sin0  *  msinp. 
For  4>  =  180°,  the  off-axis  (or  &lory)  rays  have  o  ~  t)  and  p  =  p where 0  = 
pp  +  (2g+2-p)90°,  g  is  a  non-negative  integer  (g  =  0  for  rays  in  Fig.  1)  and 
m  <  1  requires  that  p  >  3.  The  exit  plane  (dashed  line  in  Fig.  1)  touches  C' 
with  its  normal  parallel  to  the  propagation  direction  of  the  incident  wave. 

Our  description  oi  the  field  in  the  exit  plane  is  lacilitated  by  consid¬ 
ering  the  propagation  ol  a  wavelet  de  which  lies  close  to  the  hacksc.nttered 
path.  Fig,  1  shows  de  for  p  e  3;  it  emerges  as  curve  d'e'.  This  curve  appears 
to  come  from  a  ring, -like  source  at  F  known  as  the  local  circle  in  the  analo¬ 
gous  p  =  2  scattering  from  drops'  with  /?  m  •  2.  The  source  is  ring-like 
because  the  figure  may  be  rotated  around  t  lie  *  axis.  The  radius  ol  the  r  ing 
is  b  =  rtslnO.  After  the  incident  t ay  crosses  t  lie  dashed  vertical  plane  (the 
entrance  plane),  the  pro£aj»;»tJ.ori  phase  delay  for  teaching  the  exit  plane  is 
H  ■  ka[l  -  cost)  +  (1  -  cosB) sec  (0-P)  +  2mpcr>sp|.  The  ray  crosses  the  exit 
plane  at  a  radius  s  from  C'  witii  s/a  =■  slut'  -  (1  -  cosP)  t  an  (0-B)  .  The 
radius  n  of  arc  d'e*  follows  from  the  curvature  at  s  --  b:  rt.  =»  k(d  t|/ds*~)  " 

tX [  1  +  1  j ( P T  —  1 )  ^cosU)  where  i  -  tanp/tanO.  Tin-  spreading  of  the  wavelet  Is 


characterized  by  q  -  11m  d’er/de  as  de  >  0  where  tiie  bar  denotes  the  arc  length. 
An  equivalent  expression  for  q  is  |l.lm[h  -  s(0))/(h  -  asln())|  as  0  >  0; 

Its  value  from  [.'Hospital's  rule  is  a/  (a  --  a).  Vectors  (?.  *  1,2)  denote 

orthogonal  basis  vectors  in  both  the  cut  ranee  and  exit  planes;  Is  chosen 

parallel  to  the  polarization  of  t  he  incident  wave's  electric  field  E  exp(-iurt). 

n 

Tn  the  exit  plane,  the  i  ield  E'e^  of  t  ue  outgoing  ptli  glory  wave  is 
computed  by  applying  Van  de  Hu  1st’ a  method  of  first  decomposing  the  fields 
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perpendicular  and  parallel  to  the  scattering  plane.  ’  Exit-plane  polar 
coordinates  centered  on  C'  are  (s,i l>)  where  tj>  is  the  angle  relative  to 
and  s  and  4)  denote  local  basis  vectors.  We  assume  x  >>  1  and  use 
Fresnel's  coefficients  r^  for  the  internal  reflections  where  j  a  1,2  for 
fields  parallel  to  i|'  and  s,  respectively.  If  j  s— b  |  <<  CL,  the  multiple 
internal  reflections  give: 

F.J  -  Ejq"1'  F^exp [  in  +  ik(s  -  b)2/2u)  (1) 

~  ~  1  2  2  2 
where  n  ■  p  +  n(p  =  p),  F  (\p)  =  c^sin'ijj  +  c^cos  >p,  F  (i|>)  =  -  c^)sin2i^, 

and  Cj  ■  (-1)^^  “  r  j)  •  The  new  phase  term  p  accounts  for  the 

crossing  of  causLics  or  "focal  lines";  its  value  is^’^  -ir(p+g)/2.  The  r^ 

are  evaluated  at  8:  r^  *  sin(0  -  p)/sin(ti  +  p)  *  r^  “  tan(0  -  p)/tan(0  +  p) . 

The  sign  factor  in  Cj  accounts  for  a  geometrical  inversion  (present  when 

j  B  2  and  p  is  odd)  which  is  not  evident  in  descriptions  of  p  =  2  glory 

•  a  3,4 
in  drops. 

£ 

The  field  at  a  distant  point  Q  is  computed  as  follows.  The  left  exten¬ 
sion  of  the  CC1  axis  makes  an  angle  y  with  C’Q.  When  y  is  small  and  C,’Q  = 

2  12 

R'»ka  ,  scalar  diffraction  theory  and  the  Fraunhofer  approximation  give: 


MS,  ci(tR’  + 
2rr  iR '  q  J 


.,£  ik(s-b)  /2ct. 
nv  e  ds 


f£  e-iks  sinycos(l|j-f;)d(jj 


where  £  is  the  angle  between  e^  and  the  projection  of  C'Q  on  the  exit  plane. 
In  Eq.  (2),  the  approximation  given  by  Eq .  (1)  has  been  extended  beyond  its 
useful  domain  in  anticipation  of  the  stationary  phase  approximation  (SPA)  of 
the  integral.  Direct  evaluation  of  Eq.  (3)  gives  W  (Y,0  =  W  (y,f,,s  =  b) 
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“  ■n[(c1  +  c^)J^(u)  +  (CjL  -  c:  ^ )  J )  cos2fJ  and  W''  ~  it  (c^  -  c^JJ^Cu)  sin2£  where 
u  =  kbsiny.  The  SPA  oi  Eq .  (2)  Rives  the  pth  glory  contribution  to  the  scat- 
tered  field  when  kb  /u,  and  thus  x,  are  large.  In  the  experiments  to  be 
described  x  .>  A000  and  the  SPA  is  applicable. 

I 

The  total  field  may  be  approximated  hv  summing  the  K  from  Kq .  (2) 

with  the  fields  due  to  axial  reflections  and  surface  waves.  Surface  wave 

contributions  should  be  small  for  the  observed  bubbles  due  Lo  the  largeness 

of  x.  To  determine  which  glory  and  axial  terms  are  important  to  the  total 

field,  and  for  other  heuristic  reasons,  consider  the  x.- polarized  intensity 

1  of  the  pth  field  taken  alone.  The  SPA  ot  l.q .  (2)  gives: 

P 

I*  *  ( 2 / n ) x  1  p  f  I  WV'  ( ,  J.)\2  (A) 

2  2 

where  I  =  l  a  /AR  is  the  total  intensity  at.  a  distance  R  =  OQ  from  a 
K  I 

9 

perfectly  re  fleet  i  iijj  sphere  of  radios  rt  predicted  by  ray  optics,'  1  is 

2  s  2  1 

the  incident  intensity,  and  f  b "  u/rt  q  b  (o  -  i\)/d  .  in  Eq.  (A),  R 

P  *  K 

has  replaced  R'  from  (2)  and  y  becomes  180°  -  J;  because  R  (t . 

3,9  '9 

Geometrical  optics  gives  the  intensities  I  ot  separate  axial  reflections 

P 

•> 

(e.g.  p  =  0  and  2  in  rig.  1)  wh  i  -h  ate  propoi  t  i  ona  1  to  a".  The  strongest. 

reflection  has  p  ~  0  and  9.  1;  for  y  -  0,  I*  :  I  (in—  I  3  /(mil)*  wit  lie 

0  K 

“2  £  ( 

1_  “  0.  Since  f  does  not  depend  on  a,  I  "  ka  and  glory  terms  dominate 

0  P ,  g  p 

the  backseat t or i ng  when  a  is  large. 

Cons  idol  a  bubble  with  x  A00D  and  m  r-  L.AUl  The  strongest  glory 

terms  have  g  -  0  and  p  ---  3,  A,  and  3;  the  I  V I  for  y  -  0  are  respectively 

1.03,  0.A3,  and  0.16.  The  1  ^  decrease  with  increasing  p  tine  Lo  the  partial 

reflections  in  the  bubble.  The  strongest  axial  ray  gives  f^/l  =  0.028.  The 

(J  R 

interference  ot  the  fields  depends  on  a  and  our  Mio  computations  verify  that 


the  backseat t ered  intensity  is  not  simply  proportional  to  Cl  even  for  this 
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large,  value  of  x.  The  2  “  2  (cross-polarized)  scattering  is,  however, 

nearly  dominated  by  the  p  =  3  glory  term.  Due  to  symmetry,  1  -  2  scattering 

2 

vanishes  as  y  -*■  C.  The  I  (y  4  0,f.)  have  maxima  at  f,  =  f45°  and  4135°  and 

P 

they  vanish  at  f,  -  0°,  +90°,  and  180°.  Let  y  =  j  locate  the  first 

P 

2  2 
maxima  of  I^(y,C  “  45°).  The  largest  2  -  2  terms  have  I  (y  ,  (,  =  45°)/I^  = 

0.53  and  0.10  for  p  =  3  and  4.  To  the  extent  that  p  j  3  scattering  may  be 

neglected,  the  2.  ■  2  intensity  will  be  quasi-per indie  In  v. 

We  have  numerically  verified  the  validity  of  Eq .  (4)  by  using  Debye's 
3  4 

localization  principle  ’  to  modify  Mio  theory  so  that  only  partial,  waves 
associated  with  p  -  3  rays  were  included  in  the  Mie  series.  Furthermore,  when 
Eq.  (4)  is  applied  to  spheres  with  certain  m  >  1,  the  resulting  l^(y  =  0) 

p 

agree  with  the  glory  "analog"  tabulated  in  Ref.  11.  This  analog  was  derived 

by  applying  the  Watson  transformation  to  the  y  =  0  Mie  series. 

Fig.  2  diagrams  the  experiment.  A  syringe  injected  bubbles  into  the 

liquid.  The  liquid  had  a  high  kinematic  viscosity  (-600  000  cS)  and  a  single 

bubble  could  be  observed  for  hours.  The  laser's  power  output  was  5  mW  and  the 

beam  diam  was  5  mm.  The  wavelength  in  the  liquid  2n/k  was  632.8  nm/1.403; 

lay  in  the  splitter's  plane  of  incidence.  The  camera  was  iocused  on  so 

7  1  2 

the  photographs  recorded  the  far-field  intensity  pattern.  ’  Photographs 
were  made  with  a  -  0.3-0. 8  mm  corresponding  to  x  -  4000-11000.  Exposure  times 
were  typically  5s  for  TriX  film  and  a  200  nun  focal  length  camera  lens. 

Fig.  3  demonstrates  that  the  scattering  has  roughly  the  dependence  on 
£  predicted  by  Eq.  (4);  £,  =  0°  corresponds  to  scattering  toward  the  top  of 
the  photographs  and  y  =  0°  corresponds  to  the  center  of  the  symmetry.  Fig.  3(b) 
shows  that  the  9.  -  1  scattering  for  y  >  0.2°  is  significantly  stronger  for 
t,  =  ±90°  than  it  is  for  f,  =  0  .  This  agrees  with  the  following  model  resulLs: 
(i)  (Cj/c2)2  >>  1  (for  p  -  3  we.  predict  c^/c^  -5.2);  and  (ill  for  this  x. 
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vi  1 

the  1^  depend  only  weakly  on  £  and  are  dominated  by  the  I  .  One  prediction 

of  Eq.  (A)  could  be  quantitatively  checked:  when  both  siny  =  y  and  u  >>  1, 

2 

the  minima  in  I  should  be  spaced  by  Ay  rad  such  that  kbAy  =  ti  where  for 

p  *»  3,  b /a  =  0.AA7.  Fig.  A  compares  this  with  the  mean  spacing  of  ~A0  dark 

rings  lying  outside  the  9th  ring  from  the  center.  The  error  bars  combine 

uncertainties  in  measured  a  and  Ay  with  those  of  corrections  due  to  refraction 

at  the  cell-air  interface^  and  the  tilt  of  the  cell.  Fig.  A  shows  that  p  =  3 

rays  dominate  the  1  =  2  scattering.  The  modulations  of  the  intensity  along 

£  =  ±A5°  in  Fig.  3(b)  show  that  otlier  rays  contribute  to  SL  =  1  scattering 

1  2 

since  the  predicted  I  «  [J  (u))  . 

In  conclusion,  backseat tering  from  bubbles  can  be  enhanced  by  axial 

focusing.  The  number  of  significant  glory  terms  depends  on  m.  The  main 

3 

contributions  differ  from  those  for  water  drops  where  surface  waves  and 

other  diffraction  related  terms'*  play  an  essential  role.  If  focusing  were 

7-9 

not  present,  scattering  by  large  bubbles  would  be  r-<  1  in  the  region 

K 

(<J>c  +  10°)  <<{)<.  180°  where  t})^  =  2  cos  ^m  =  89°  for  m  ^  =  1.A03.  We  also 
find  evidence  of  p  =  3  glory  In  Mie  computations  for  bubbles  in  water . 

This  work  was  supported  by  the  Office  of  Naval  Research.  P.  L.  Marston 
is  an  Alfred  P.  Sloan  Research  Follow. 
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Fig.  Captions 

Fig.  1.  Rays  which  contribute  to  hackscattering.  The  local  angle  of  inci¬ 
dence  is  0  and  C  is  the  bubble's  center. 

Fig.  2.  Apparatus  for  observing  backscattering  from  bubbles. 

Fig.  3.  Photographs  for:  (a)  crossed  polarizer  (2  =  2  scattering); 

(b)  uncrossed  polarize  (2  ■  1);  and  (c)  no  polarizer.  The  incident  polar¬ 
ization  was  vertical,  a.  =  0.49  mm  and  x  =  6830. 

Fig.  4.  Measurement  and  model  for  the  angul.,r  separation  of  the  dark  rings 
in  the  2=2  scattering. 
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